Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,167,915 C→A 5.4% intergenic (‑101/+74) yegZ ← / ← yegR pseudogene, gpD phage P2‑like protein D/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,167,9150CA 5.4% 270.1 / 2.4 86intergenic (‑101/+74)yegZ/yegRpseudogene, gpD phage P2‑like protein D/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (42/39);  new base A (3/2);  total (45/41)
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

AAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACTATTTG  >  NC_000913/2167769‑2168063
                                                                                                                                                  |                                                                                                                                                    
aaaTATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGc                                                                                                                                                   <  2:61224/150‑1 (MQ=255)
   tatCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATc                                                                                                                                                >  2:662891/1‑150 (MQ=255)
     tCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGt                                                                                                                                              <  1:562561/150‑1 (MQ=255)
        gCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCt                                                                                                                                           <  1:1683900/150‑1 (MQ=255)
               cGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCCTAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAg                                                                                                                                    <  1:1233355/150‑1 (MQ=255)
               cGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAg                                                                                                                                    <  2:1507023/150‑1 (MQ=255)
                ggTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAgg                                                                                                                                   <  1:949774/150‑1 (MQ=255)
                ggTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTATTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAgg                                                                                                                                   >  2:681018/1‑150 (MQ=255)
                   tgCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTc                                                                                                                                <  1:1496529/149‑1 (MQ=255)
                   aGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTc                                                                                                                                >  2:1298996/1‑150 (MQ=255)
                     ccAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCg                                                                                                                              >  1:1591192/1‑150 (MQ=255)
                        ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCt                                                                                                                           <  2:1568663/150‑1 (MQ=255)
                               ggAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGcaca                                                                                                                    >  1:625368/1‑150 (MQ=255)
                               ggAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCATCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACCGGTTCCGGCTACTcacc                                                                                                                    >  1:216738/1‑149 (MQ=255)
                                gagaACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCGTCATTCATAGCACGATGCCTTTCACTTATCTGCATCGTGCTCACACAGGTTCCGTCTACGCACAg                                                                                                                   >  1:829842/1‑150 (MQ=255)
                                gagaACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAg                                                                                                                   <  1:1414930/150‑1 (MQ=255)
                                gagaACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAg                                                                                                                   <  1:475868/150‑1 (MQ=255)
                                       ctGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAAc                                                                                                            >  1:421171/1‑150 (MQ=255)
                                        tGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCTGGCTGTGCGTAGCCGGAAcc                                                                                                           <  2:147614/150‑2 (MQ=255)
                                         gCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAAcgc                                                                                                          >  2:1185421/1‑150 (MQ=255)
                                              acGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACTATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCTGGCTGTGCGTAGCCGGAAcctgtggg                                                                                                     <  2:1667815/150‑8 (MQ=255)
                                                  ggttATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTg                                                                                                 <  2:558908/148‑1 (MQ=255)
                                                   cTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGa                                                                                                >  1:478113/1‑150 (MQ=255)
                                                   cTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGa                                                                                                >  1:338602/1‑150 (MQ=255)
                                                    ttATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGAc                                                                                               <  1:1400755/150‑1 (MQ=255)
                                                      aTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCTGGCTCCGCACAGCCAGAACGCGCATATTTGACGc                                                                                             >  2:1659629/1‑150 (MQ=255)
                                                      aTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGc                                                                                             >  1:1204837/1‑150 (MQ=255)
                                                      aTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGc                                                                                             >  2:797635/1‑150 (MQ=255)
                                                           tttATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTTACCACATCGATAGAATCCTCCTTCATACCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATCTGACGCTTAcc                                                                                        <  1:356032/150‑1 (MQ=255)
                                                           tttATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCAGAGGTTCCGGCTACGCACAGCCAGAACGCGCATTTTTGACGCTTAcc                                                                                        <  1:1331935/150‑1 (MQ=255)
                                                           tttATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTAcc                                                                                        <  2:759579/150‑1 (MQ=255)
                                                           tttATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTAcc                                                                                        <  2:1112321/150‑1 (MQ=255)
                                                                 cATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGCACGCGCATATTTGACGCTTACCaaaaaa                                                                                  >  1:682637/1‑150 (MQ=255)
                                                                  aTAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAat                                                                                 >  1:1429221/1‑150 (MQ=255)
                                                                  aTAAGCCCTATCTCTGGTAACCGTCTTACATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCAACAGGTTCCGGCTACGCACACCCAGAACGCGCATATTTGACGCATACCAAAAAat                                                                                 >  1:27220/1‑150 (MQ=255)
                                                                    aaGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTCCGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAatat                                                                               <  2:210655/150‑1 (MQ=255)
                                                                     aGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGGATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAatatg                                                                              >  2:953547/1‑149 (MQ=255)
                                                                      gCCCTATCTCTGGTAAGCGTCTTCCATTGACCACATCGTTAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATAGTCTTCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATtc                                                                             <  1:269611/150‑1 (MQ=255)
                                                                      gCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATtc                                                                             >  2:168258/1‑150 (MQ=255)
                                                                        ccTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATtctc                                                                           <  2:338602/150‑1 (MQ=255)
                                                                            tctcTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAATATTCTCActct                                                                      >  1:107035/1‑150 (MQ=255)
                                                                             ctctGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCActct                                                                      >  1:1225393/1‑150 (MQ=255)
                                                                               ctcgTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACATATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCAAATTTGACGCTTACCAAAAAATATTCTCACTCTcc                                                                    >  1:1232366/4‑150 (MQ=255)
                                                                               ctGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTcc                                                                    <  1:862054/150‑1 (MQ=255)
                                                                                    aaCCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACAtt                                                               <  1:30118/150‑1 (MQ=255)
                                                                                      ccGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTg                                                             <  2:990136/150‑1 (MQ=255)
                                                                                        gTCTTCCATTGACCACATCGATAGAATCCTCCTTCATATCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGaa                                                           >  2:323161/1‑150 (MQ=255)
                                                                                        gTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGaa                                                           <  1:662891/150‑1 (MQ=255)
                                                                                         tCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAAt                                                          <  1:1419081/150‑1 (MQ=255)
                                                                                          cTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATg                                                         <  2:1721198/150‑1 (MQ=255)
                                                                                          cTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATg                                                         <  2:422988/150‑1 (MQ=255)
                                                                                           ttCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGt                                                        >  1:654377/1‑150 (MQ=255)
                                                                                            tCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTc                                                       <  2:625927/150‑1 (MQ=255)
                                                                                            tCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTc                                                       >  2:824/1‑150 (MQ=255)
                                                                                                  gACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGa                                                 >  1:508643/1‑150 (MQ=255)
                                                                                                   aCCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAg                                                >  1:229188/1‑150 (MQ=255)
                                                                                                     cacaTCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCg                                              >  2:323476/1‑150 (MQ=255)
                                                                                                          cGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACAcc                                         <  2:1441174/150‑1 (MQ=255)
                                                                                                          cGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACAcc                                         >  2:794190/1‑150 (MQ=255)
                                                                                                          cGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACAcc                                         >  2:274779/1‑150 (MQ=255)
                                                                                                               gAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTa                                    >  1:247821/1‑150 (MQ=255)
                                                                                                                 aTCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTaaa                                  >  2:1441802/1‑149 (MQ=255)
                                                                                                                 aTCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAAt                                  >  1:1590128/1‑150 (MQ=255)
                                                                                                                  tcctccTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGCCGCTTACCAAAAAATATTATCACTCTCCACATTTGAATGTCAGACGAGCGACCCCATGTAATc                                 >  1:1403319/1‑150 (MQ=255)
                                                                                                                   cctcctTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATcc                                <  2:1160778/150‑1 (MQ=255)
                                                                                                                     tcctTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTa                              <  1:824/150‑1 (MQ=255)
                                                                                                                        ttCATAGCACGATGCCTTTCACTTATCTGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTacac                           <  1:1228400/150‑1 (MQ=255)
                                                                                                                        ttCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTacac                           <  1:797635/150‑1 (MQ=255)
                                                                                                                          cATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCt                         >  2:1213634/1‑150 (MQ=255)
                                                                                                                          cATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTACGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACACTCCACATTTGAATGTCAGACGAGCGACACAATGTAAACCTACACCt                         >  2:254517/1‑150 (MQ=255)
                                                                                                                           aTAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCtt                        >  2:1507239/1‑150 (MQ=255)
                                                                                                                            gaGCACGATGCCGATAAGTGAAAGGCATCGTGCTCCCACAGGTTCCGGCTACGCACCGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTc                       >  1:479865/2‑150 (MQ=255)
                                                                                                                            gaGCACGATGCCGATAAGTGAAAGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACCCCTTc                       >  1:1652730/2‑150 (MQ=255)
                                                                                                                            gaGCACGATGCCGATAAGTGAAAGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCAGGTAATCCTACACCTTc                       >  2:1039747/2‑150 (MQ=255)
                                                                                                                              gCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCTCATATTTGACTCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTg                     <  1:581855/150‑1 (MQ=255)
                                                                                                                              gCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTg                     <  2:748692/150‑1 (MQ=255)
                                                                                                                                 cGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCt                  >  1:1353257/1‑150 (MQ=255)
                                                                                                                                    ttccTTTAACTTATAGGCATAGTGCTCCCACAGTTTAAGGATTCGTAGAGCCAGAACTAGAATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATATGAATGTCCGACGTGCGACACCATGTAATCCTACACCTTCTGTCTTCa               <  1:156485/148‑1 (MQ=255)
                                                                                                                                     gCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAg              <  2:507142/150‑1 (MQ=255)
                                                                                                                                      ccTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGc             >  2:922624/1‑150 (MQ=255)
                                                                                                                                      ccTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGc             <  1:471641/150‑1 (MQ=255)
                                                                                                                                           cACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAAc        >  2:1670878/1‑150 (MQ=255)
                                                                                                                                            aCTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACt       >  1:446335/1‑150 (MQ=255)
                                                                                                                                               tATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACTAtt    >  1:865890/1‑150 (MQ=255)
                                                                                                                                                 tCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACTATTTg  <  2:654377/150‑1 (MQ=255)
                                                                                                                                                 tCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACTATTTg  <  1:168258/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTCCTTCATAGCACGATGCCTTTCACTTATCGGCATCGTGCTCCCACAGGTTCCGGCTACGCACAGCCAGAACGCGCATATTTGACGCTTACCAAAAAATATTCTCACTCTCCACATTTGAATGTCAGACGAGCGACACCATGTAATCCTACACCTTCTGTCTTCAGCTCAACTATTTG  >  NC_000913/2167769‑2168063

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