Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 468,271 (T)8→7 intergenic (+41/‑112) cof → / → ybaO thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913468,2640T.98.1% 464.9 / 2.0 104intergenic (+34/‑119)cof/ybaOthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (64/38);  total (66/38)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 109 total aligned reads displayed.

GCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGA  >  NC_000913/468168‑468359
                                                                                                |                                                                                               
gCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttt                                                                                              >  2:1518951/1‑100 (MQ=255)
   gTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGttttttt                                                                                           <  1:1353294/100‑1 (MQ=255)
       ctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttttttgcgt                                                                                       >  1:1928565/1‑96 (MQ=255)
       ctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttttttgcgt                                                                                       <  1:272154/100‑5 (MQ=255)
        tcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGtt                                                                                     >  1:424364/1‑100 (MQ=255)
         cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCTCTTCCAGCACCGGGCCATTGCCCGGttgtttgcgtt                                                                                       >  2:203127/1‑93 (MQ=255)
         cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTg                                                                                    <  2:1266439/100‑1 (MQ=255)
         cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTg                                                                                    <  2:3053477/100‑1 (MQ=255)
         cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTg                                                                                    >  1:1293013/1‑100 (MQ=255)
            tATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAAt                                                                                 <  1:2000770/100‑1 (MQ=255)
             aTTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAAtt                                                                                >  2:2946188/1‑100 (MQ=255)
               ttGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAAtt                                                                                <  2:1335458/98‑1 (MQ=255)
               ttGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTg                                                                              <  1:3010940/100‑1 (MQ=255)
               ttGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTg                                                                              >  2:2577292/1‑100 (MQ=255)
               ttGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTg                                                                              <  2:1873122/100‑1 (MQ=255)
               ttGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTg                                                                              <  1:628067/100‑1 (MQ=255)
                 gACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTc                                                                            >  2:2572231/1‑100 (MQ=255)
                 gACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑‑TTTTTTGCGTTGAATTTGTCa                                                                           >  1:1429930/1‑100 (MQ=255)
                   cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAt                                                                          >  2:365551/1‑100 (MQ=255)
                   cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAt                                                                          <  1:2527573/100‑1 (MQ=255)
                   cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAt                                                                          >  2:365541/1‑100 (MQ=255)
                   cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAt                                                                          >  2:1123216/1‑100 (MQ=255)
                   cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTACCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAt                                                                          >  2:574266/1‑100 (MQ=255)
                    gCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAtt                                                                         >  1:1740336/1‑100 (MQ=255)
                     cACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAttt                                                                        >  2:1709718/1‑100 (MQ=255)
                     cACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCAttt                                                                        >  1:2844384/1‑100 (MQ=255)
                       ctggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTtg                                                                      <  2:1274939/100‑1 (MQ=255)
                         ggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTtgtg                                                                    >  1:1121088/1‑100 (MQ=255)
                         ggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTtg                                                                      >  1:990097/1‑98 (MQ=255)
                         ggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAAATGCCCGGTTTTTTTGCGTTGAATTTGTCATTTtgtg                                                                    >  2:2154488/1‑100 (MQ=255)
                          gctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTtgtg                                                                    <  1:2982957/99‑1 (MQ=255)
                           ctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGTGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGcc                                                                  >  1:1268148/1‑100 (MQ=255)
                           ctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGcc                                                                  >  1:2358910/1‑100 (MQ=255)
                             ggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCgt                                                                >  1:1312445/1‑100 (MQ=255)
                             ggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCgt                                                                >  2:298317/1‑100 (MQ=255)
                              gACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtg                                                               <  2:513334/100‑1 (MQ=255)
                                 tATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtggtg                                                            >  2:2856085/1‑100 (MQ=255)
                                 tATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtggtg                                                            <  2:1079219/100‑1 (MQ=255)
                                  aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGgtgt                                                           >  2:1229719/1‑100 (MQ=255)
                                  aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGgtgt                                                           >  2:556982/1‑100 (MQ=255)
                                  aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGgtgt                                                           >  2:1661649/1‑100 (MQ=255)
                                   tCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGtt                                                          >  2:837871/1‑100 (MQ=255)
                                    ccACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGttt                                                         >  1:2918880/1‑100 (MQ=255)
                                     cacaTCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTa                                                        >  2:1358894/1‑100 (MQ=255)
                                     cacaTCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTa                                                        >  1:491877/1‑100 (MQ=255)
                                         tCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAAcc                                                    >  2:1721352/1‑100 (MQ=255)
                                          cTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTTCTTTTTTTTTTTCATTTTGCGGCGTGGTGTTTAAAc                                                     >  1:1847051/1‑98 (MQ=255)
                                          cTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAAcc                                                    <  2:838624/99‑1 (MQ=255)
                                           tACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGc                                                  >  1:2786349/1‑100 (MQ=255)
                                            aCCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGca                                                 >  2:2173144/1‑100 (MQ=255)
                                            aCCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGc                                                  >  1:814895/1‑99 (MQ=255)
                                             ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGcac                                                <  2:1282711/100‑1 (MQ=255)
                                             ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGcac                                                <  1:1721345/100‑1 (MQ=255)
                                             ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGca                                                 <  1:1769798/99‑1 (MQ=255)
                                              cTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGcaca                                               <  1:6810/100‑1 (MQ=255)
                                                tATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGa                                             <  1:1588048/100‑1 (MQ=255)
                                                 aTTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGaa                                            >  1:2397187/1‑100 (MQ=255)
                                                    cccccGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATaa                                         <  1:1494546/100‑1 (MQ=255)
                                                    cccccGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATaa                                         >  2:788395/1‑100 (MQ=255)
                                                    cccccGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATa                                          >  1:502978/1‑99 (MQ=255)
                                                      cccGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAAt                                       >  1:1996574/1‑100 (MQ=255)
                                                      cccGAATAACGAGATCCCTTCCAGCACCGGGCAATTACCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAAt                                       >  2:2163590/1‑100 (MQ=255)
                                                            tAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCTTTTAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGt                                 >  1:232840/1‑100 (MQ=255)
                                                             naCGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTg                                >  2:851922/2‑100 (MQ=255)
                                                             cgatctaTCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTg                                <  1:1270692/94‑1 (MQ=255)
                                                             aaCGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTg                                >  2:2798033/1‑100 (MQ=255)
                                                             aaCGAGATCCCTTCCAGCACCGGGCAAGTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAGCCGCACAGAATAAATTGTCGTg                                <  2:179295/100‑1 (MQ=255)
                                                              aCGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGa                               >  2:2629768/1‑100 (MQ=255)
                                                               cGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGAt                              >  2:386884/1‑100 (MQ=255)
                                                                  gaTCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTc                           >  2:1690623/1‑100 (MQ=255)
                                                                   aTCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGagatcg                          >  2:1270696/1‑95 (MQ=255)
                                                                    tCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAc                         >  2:2196562/1‑100 (MQ=255)
                                                                    tCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAc                         >  2:813783/1‑100 (MQ=255)
                                                                     cccTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAcc                        >  2:533099/1‑100 (MQ=255)
                                                                     cccTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTGTAAACCGCACAGAATAAATTGTCGTGATTTCAcc                        <  1:545813/100‑1 (MQ=255)
                                                                       cTTCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCtt                      >  1:336563/1‑100 (MQ=255)
                                                                        ttCCAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCttt                     >  1:419474/1‑100 (MQ=255)
                                                                           cAGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaa                  >  2:2972537/1‑100 (MQ=255)
                                                                           aaaaaCCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaa                  >  2:505867/5‑100 (MQ=255)
                                                                            aGCACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaa                 <  2:2918899/100‑1 (MQ=255)
                                                                              cACCGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaata               >  2:810876/1‑100 (MQ=255)
                                                                                ccGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaa             >  1:2882329/1‑100 (MQ=255)
                                                                                ccGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaa             <  1:164928/100‑1 (MQ=255)
                                                                                ccGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaa             <  2:2196034/100‑1 (MQ=255)
                                                                                ccGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTCTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaa             <  1:2154478/100‑1 (MQ=255)
                                                                                 cGGGCAATTGCCCGGTTTTTTTTTGTTTAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaaa            >  1:1914847/1‑100 (MQ=255)
                                                                                 cGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaaa            <  1:2242182/100‑1 (MQ=255)
                                                                                 cGGGCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaaa            >  1:1909708/1‑100 (MQ=255)
                                                                                    gCAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTa         >  2:2455541/1‑100 (MQ=255)
                                                                                     cAATTGCCCGGTTTTTTTGCTTTTAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCCCCTTTAAAATAAAATTa         >  1:409493/1‑99 (MQ=255)
                                                                                     cAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaa        >  1:1025993/1‑100 (MQ=255)
                                                                                     cAATTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaa        >  2:1522321/1‑100 (MQ=255)
                                                                                      aaTTGCCTGGTTGTGTTGCGTTGAATTTCTCATATTATCACTTGGTGTTTAAAACGCACAAAATAAATTGTCGTGATTTCACCTTTTAAATAAAAtt          <  1:1525342/97‑1 (MQ=255)
                                                                                       aTTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaaaa      <  2:30749/100‑1 (MQ=255)
                                                                                       aTTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaaaa      <  1:772423/100‑1 (MQ=255)
                                                                                       aTTGCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaaaa      <  1:1119651/100‑1 (MQ=255)
                                                                                          gCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAagag   >  2:362170/1‑100 (MQ=255)
                                                                                          gCCCGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAagag   >  1:1231291/1‑100 (MQ=255)
                                                                                           cccGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAgaga  <  1:2577278/100‑1 (MQ=255)
                                                                                           cccGGTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAgaga  >  2:2684030/1‑100 (MQ=255)
                                                                                                |                                                                                               
GCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGA  >  NC_000913/468168‑468359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: