Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,023,525 G→T 7.0% intergenic (+13/‑29) tatC → / → tatD twin arginine protein translocation system ‑ TatC protein/3' ‑> 5' ssDNA/RNA exonuclease TatD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,023,5250GT 7.0% 563.8 / 16.6 170intergenic (+13/‑29)tatC/tatDtwin arginine protein translocation system ‑ TatC protein/3' ‑> 5' ssDNA/RNA exonuclease TatD
Reads supporting (aligned to +/- strand):  ref base G (74/82);  new base T (7/5);  total (83/87)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 177 total aligned reads displayed.

CGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTG  >  NC_000913/4023376‑4023674
                                                                                                                                                     |                                                                                                                                                     
cGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCg                                                                                                                                                       >  2:1681808/1‑150 (MQ=255)
cGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCg                                                                                                                                                       >  1:831847/1‑150 (MQ=255)
     aCGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAgg                                                                                                                                                  >  2:39328/1‑150 (MQ=255)
      cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAggg                                                                                                                                                 >  2:658080/1‑150 (MQ=255)
       gCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGc                                                                                                                                                <  1:2193257/150‑1 (MQ=255)
        cTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCg                                                                                                                                               <  2:1347351/150‑1 (MQ=255)
            tGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTg                                                                                                                                           <  1:120567/150‑1 (MQ=255)
              gCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTc                                                                                                                                         <  2:1890610/150‑1 (MQ=255)
               cGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTgaat                                                                                                                                        <  2:356081/150‑4 (MQ=255)
                gATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCTAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCat                                                                                                                                       <  1:1304666/150‑1 (MQ=255)
                 aTCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGCAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTgaattt                                                                                                                                      <  1:1554057/150‑6 (MQ=255)
                  tCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCatat                                                                                                                                     >  1:950240/1‑150 (MQ=255)
                       aTGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGt                                                                                                                                <  2:1523596/150‑1 (MQ=255)
                       aTGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGt                                                                                                                                <  2:1167043/150‑1 (MQ=255)
                         gTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGTAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTAc                                                                                                                              <  2:301811/150‑1 (MQ=255)
                          tACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTgaatt                                                                                                                                       <  1:51240/140‑5 (MQ=255)
                          tACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACa                                                                                                                             >  2:276071/1‑150 (MQ=255)
                          tACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACa                                                                                                                             >  2:879925/1‑150 (MQ=255)
                           aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            >  2:2148339/1‑150 (MQ=255)
                              gtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGAt                                                                                                                         >  2:1854162/1‑150 (MQ=255)
                              gtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGAt                                                                                                                         >  2:722061/1‑150 (MQ=255)
                              gtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGAt                                                                                                                         <  2:1715913/150‑1 (MQ=255)
                               tctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATg                                                                                                                        <  2:2250085/150‑1 (MQ=255)
                                 tgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGtt                                                                                                                      >  1:2155858/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  1:1591113/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATAAGGAGTACAGGATGttt                                                                                                                     >  2:2267365/1‑150 (MQ=255)
                                    ttGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGTTAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGATGAATTAATTCAACCGCGGGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGa                                                                                                                   <  1:2424434/150‑1 (MQ=255)
                                     tGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGat                                                                                                                  >  2:1691808/1‑150 (MQ=255)
                                        aaTCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGATGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATc                                                                                                               >  2:1981387/1‑150 (MQ=255)
                                          tCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCgg                                                                                                             <  1:2087504/150‑1 (MQ=255)
                                                tcttcttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTaa                                                                                                       <  2:1302793/150‑1 (MQ=255)
                                                tcttcttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTaa                                                                                                       >  2:2231616/1‑150 (MQ=255)
                                                   tcttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAAttt                                                                                                    >  2:1408350/1‑150 (MQ=255)
                                                    cttctCACGCTTTTATGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTg                                                                                                   >  2:1788975/1‑150 (MQ=255)
                                                    cttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAATACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTg                                                                                                   <  1:2486226/150‑1 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                <  2:1542078/150‑1 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                >  2:2379090/1‑150 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                >  1:1864786/1‑150 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGAGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                <  1:888305/150‑1 (MQ=255)
                                                        tcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCa                                                                                               >  1:2272161/1‑150 (MQ=255)
                                                         cACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAg                                                                                              >  1:266364/1‑150 (MQ=255)
                                                           cGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGtt                                                                                            >  2:144374/1‑150 (MQ=255)
                                                           cGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGtt                                                                                            >  1:1046580/1‑150 (MQ=255)
                                                           cGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGtt                                                                                            >  1:1523639/1‑150 (MQ=255)
                                                             cTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCg                                                                                          >  1:1481107/1‑150 (MQ=255)
                                                             cTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACAGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCg                                                                                          >  2:1832780/1‑150 (MQ=255)
                                                                  aCGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAAtt                                                                                     >  2:2180764/1‑150 (MQ=255)
                                                                   cGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAAttt                                                                                    >  1:1717034/1‑150 (MQ=255)
                                                                    gTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTg                                                                                   >  2:2417074/1‑150 (MQ=255)
                                                                       ggTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGcaa                                                                                <  1:619727/150‑3 (MQ=255)
                                                                       ggTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGa                                                                                <  2:19690/150‑1 (MQ=255)
                                                                         tAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGaaa                                                                              >  2:1491568/1‑150 (MQ=255)
                                                                                  aaaTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgat                                                                     <  1:725162/150‑1 (MQ=255)
                                                                                  aaaTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgat                                                                     <  1:2330966/150‑1 (MQ=255)
                                                                                    aTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgatga                                                                   <  2:1317369/150‑1 (MQ=255)
                                                                                    aTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgatga                                                                   >  2:2354600/1‑150 (MQ=255)
                                                                                         gAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGAtgttg                                                              >  1:1082905/1‑150 (MQ=255)
                                                                                         gAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGAtgttg                                                              <  1:234261/150‑1 (MQ=255)
                                                                                             agGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGc                                                          >  2:1873275/1‑150 (MQ=255)
                                                                                              ggAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCg                                                         <  1:2379090/150‑1 (MQ=255)
                                                                                                  aaCGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTgcg                                                     >  1:1924886/1‑150 (MQ=255)
                                                                                                    cgacgCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCt                                                   <  2:1647223/150‑1 (MQ=255)
                                                                                                      ccgcgGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCttt                                                 <  2:615084/149‑1 (MQ=255)
                                                                                                      acgCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCttt                                                 >  1:1324169/1‑150 (MQ=255)
                                                                                                      acgCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCttt                                                 <  1:1768367/150‑1 (MQ=255)
                                                                                                        gCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTg                                               <  1:1513107/150‑1 (MQ=255)
                                                                                                          tGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGAc                                             <  1:1691808/150‑1 (MQ=255)
                                                                                                               cAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTTATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCggg                                        <  1:2175698/150‑1 (MQ=255)
                                                                                                               cAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCggg                                        <  1:1698332/150‑1 (MQ=255)
                                                                                                               cAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCggg                                        <  2:2042576/150‑1 (MQ=255)
                                                                                                                 gAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAg                                      <  1:971000/150‑1 (MQ=255)
                                                                                                                 gAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAg                                      <  2:795115/150‑1 (MQ=255)
                                                                                                                 gAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAg                                      <  2:831847/150‑1 (MQ=255)
                                                                                                                  aaaGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGt                                     <  2:449230/150‑1 (MQ=255)
                                                                                                                   aaGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGtt                                    <  2:1864786/150‑1 (MQ=255)
                                                                                                                   aaGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGtt                                    >  2:924659/1‑150 (MQ=255)
                                                                                                                      gtaaaaaaCTGAATAATAAATTCAACCGCCCGTAAGTGCGGTTGTCATATGGAGTACTTTATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCTAAAGACCGTTATGATGTTGTAGCTTGCGCTTTTGACGCTGGAGTTAAt                                 <  2:974575/148‑1 (MQ=255)
                                                                                                                       gAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCTTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATg                                <  2:1530920/150‑1 (MQ=255)
                                                                                                                         aaaaaCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATggg                              >  1:1291065/1‑150 (MQ=255)
                                                                                                                         aaaaaCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATggg                              >  1:447526/1‑150 (MQ=255)
                                                                                                                           aaaCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCt                            >  1:438874/1‑150 (MQ=255)
                                                                                                                             aCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTAc                          <  2:2096199/150‑1 (MQ=255)
                                                                                                                             aCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGGCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTAc                          >  2:492327/1‑150 (MQ=255)
                                                                                                                               tGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtc                        <  2:645136/150‑1 (MQ=255)
                                                                                                                                gaagaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtca                       <  2:2289482/150‑1 (MQ=255)
                                                                                                                                gaagaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtca                       <  2:2344936/150‑1 (MQ=255)
                                                                                                                                gaagaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtca                       <  1:317020/150‑1 (MQ=255)
                                                                                                                                gaagaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtca                       <  1:2258268/150‑1 (MQ=255)
                                                                                                                                  agaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtcatc                     >  2:557979/1‑150 (MQ=255)
                                                                                                                                      tAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCg                 >  1:400979/1‑150 (MQ=255)
                                                                                                                                        aaTTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGTGCTACTCATCACCGGc               >  2:327718/1‑150 (MQ=255)
                                                                                                                                         tgacAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCa              >  1:845558/4‑150 (MQ=255)
                                                                                                                                         aTTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCa              <  1:69199/150‑1 (MQ=255)
                                                                                                                                          ttCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCAc             <  1:1676860/150‑1 (MQ=255)
                                                                                                                                            cAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTg                                               <  2:1700367/114‑1 (MQ=255)
                                                                                                                                            cAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTa           >  1:1347275/1‑150 (MQ=255)
                                                                                                                                              aCCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAAc         >  2:1570621/1‑150 (MQ=255)
                                                                                                                                               ccGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAAcc        <  1:978046/150‑1 (MQ=255)
                                                                                                                                               ccGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAAcc        <  1:1981387/150‑1 (MQ=255)
                                                                                                                                                    cGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGt   <  2:1046580/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTG  >  NC_000913/4023376‑4023674

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: