Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,023,525 G→T 7.3% 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.3% 348.0 / 10.0 108intergenic (+13/‑29)tatC/tatDtwin arginine protein translocation system ‑ TatC protein/3' ‑> 5' ssDNA/RNA exonuclease TatD
Reads supporting (aligned to +/- strand):  ref base G (54/46);  new base T (5/3);  total (59/49)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 113 total aligned reads displayed.

CGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTG  >  NC_000913/4023376‑4023674
                                                                                                                                                     |                                                                                                                                                     
cGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCg                                                                                                                                                       <  2:604562/150‑1 (MQ=255)
cGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCg                                                                                                                                                       >  2:1354638/1‑150 (MQ=255)
cGCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCg                                                                                                                                                       >  2:1628300/1‑150 (MQ=255)
 gCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGt                                                                                                                                                      >  2:158389/1‑150 (MQ=255)
 gCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGt                                                                                                                                                      >  1:1007942/1‑150 (MQ=255)
   aaaCGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCa                                                                                                                                                    >  2:353403/1‑150 (MQ=255)
      cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAggg                                                                                                                                                 >  2:820512/1‑150 (MQ=255)
      cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAggg                                                                                                                                                 <  1:1036622/150‑1 (MQ=255)
       gCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGTGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGc                                                                                                                                                <  1:582221/150‑1 (MQ=255)
         tGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCgg                                                                                                                                              >  2:624519/1‑150 (MQ=255)
          ggtGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTAAACCGCCCGACAGGGCGGt                                                                                                                                             >  1:152086/3‑150 (MQ=255)
          gTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGt                                                                                                                                             >  1:1162652/1‑150 (MQ=255)
            tGGCGATCCCGATGTACTGTCTGTTTGAAATCGTTGTCTTCTTCTCTCGCTTTTAGGTTGGTAAAGGGCGAAATCGGGAAGTGGAAATCGACGCTGAAGCAGTAAGGGAAAAATCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTg                                                                                                                                           <  1:1227174/150‑1 (MQ=255)
            tGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTg                                                                                                                                           >  1:44022/1‑150 (MQ=255)
             ggCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTga                                                                                                                                          <  1:457477/150‑2 (MQ=255)
              gCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTc                                                                                                                                         <  2:1537062/150‑1 (MQ=255)
               cGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTgaat                                                                                                                                        <  1:156051/150‑4 (MQ=255)
               cGATCCCGATGTACTGTCTGTATGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTTGTAAAGGGCGAAATCGGTAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCa                                                                                                                                        <  2:1441225/150‑1 (MQ=255)
                gATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCat                                                                                                                                       <  2:1292037/150‑1 (MQ=255)
                 aTCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGTAAGAGGAAAACTACGCTTAAGCAGAAATCGAAAAAACTTAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCata                                                                                                                                      <  2:1371640/150‑1 (MQ=255)
                       aTGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGt                                                                                                                                <  2:1558434/150‑1 (MQ=255)
                       aTGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGt                                                                                                                                <  1:764566/150‑1 (MQ=255)
                           aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  1:1235436/150‑1 (MQ=255)
                           aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  1:1628300/150‑1 (MQ=255)
                           aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  2:262138/150‑1 (MQ=255)
                           aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  1:965043/150‑1 (MQ=255)
                             tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          <  2:1474523/150‑1 (MQ=255)
                             tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          <  1:1107843/150‑1 (MQ=255)
                             tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          >  2:869690/1‑150 (MQ=255)
                              gtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGAt                                                                                                                         >  2:1488674/1‑150 (MQ=255)
                                 tgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGtt                                                                                                                      >  1:1316111/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  2:558235/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  2:895480/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  1:1024739/1‑150 (MQ=255)
                                  gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGGCAGATGGAGTACAGGAGGttt                                                                                                                     >  2:94755/1‑150 (MQ=255)
                                   tttGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTg                                                                                                                    <  2:296120/150‑1 (MQ=255)
                                   tttGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTg                                                                                                                    <  2:1619674/150‑1 (MQ=255)
                                     tGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGat                                                                                                                  >  2:920500/1‑150 (MQ=255)
                                      gAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGata                                                                                                                 >  1:635473/1‑150 (MQ=255)
                                      gAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGata                                                                                                                 >  1:343151/1‑150 (MQ=255)
                                         aTCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCg                                                                                                              <  1:662678/150‑1 (MQ=255)
                                          tCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATTGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCgg                                                                                                             <  2:820379/150‑1 (MQ=255)
                                           cGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGc                                                                                                            <  1:716997/150‑1 (MQ=255)
                                            ggTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCg                                                                                                           <  1:91739/150‑1 (MQ=255)
                                              tgtCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGtt                                                                                                         <  1:675654/150‑1 (MQ=255)
                                                   tcttctCACGCTTTTACGTTGGAAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAAttt                                                                                                    <  1:1384300/150‑1 (MQ=255)
                                                    cttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTg                                                                                                   <  1:820512/150‑1 (MQ=255)
                                                    cttctCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTg                                                                                                   >  2:232099/1‑150 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                >  2:1442373/1‑150 (MQ=255)
                                                       ctcACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGAcc                                                                                                <  2:1291993/150‑1 (MQ=255)
                                                         cACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAg                                                                                              >  2:930159/1‑150 (MQ=255)
                                                         cACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAg                                                                                              >  2:1365096/1‑150 (MQ=255)
                                                           cGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGtt                                                                                            >  2:442644/1‑150 (MQ=255)
                                                              ttttACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGc                                                                                         <  1:717575/150‑1 (MQ=255)
                                                               tttACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCa                                                                                        >  1:1077148/1‑150 (MQ=255)
                                                               tttACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCa                                                                                        <  1:869690/150‑1 (MQ=255)
                                                                    gTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTg                                                                                   <  1:1442373/150‑1 (MQ=255)
                                                                      tGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCg                                                                                 <  2:178195/150‑1 (MQ=255)
                                                                      tGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCg                                                                                 <  2:1503969/150‑1 (MQ=255)
                                                                       ggTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATAAGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGa                                                                                >  2:685399/1‑150 (MQ=255)
                                                                         tAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGaaa                                                                              >  2:869590/1‑150 (MQ=255)
                                                                            aGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGAc                                                                           <  2:343151/150‑1 (MQ=255)
                                                                              ggCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCg                                                                         >  1:1455764/1‑150 (MQ=255)
                                                                               gCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGt                                                                        >  1:1038383/1‑150 (MQ=255)
                                                                               gCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGt                                                                        >  1:1325632/1‑150 (MQ=255)
                                                                                cGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtg                                                                       <  2:1024739/150‑1 (MQ=255)
                                                                                cGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtg                                                                       <  2:740710/150‑1 (MQ=255)
                                                                                 gAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtga                                                                      <  2:1007942/150‑1 (MQ=255)
                                                                                  aaaTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgat                                                                     <  1:1034004/150‑1 (MQ=255)
                                                                                   aaTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgatg                                                                    >  1:1531313/1‑150 (MQ=255)
                                                                                    aTCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGtgatga                                                                   >  1:847413/1‑150 (MQ=255)
                                                                                        ggAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGAtgtt                                                               >  2:355622/1‑150 (MQ=255)
                                                                                          aaGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGAtgttgt                                                             <  2:107507/150‑1 (MQ=255)
                                                                                              ggAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCg                                                         >  2:1370368/1‑150 (MQ=255)
                                                                                               gAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGt                                                        >  1:336002/1‑150 (MQ=255)
                                                                                                 aaaCGACGCTGAAGCAGCAAGAGAAAAAACTGAAGAATAAATTCAAACGCCCGTCAGGGCGGTTGTCATATCGAGTACAGGATGTTTGATATCGGCGTTAAGTTGACCAGATCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGAgc                                                      >  1:777889/1‑150 (MQ=255)
                                                                                                    cgacgCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCt                                                   <  2:1316111/150‑1 (MQ=255)
                                                                                                       cgCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCtttt                                                >  1:738022/1‑150 (MQ=255)
                                                                                                        gCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTg                                               <  2:1455764/150‑1 (MQ=255)
                                                                                                          tGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGAc                                             <  1:999147/150‑1 (MQ=255)
                                                                                                              gCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCgg                                         >  2:218832/1‑150 (MQ=255)
                                                                                                                   aaGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGtt                                    <  1:1354638/150‑1 (MQ=255)
                                                                                                                      cGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTTCGAAAGAACGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAAt                                 >  1:1271412/1‑150 (MQ=255)
                                                                                                                        aaaaaaCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGa                                                                                <  2:241273/101‑1 (MQ=255)
                                                                                                                        aaaaaaCTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGa                                                                                >  1:241273/1‑101 (MQ=255)
                                                                                                                              cTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACt                         >  2:1459349/1‑150 (MQ=255)
                                                                                                                                  agaaTAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACtcatc                     >  2:990178/1‑150 (MQ=255)
                                                                                                                                      tAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCg                 <  1:558235/150‑1 (MQ=255)
                                                                                                                                       aaaTTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCgg                >  1:65535/1‑150 (MQ=255)
                                                                                                                                        atgacAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGc               >  1:1320048/5‑150 (MQ=255)
                                                                                                                                        atgacAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGCATGTTTGATATCGGCGTTAATTTGACCACTTCGCAATTTGCGAAAGACCGTGATGATCTTGTCGCGAGCGCTTTTCACGCGGGAGTTCATGCGCTACTCATCACCGGc               >  2:1472576/5‑150 (MQ=255)
                                                                                                                                        aaTTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGc               <  1:1531710/150‑1 (MQ=255)
                                                                                                                                         tgacAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCa              >  1:1124617/4‑150 (MQ=255)
                                                                                                                                          ttCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCAc             <  2:1069180/150‑1 (MQ=255)
                                                                                                                                           acAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACt            >  1:288173/2‑150 (MQ=255)
                                                                                                                                              aCCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAAc         >  1:512796/1‑150 (MQ=255)
                                                                                                                                                cGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCt       >  2:537652/1‑150 (MQ=255)
                                                                                                                                                cGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCt       >  2:1369416/1‑150 (MQ=255)
                                                                                                                                                    cGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGt   >  1:419365/1‑150 (MQ=255)
                                                                                                                                                     gTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTg  <  1:537652/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 < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: