Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,023,526 T→G 7.3% intergenic (+14/‑28) 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,5260TG 7.3% 310.1 / 13.5 108intergenic (+14/‑28)tatC/tatDtwin arginine protein translocation system ‑ TatC protein/3' ‑> 5' ssDNA/RNA exonuclease TatD
Reads supporting (aligned to +/- strand):  ref base T (51/48);  new base G (5/3);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 110 total aligned reads displayed.

GCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTG  >  NC_000913/4023377‑4023674
                                                                                                                                                     |                                                                                                                                                    
gCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGt                                                                                                                                                      >  2:158389/1‑150 (MQ=255)
gCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGt                                                                                                                                                      >  1:1007942/1‑150 (MQ=255)
  aaaCGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCa                                                                                                                                                    >  2:353403/1‑150 (MQ=255)
     cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAggg                                                                                                                                                 <  1:1036622/150‑1 (MQ=255)
     cGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAggg                                                                                                                                                 >  2:820512/1‑150 (MQ=255)
      gCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGTGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGc                                                                                                                                                <  1:582221/150‑1 (MQ=255)
        tGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCgg                                                                                                                                              >  2:624519/1‑150 (MQ=255)
         gTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGt                                                                                                                                             >  1:1162652/1‑150 (MQ=255)
           tGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCTGACGGGCGGTTg                                                                                                                                           <  1:470528/150‑1 (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                                                                                                                                <  1:764566/150‑1 (MQ=255)
                      aTGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGt                                                                                                                                <  2:1558434/150‑1 (MQ=255)
                          aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  1:1235436/150‑1 (MQ=255)
                          aCTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAg                                                                                                                            <  1:1628300/150‑1 (MQ=255)
                           ctgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAgg                                                                                                                           >  2:1531710/1‑150 (MQ=255)
                            tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          >  2:869690/1‑150 (MQ=255)
                            tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          <  2:1474523/150‑1 (MQ=255)
                            tgtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGa                                                                                                                          <  1:1107843/150‑1 (MQ=255)
                             gtctgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGAt                                                                                                                         >  2:1488674/1‑150 (MQ=255)
                                tgtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGtt                                                                                                                      >  1:1316111/1‑150 (MQ=255)
                                 gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  2:895480/1‑150 (MQ=255)
                                 gtTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGttt                                                                                                                     >  2:558235/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:343151/1‑150 (MQ=255)
                                     gAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGata                                                                                                                 >  1:635473/1‑150 (MQ=255)
                                      aaaTCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGatat                                                                                                                >  1:885004/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:1365096/1‑150 (MQ=255)
                                                        cACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAg                                                                                              >  2:930159/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)
                                                                      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)
                                                                                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)
                                                                                             ggAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCg                                                         >  2:451972/1‑150 (MQ=255)
                                                                                              gAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGt                                                        >  1:336002/1‑150 (MQ=255)
                                                                                                aaaCGACGCTGAAGCAGCAAGAGAAAAAACTGAAGAATAAATTCAAACGCCCGTCAGGGCGGTTGTCATATCGAGTACAGGATGTTTGATATCGGCGTTAAGTTGACCAGATCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGAgc                                                      >  1:777889/1‑150 (MQ=255)
                                                                                                aaaCGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTgc                                                      >  2:849543/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)
                                                                                                          gAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGAcg                                            >  1:1152357/1‑150 (MQ=255)
                                                                                                             gCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCgg                                         >  2:218832/1‑150 (MQ=255)
                                                                                                              cAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCggg                                        <  2:1182174/150‑1 (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:200813/1‑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)
                                                                                                                                           cAACCGCCCTGACGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTAGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTa           >  2:142722/1‑150 (MQ=255)
                                                                                                                                             aCCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAAc         >  1:512796/1‑150 (MQ=255)
                                                                                                                                               cGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCt       >  2:537652/1‑150 (MQ=255)
                                                                                                                                                 cccGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGc     >  2:226539/1‑150 (MQ=255)
                                                                                                                                                   cGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGt   >  1:419365/1‑150 (MQ=255)
                                                                                                                                                    gTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTg  <  1:537652/150‑1 (MQ=255)
                                                                                                                                                    gTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTg  <  1:1369416/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAACGCTGTTGGCGATCCCGATGTACTGTCTGTTTGAAATCGGTGTCTTCTTCTCACGCTTTTACGTTGGTAAAGGGCGAAATCGGGAAGAGGAAAACGACGCTGAAGCAGAAAGCGAAAAAACTGAAGAATAAATTCAACCGCCCGTCAGGGCGGTTGTCATATGGAGTACAGGATGTTTGATATCGGCGTTAATTTGACCAGTTCGCAATTTGCGAAAGACCGTGATGATGTTGTAGCGTGCGCTTTTGACGCGGGAGTTAATGGGCTACTCATCACCGGCACTAACCTGCGTG  >  NC_000913/4023377‑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: