Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,477 A→C 5.2% R98R (CGT→CGG yciF ← DUF892 domain‑containing protein YciF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,4770AC 5.2% 290.2 / 4.2 91R98R (CGT→CGGyciFDUF892 domain‑containing protein YciF
Reads supporting (aligned to +/- strand):  ref base A (49/37);  new base C (2/3);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTC  >  NC_000913/1315328‑1315626
                                                                                                                                                     |                                                                                                                                                     
aTGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCa                                                                                                                                                       <  1:541129/150‑1 (MQ=255)
aTGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCa                                                                                                                                                       <  1:392115/150‑1 (MQ=255)
 tGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCAc                                                                                                                                                      <  2:880136/150‑1 (MQ=255)
  gTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGCCTTTCTGTGCTGCGGCAATCAGTGCGGCATCACg                                                                                                                                                     >  2:696886/1‑150 (MQ=255)
   tCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGc                                                                                                                                                    >  2:125947/1‑150 (MQ=255)
      gTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACt                                                                                                                                                 >  2:756106/1‑150 (MQ=255)
        ggCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtga                                                                                                                                               >  1:405848/1‑148 (MQ=255)
         gCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgat                                                                                                                                              <  1:856327/150‑4 (MQ=255)
          ccTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGt                                                                                                                                             <  1:974831/150‑1 (MQ=255)
          ccTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGt                                                                                                                                             >  1:82626/1‑150 (MQ=255)
          ccTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGt                                                                                                                                             <  1:1405013/150‑1 (MQ=255)
           cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgattg                                                                                                                                            <  2:405848/150‑6 (MQ=255)
           cTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTGGtt                                                                                                                                            <  1:800551/150‑1 (MQ=255)
            ttttCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACAGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGTGATGCCGCACtgattgc                                                                                                                                           <  2:972091/150‑7 (MQ=255)
               tcttcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGTGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTtc                                                                                                                                        >  1:906419/1‑150 (MQ=255)
                  tcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGAAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCg                                                                                                                                     >  2:1395454/1‑150 (MQ=255)
                   cttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCgg                                                                                                                                    <  2:620342/150‑1 (MQ=255)
                   cttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCgg                                                                                                                                    <  2:369432/150‑1 (MQ=255)
                   cttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCgg                                                                                                                                    >  1:1173208/1‑150 (MQ=255)
                   cttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCgg                                                                                                                                    >  1:1168205/1‑150 (MQ=255)
                     tcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTa                                                                                                                                  >  2:239369/1‑150 (MQ=255)
                               ctttcAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGa                                                                                                                        >  1:640015/1‑150 (MQ=255)
                                 ttcAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGAcc                                                                                                                      >  1:425606/1‑150 (MQ=255)
                                  tcAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATATCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCt                                                                                                                     <  1:16233/150‑1 (MQ=255)
                                   cAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTc                                                                                                                    >  1:1158501/1‑150 (MQ=255)
                                    aGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCAATGACCTCa                                                                                                                   <  2:724636/150‑1 (MQ=255)
                                     gAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCAt                                                                                                                  >  2:1201004/1‑150 (MQ=255)
                                          ttCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGct                                                                                                             <  1:1261985/150‑1 (MQ=255)
                                          ttCGCTGCTTTACCGTAACCTAATTGTTCAGCCAGCGTCGCTAAAGTCCCGAAACTGGCAATCTCATAATGCTCGACAATCTCTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGct                                                                                                             >  2:1336023/1‑150 (MQ=255)
                                                gctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCa                                                                                                       <  2:554321/150‑1 (MQ=255)
                                                gctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCa                                                                                                       <  2:640015/150‑1 (MQ=255)
                                                 ctTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCaa                                                                                                      <  1:183594/150‑1 (MQ=255)
                                                    tACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGAt                                                                                                                          <  1:546499/127‑1 (MQ=255)
                                                    tACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGAt                                                                                                                          >  2:546499/1‑127 (MQ=255)
                                                     aCGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGAGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAg                                                                                                  >  1:775293/1‑150 (MQ=255)
                                                         tAACCTAATTGTTCAGCCAGTGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCt                                                                                              <  2:906419/150‑1 (MQ=255)
                                                            ccTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTcc                                                                                           >  2:1408238/1‑150 (MQ=255)
                                                              tAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCAt                                                                                         >  1:833748/1‑150 (MQ=255)
                                                              tAATTGTTCAGCCAGCGTCGCTAATGTCCCGAAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAg                                                                                                  >  2:902374/1‑141 (MQ=255)
                                                              tAATTGTTCAGCCAGCGTCGCTAATGTCCCGAAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAg                                                                                                  <  1:902374/141‑1 (MQ=255)
                                                                  tGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGcc                                                                                     <  1:239369/150‑1 (MQ=255)
                                                                  tGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGcc                                                                                     <  1:842447/150‑1 (MQ=255)
                                                                  tGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATATCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTAAATAAGACCTTCCATTGcc                                                                                     <  2:943740/150‑1 (MQ=255)
                                                                    ttCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcac                                                                                   <  1:1322402/150‑1 (MQ=255)
                                                                     tCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                  <  1:1408238/150‑1 (MQ=255)
                                                                                cGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTAATgcg                                                                       <  1:756106/150‑1 (MQ=255)
                                                                                  cTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCt                                                                     <  1:1409440/150‑1 (MQ=255)
                                                                                        tCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAttt                                                               >  1:840322/1‑150 (MQ=255)
                                                                                         cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                              >  1:46887/1‑150 (MQ=255)
                                                                                          ccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTc                                                             >  2:74559/1‑150 (MQ=255)
                                                                                           cGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCa                                                            >  2:897784/1‑150 (MQ=255)
                                                                                             tAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCa                                                            >  2:751251/1‑148 (MQ=255)
                                                                                             tAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCa                                                            <  1:751251/148‑1 (MQ=255)
                                                                                                  ggCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                     >  1:998884/1‑150 (MQ=255)
                                                                                                    cAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCCCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtg                                                   >  2:1163972/1‑149 (MQ=255)
                                                                                                        ctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGa                                               >  2:191460/1‑150 (MQ=255)
                                                                                                          cATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACt                                             >  2:1083692/1‑150 (MQ=255)
                                                                                                          cATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACt                                             <  2:425606/150‑1 (MQ=255)
                                                                                                            tAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTcc                                           >  2:311739/1‑150 (MQ=255)
                                                                                                             aaTGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCca                                          <  2:1168205/150‑1 (MQ=255)
                                                                                                             aaTGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCca                                          >  2:472056/1‑150 (MQ=255)
                                                                                                             aaTGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCca                                          <  2:1173208/150‑1 (MQ=255)
                                                                                                                gCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACaa                                       >  1:1185523/1‑150 (MQ=255)
                                                                                                                gCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACaa                                       >  2:310246/1‑150 (MQ=255)
                                                                                                                 cTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                      >  2:702730/1‑150 (MQ=255)
                                                                                                                 cTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                      >  2:497267/1‑150 (MQ=255)
                                                                                                                   cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACtt                                    <  2:766149/150‑1 (MQ=255)
                                                                                                                     aCTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGa                                  <  1:191460/150‑1 (MQ=255)
                                                                                                                       tttCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                <  1:74559/150‑1 (MQ=255)
                                                                                                                       tttCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                <  1:1395454/150‑1 (MQ=255)
                                                                                                                         tCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGCCTCCACAACTTGATCaa                              >  1:1122399/1‑150 (MQ=255)
                                                                                                                           tgtgCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATa                            <  2:82626/150‑1 (MQ=255)
                                                                                                                             tgctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACg                          <  1:125947/150‑1 (MQ=255)
                                                                                                                             tgctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACg                          >  2:758067/1‑150 (MQ=255)
                                                                                                                               ctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGtt                        >  2:201968/1‑150 (MQ=255)
                                                                                                                                tgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTc                       >  2:689946/1‑150 (MQ=255)
                                                                                                                                tgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTc                       <  1:497267/150‑1 (MQ=255)
                                                                                                                                tgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTc                       >  2:851698/1‑150 (MQ=255)
                                                                                                                                 gcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCa                      >  2:1130318/1‑150 (MQ=255)
                                                                                                                                   ggCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCTGAGTGCACAACTTGATCAATACGTTCAAt                    <  2:833748/150‑1 (MQ=255)
                                                                                                                                   ggCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAAt                    >  1:100773/1‑150 (MQ=255)
                                                                                                                                    gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                   >  2:292940/1‑150 (MQ=255)
                                                                                                                                    gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                   <  1:702730/150‑1 (MQ=255)
                                                                                                                                     cAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAg                                                           <  2:441761/109‑1 (MQ=255)
                                                                                                                                     cAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAg                                                           >  1:441761/1‑109 (MQ=255)
                                                                                                                                      aaTCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTg                 >  2:524448/1‑150 (MQ=255)
                                                                                                                                      aaTCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTg                 >  2:1358791/1‑150 (MQ=255)
                                                                                                                                       aTCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGt                >  1:1273749/1‑150 (MQ=255)
                                                                                                                                         cAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTcc              >  2:732226/1‑150 (MQ=255)
                                                                                                                                           gTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGCCTCCACAACATGATCAATACGTTCAATCTGTCCAt            >  1:1300829/1‑150 (MQ=255)
                                                                                                                                           gTGCGGCATCACGCACTTCGTTTTTCTAGGTACTTTCGATGAACTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTAATGCGCTTAATTTTCAGGTTAGATTCGGACTCAACAACTTGATCAATAAGTTCAATCTGTCAAt            >  1:1306016/1‑150 (MQ=255)
                                                                                                                                            tGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATg           <  1:1083692/150‑1 (MQ=255)
                                                                                                                                                    cACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCt   >  1:555253/1‑150 (MQ=255)
                                                                                                                                                     aCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTc  <  2:775293/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTC  >  NC_000913/1315328‑1315626

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