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.99.3% 647.4 / ‑4.3 135intergenic (+34/‑119)cof/ybaOthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (57/77);  minor base G (1/0);  total (58/77)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CAACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTGACCGTAAG  >  NC_000913/468114‑468406
                                                                                                                                                      |                                                                                                                                              
cAACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGt                                                                                                                                                <  2:83768/151‑1 (MQ=255)
 aaCTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtt                                                                                                                                               >  1:1036772/1‑151 (MQ=255)
     gcgcgcGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttt                                                                                                                                             <  1:247471/149‑1 (MQ=255)
     gcgcgcGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttt                                                                                                                                             <  1:250555/149‑1 (MQ=255)
     gcgcgcGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttt                                                                                                                                             <  2:92309/149‑1 (MQ=255)
      cgcgggGAGCGCCCGCTTTTAGCGGTGAGTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGttttttt                                                                                                                                          >  1:494345/1‑151 (MQ=255)
         gcgGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttttttgcg                                                                                                                                       >  2:1641568/1‑148 (MQ=255)
         gcgGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGtttttttgcg                                                                                                                                       >  2:1532657/1‑148 (MQ=255)
             aGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTg                                                                                                                                   >  4:399021/1‑150 (MQ=255)
             aGCTCCCGCATTTACCGGGGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTg                                                                                                                                   >  1:1360735/1‑150 (MQ=255)
               cTCCCGCATTTACCGGTGATTGGACATTGCCGACATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCTCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAAt                                                                                                                                >  1:1241205/1‑151 (MQ=255)
               cTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAAt                                                                                                                                >  2:1735511/1‑151 (MQ=255)
                 cccGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAAttt                                                                                                                              >  2:1000236/1‑151 (MQ=255)
                 cccGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAAttt                                                                                                                              <  3:136500/151‑1 (MQ=255)
                  ccGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTg                                                                                                                             <  2:421797/151‑1 (MQ=255)
                     cATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTc                                                                                                                           >  2:449483/1‑150 (MQ=255)
                      aTTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCAt                                                                                                                         <  1:629673/151‑1 (MQ=255)
                        ttACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCAtt                                                                                                                        <  1:811024/150‑1 (MQ=255)
                        ttACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCAtt                                                                                                                        <  2:354638/150‑1 (MQ=255)
                         tACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCAtttt                                                                                                                      >  2:354745/1‑151 (MQ=255)
                          aCCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTtg                                                                                                                     <  2:1596354/151‑1 (MQ=255)
                            cGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTtgtg                                                                                                                   <  3:39972/151‑1 (MQ=255)
                              gTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTCCCGCCACC‑GGCCATTTCCCGGGTTTTTTTCCGTTGATTTTG‑CATTTTTTGcc                                                                                                                 >  3:528441/1‑150 (MQ=255)
                               tGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑‑TTTTTTGCGTTGAATTTGTCATTTTGTGCCg                                                                                                                >  1:13244/1‑150 (MQ=255)
                                 aTTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtg                                                                                                              <  2:573101/151‑1 (MQ=255)
                                   tGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtgg                                                                                                             <  1:1660455/150‑1 (MQ=255)
                                    ggACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCgtggtg                                                                                                           <  2:1319474/151‑1 (MQ=255)
                                      aCATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGtt                                                                                                         <  4:433464/151‑1 (MQ=255)
                                           gCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAAc                                                                                                    >  3:562029/1‑151 (MQ=255)
                                                aaTCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGcac                                                                                               <  1:457528/151‑1 (MQ=255)
                                                     ggCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAAt                                                                                          <  2:582904/151‑1 (MQ=255)
                                                         gTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATaaa                                                                                       <  1:1702999/150‑1 (MQ=255)
                                                         gTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATaa                                                                                        <  1:982460/149‑1 (MQ=255)
                                                          tctctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAAtt                                                                                     >  1:1710653/1‑151 (MQ=255)
                                                            tctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGt                                                                                   >  2:448971/1‑151 (MQ=255)
                                                             ctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAAtta                                                                                    <  1:158043/149‑2 (MQ=255)
                                                             ctcACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAAtt                                                                                     <  2:1599199/148‑1 (MQ=255)
                                                               cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCg                                                                                 <  1:591493/150‑1 (MQ=255)
                                                               cACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGt                                                                                <  2:342637/151‑1 (MQ=255)
                                                                 cTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTg                                                                               <  1:1685523/150‑1 (MQ=255)
                                                                   aTTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGa                                                                              <  2:21898/149‑1 (MQ=255)
                                                                   aTTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGAtt                                                                            <  2:305497/151‑1 (MQ=255)
                                                                    tttGACGCACTGGCTTGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGAGGTGTTTAAACCGCACAGAATAAATTGTCGTGAttt                                                                           >  1:330847/1‑151 (MQ=255)
                                                                       gACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCa                                                                         <  3:359620/150‑1 (MQ=255)
                                                                       gACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAc                                                                        <  1:730112/151‑1 (MQ=255)
                                                                        aCGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAcc                                                                       >  1:1631426/1‑151 (MQ=255)
                                                                         cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCt                                                                      <  2:1074749/151‑1 (MQ=255)
                                                                         cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCt                                                                      >  1:115997/1‑151 (MQ=255)
                                                                         cGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCt                                                                      >  2:1630701/1‑151 (MQ=255)
                                                                           cACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCtt                                                                     <  2:1299399/150‑1 (MQ=255)
                                                                             ctggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTta                                                                   >  1:1602182/1‑150 (MQ=255)
                                                                               ggctggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaa                                                                 >  4:486972/1‑150 (MQ=255)
                                                                                   ggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataaa                                                            <  2:789503/151‑1 (MQ=255)
                                                                                   ggACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataa                                                             <  2:696416/150‑1 (MQ=255)
                                                                                     aCTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTaaaataaaat                                                          <  2:1458244/151‑1 (MQ=255)
                                                                                       tATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAAtt                                                         >  3:556087/1‑150 (MQ=255)
                                                                                        aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaa                                                       <  1:1073841/151‑1 (MQ=255)
                                                                                        aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTa                                                        >  2:954765/1‑150 (MQ=255)
                                                                                        aTCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTa                                                        <  1:705673/150‑1 (MQ=255)
                                                                                           ccccatttccctaatCCCCCGGATTACGAAAATCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTaa                                                       <  4:528425/134‑1 (MQ=255)
                                                                                           cacaTCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAag                                                    <  2:1515248/151‑1 (MQ=255)
                                                                                            acaTCTACCTTATTACCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCTTTGAATTTGTCATTTTGAGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGCTTTCACCTTTAAAATAAAATTaaa                                                      >  1:809870/1‑148 (MQ=255)
                                                                                             caTCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAaga                                                   <  2:1064927/150‑1 (MQ=255)
                                                                                              aTCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAag                                                    >  2:1151970/1‑148 (MQ=255)
                                                                                                cTACCTTACTCCCCCGAATAACAAGAGCCCTTGCAGCACCGGGCAATTGCCCGG‑ATTTTTTGCGTTGAATTTGTCATATTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTGCACCTTTAAAATAAAATTAAAag                                                    <  1:338425/146‑1 (MQ=255)
                                                                                                  aCCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGaaaaa                                             <  1:1219539/151‑1 (MQ=255)
                                                                                                   ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGaaaaaa                                            <  1:519120/151‑1 (MQ=255)
                                                                                                   ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGaaaaaa                                            <  1:1532496/151‑1 (MQ=255)
                                                                                                   ccTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGATGTTTAAACCGCACAGAATAAATTGTCGTTATTTCACCTTTAAAATAAAATTAAAAGAGgaa                                               <  1:108584/148‑1 (MQ=255)
                                                                                                    cTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTtaaaataa                                                             >  4:145680/1‑133 (MQ=255)
                                                                                                      tATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGaaaaaaa                                           <  2:1595271/149‑1 (MQ=255)
                                                                                                      tATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGaaaaaaa                                           <  1:1630537/149‑1 (MQ=255)
                                                                                                      tATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAAtt                                         >  2:86732/1‑151 (MQ=255)
                                                                                                         tCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATtctc                                      >  1:118899/1‑151 (MQ=255)
                                                                                                         tCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATtctc                                      >  2:1283123/1‑151 (MQ=255)
                                                                                                             ccGAATAACGCTATCCCTACCAACACCGGGCAATTGCCCGG‑GTGTTTTGCGTTGAAGTTGTCATTTTGTGCCGTTGTGTTTAAACCGCACAGAATAAATTGTCGTGGTTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCtgtg                                  >  1:1353001/1‑151 (MQ=255)
                                                                                                                aaTAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGaa                               <  2:247663/151‑1 (MQ=255)
                                                                                                                   aaCGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAgg                             <  2:1729937/150‑1 (MQ=255)
                                                                                                                    aCGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGc                           <  1:1313496/151‑1 (MQ=255)
                                                                                                                     cGAGATCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCt                          >  2:167264/1‑151 (MQ=255)
                                                                                                                       agaTCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTAt                        >  1:6375/1‑151 (MQ=255)
                                                                                                                       agaTCCCTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTAt                        <  2:319318/151‑1 (MQ=255)
                                                                                                                            ccTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGtt                     >  1:1368025/1‑149 (MQ=255)
                                                                                                                            ccTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAg                   <  2:565354/151‑1 (MQ=255)
                                                                                                                             cTTCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGtt                     >  2:429644/1‑148 (MQ=255)
                                                                                                                               tCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGa                  <  2:250040/149‑1 (MQ=255)
                                                                                                                               tCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGAt                 >  2:549480/1‑150 (MQ=255)
                                                                                                                               tCCAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGAt                 >  2:111925/1‑150 (MQ=255)
                                                                                                                                ccAGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATa                <  2:1230920/150‑1 (MQ=255)
                                                                                                                                  aGCACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCCATGTTAGATaaa              >  2:1149957/1‑150 (MQ=255)
                                                                                                                                    cACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCTCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAAtt           <  2:1670586/151‑1 (MQ=255)
                                                                                                                                    cACCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAAtt           >  1:444339/1‑151 (MQ=255)
                                                                                                                                     aCCGGGCAATTGCCTGG‑TTTTTTTGCGTTGAATTTATCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCAACTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTg          <  1:226142/151‑1 (MQ=255)
                                                                                                                                     aCCGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTg          >  1:143193/1‑151 (MQ=255)
                                                                                                                                      ccGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTg          >  2:1453913/1‑150 (MQ=255)
                                                                                                                                      ccGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTg          >  1:477184/1‑150 (MQ=255)
                                                                                                                                      ccGGGCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTGa         <  1:1584413/151‑1 (MQ=255)
                                                                                                                                        gggCAATTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTGAcc       >  4:157313/1‑151 (MQ=255)
                                                                                                                                            aaTTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTGACCGTa    >  3:219243/1‑150 (MQ=255)
                                                                                                                                             aTTGCCCGG‑TTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAACTCTCTGTGGAAGGGCTATGTTAGATAAAATTGACCGTAAg  <  1:954644/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CAACTGCGCGCGGAGCTCCCGCATTTACCGGTGATTGGACATTGCCGAAATCAGGCTGTCTCTCACTATTTGACGCACTGGCTGGACTATCCACATCTACCTTATTCCCCCGAATAACGAGATCCCTTCCAGCACCGGGCAATTGCCCGGTTTTTTTTGCGTTGAATTTGTCATTTTGTGCCGTGGTGTTTAAACCGCACAGAATAAATTGTCGTGATTTCACCTTTAAAATAAAATTAAAAGAGAAAAAAATTCTCTGTGGAAGGGCTATGTTAGATAAAATTGACCGTAAG  >  NC_000913/468114‑468406

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