Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,357,593 C→A A357D (GCT→GAT)  SAUSA300_RS06725 → transketolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,357,5930CA98.9% 264.4 / NA 93A357D (GCT→GAT) SAUSA300_RS06725transketolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (90/2);  minor base G (1/0);  total (91/2)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAAGAATTTAAATTAGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGCTGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTAA  >  NC_007793/1357497‑1357688
                                                                                                |                                                                                               
tagaagaATTTAAATTAGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGtgatga                                                                                              >  1:755219/2‑100 (MQ=255)
           aaaTTAGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTg                                                                                   >  1:4867147/1‑100 (MQ=255)
               aaGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTAt                                                                               >  1:5126961/2‑100 (MQ=255)
                aGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTAtt                                                                              >  1:3882788/1‑100 (MQ=255)
                             aaaTTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTa                                                                 >  1:5484089/1‑100 (MQ=255)
                               aatGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTaaa                                                               >  1:2133658/3‑100 (MQ=255)
                                 tGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAgatcgga                                                             >  1:4069331/1‑94 (MQ=255)
                                 tGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAAc                                                             >  1:623663/1‑100 (MQ=255)
                                 tGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAAc                                                             >  1:5697293/1‑100 (MQ=255)
                                 tGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAAc                                                             >  1:3107653/1‑100 (MQ=255)
                                 tGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAACAACCAgatcgga                                                             >  1:2266356/1‑94 (MQ=255)
                                  gCCTAAAAATTATAAGGATGAATTACCACGTTTTGCACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACt                                                            >  1:1428354/1‑100 (MQ=255)
                                  gCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACt                                                            >  1:149570/1‑100 (MQ=255)
                                  gCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACt                                                            >  1:605986/1‑100 (MQ=255)
                                  gCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACt                                                            >  1:4993857/1‑100 (MQ=255)
                                      aaaaaTTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTcc                                                        >  1:1084143/1‑100 (MQ=255)
                                      aaaaaTTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTca                                                        <  1:805539/100‑2 (MQ=255)
                                                      aaTTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTtggtgg                                        <  1:1976416/100‑1 (MQ=255)
                                                        ttACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGAt                                      >  1:5424934/1‑100 (MQ=255)
                                                        ttACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGGGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTAGAt                                      >  1:1715425/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGaga                                     >  1:5263425/1‑98 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATc                                     >  1:4626623/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATc                                     >  1:2072307/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATc                                     >  1:1392325/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATc                                     >  1:1091484/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATc                                     >  1:2815480/1‑100 (MQ=255)
                                                         tACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTCGGTGGATc                                     >  1:2817184/1‑100 (MQ=255)
                                                          gccACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:4275863/2‑100 (MQ=255)
                                                          cccACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:3918690/2‑100 (MQ=255)
                                                          aCCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGTATca                                    >  1:5596432/1‑100 (MQ=255)
                                                          aCCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:5250578/1‑100 (MQ=255)
                                                          aCCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:1878169/1‑100 (MQ=255)
                                                          aCCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:461364/1‑100 (MQ=255)
                                                          aCCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATca                                    >  1:4224266/1‑100 (MQ=255)
                                                           gcACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATcag                                   >  1:2422511/2‑100 (MQ=255)
                                                           gcACGTTTTGAACTGGGTCACAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTGTGGTAGATcgg                                   >  1:2453685/2‑98 (MQ=255)
                                                            cccGTTTTGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATcagc                                  >  1:4455273/3‑100 (MQ=255)
                                                                 tttGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGAcc                             >  1:1052474/1‑100 (MQ=255)
                                                                  atGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCt                            >  1:2442334/2‑100 (MQ=255)
                                                                   tGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCGGCAGACCtt                           >  1:2466921/1‑100 (MQ=255)
                                                                   tGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCtt                           >  1:1534878/1‑100 (MQ=255)
                                                                   tGAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGAACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCtt                           >  1:820696/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGGTGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:830208/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:175427/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:1118075/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:1467436/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:665276/1‑100 (MQ=255)
                                                                    gAACTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGATATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:2754069/1‑100 (MQ=255)
                                                                    ctaCTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTg                          >  1:4483252/3‑100 (MQ=255)
                                                                     gaCTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGc                         >  1:2879522/2‑100 (MQ=255)
                                                                     caCTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGc                         >  1:3850201/2‑100 (MQ=255)
                                                                     aaCTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGc                         >  1:3858555/1‑100 (MQ=255)
                                                                      cttGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCt                        >  1:2294437/3‑100 (MQ=255)
                                                                       gtGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTg                       >  1:2197671/2‑100 (MQ=255)
                                                                       cTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTg                       >  1:1186049/1‑100 (MQ=255)
                                                                       cTGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTg                       >  1:1598557/1‑100 (MQ=255)
                                                                        tGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTgg                      >  1:5761213/1‑100 (MQ=255)
                                                                        tGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTgg                      >  1:5152789/1‑100 (MQ=255)
                                                                        tGGGTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGCAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTgg                      >  1:876151/1‑100 (MQ=255)
                                                                         gggTCTTAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTCTGGTAGATCGGAAGagcacacgtct                     >  1:3551503/1‑90 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:4285393/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:960746/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:3755743/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:1036549/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:2937119/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:1479390/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:1612252/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:2200209/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:3246001/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:2206062/1‑100 (MQ=255)
                                                                         gggTCATAATGGTGCATCTCGAGATGATTCTGGAACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGt                     >  1:1324546/1‑100 (MQ=255)
                                                                             cATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTcaa                 >  1:3795705/1‑100 (MQ=255)
                                                                             cATAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTcaa                 >  1:854446/1‑100 (MQ=255)
                                                                              ctAATGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTcaaa                >  1:274702/2‑100 (MQ=255)
                                                                                aaTGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTcaaaca              >  1:5765323/1‑100 (MQ=255)
                                                                                 ctGGTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTcaaacaa             >  1:1273122/2‑100 (MQ=255)
                                                                                   ggTGCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAAt           >  1:5012358/1‑100 (MQ=255)
                                                                                      gCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCa        >  1:2540246/1‑100 (MQ=255)
                                                                                      gCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCa        >  1:993193/1‑100 (MQ=255)
                                                                                      gCATCTCGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCa        >  1:3191833/1‑100 (MQ=255)
                                                                                          ctcGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGt    >  1:2556057/1‑100 (MQ=255)
                                                                                          ctcGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGt    >  1:5723199/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAACTCCAATGTaa  >  1:4514376/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTag  >  1:2692021/1‑99 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:579238/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:3285941/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:1952618/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:3363868/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:3415305/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:5000766/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:3593938/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:959277/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:4045030/1‑100 (MQ=255)
                                                                                            cGTGATGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTaa  >  1:4369397/1‑100 (MQ=255)
                                                                                                |                                                                                               
CAGAAGAATTTAAATTAGCGATTAGTGGTAAATTGCCTAAAAATTATAAGGATGAATTACCACGTTTTGAACTGGGTCATAATGGTGCATCTCGTGCTGATTCTGGTACTGTTATTCAAGCAATCAGTAAAACTGTCCCTTCATTCTTTGGTGGATCAGCAGACCTTGCTGGTTCAAACAAATCCAATGTAA  >  NC_007793/1357497‑1357688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑