Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 777,624 T→C 5.2% intergenic (+17/‑116) tolA → / → tolB Tol‑Pal system protein TolA/Tol‑Pal system periplasmic protein TolB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913777,6240TC 5.2% 354.1 / 6.3 114intergenic (+17/‑116)tolA/tolBTol‑Pal system protein TolA/Tol‑Pal system periplasmic protein TolB
Reads supporting (aligned to +/- strand):  ref base T (50/58);  new base C (3/3);  total (53/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 117 total aligned reads displayed.

CAAACCTGAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGTTT  >  NC_000913/777475‑777771
                                                                                                                                                     |                                                                                                                                                   
cAAACCTGAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTTGACTTCAAACCGTAATCGCGATGTTGACTGtt                                                                                                                                                     >  2:547388/1‑150 (MQ=255)
   aCCTGAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCgg                                                                                                                                                  <  2:71516/150‑1 (MQ=255)
       gAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACgg                                                                                                                                              <  1:1152894/150‑1 (MQ=255)
       gAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCCCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACgg                                                                                                                                              <  1:77606/150‑1 (MQ=255)
       gAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGATCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACgg                                                                                                                                              <  1:1450740/150‑1 (MQ=255)
            tGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAAc                                                                                                                                         <  1:421528/150‑1 (MQ=255)
             ggCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACa                                                                                                                                        >  1:583937/1‑150 (MQ=255)
               cGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATc                                                                                                                                      <  1:1191024/150‑1 (MQ=255)
                gATCCTGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACCGTCCGAACAGTCAACATcg                                                                                                                                     <  1:790631/150‑2 (MQ=255)
                gATCCCGCACTTTGTGATGCTGAGTTGGCAGCAGCTAATCTTGCGAAGTTCGCGTAAGGACGATGCGAGGCAGTATATGAAGTGTTCAAAAACGCGCGATTGGATTTCAAAGCGTAATCGCGATGTTGACTGTTGGGACGGTCAACATCa                                                                                                                                     <  2:876198/150‑1 (MQ=255)
                gATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCa                                                                                                                                     <  1:1010972/150‑1 (MQ=255)
                gATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACCGTCCGAACAGTCAACATcg                                                                                                                                     <  1:1282415/150‑2 (MQ=255)
                 aTCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACCGTCCGAACAGTCAACATcgc                                                                                                                                    <  1:1638624/150‑3 (MQ=255)
                  tCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAgg                                                                                                                                   <  1:642100/150‑1 (MQ=255)
                   cccGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGc                                                                                                                                  >  1:1180750/1‑150 (MQ=255)
                   cccGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGc                                                                                                                                  >  2:1358534/1‑150 (MQ=255)
                      gCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGGAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCAcc                                                                                                                               <  2:1449683/150‑1 (MQ=255)
                             gTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGcc                                                                                                                        >  2:296857/1‑150 (MQ=255)
                               cAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCAc                                                                                                                      >  2:416468/1‑150 (MQ=255)
                                aGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACg                                                                                                                     >  1:807490/1‑150 (MQ=255)
                                 ggCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACgg                                                                                                                    >  1:972869/1‑150 (MQ=255)
                                 ggCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACgg                                                                                                                    >  1:384978/1‑150 (MQ=255)
                                   cTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACgggg                                                                                                                  <  1:352466/150‑1 (MQ=255)
                                        ttGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTg                                                                                                             <  2:1556085/150‑1 (MQ=255)
                                        ttGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTg                                                                                                             >  2:919613/1‑150 (MQ=255)
                                         tGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTgg                                                                                                            <  1:90650/150‑1 (MQ=255)
                                         tGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTgg                                                                                                            <  1:296857/150‑1 (MQ=255)
                                          ggCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGt                                                                                                           <  2:1514509/150‑1 (MQ=255)
                                           gcagcaGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTa                                                                                                          >  2:875890/1‑150 (MQ=255)
                                            cagcagCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAg                                                                                                         >  1:1635481/1‑150 (MQ=255)
                                            cagcagCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAg                                                                                                         >  2:544818/1‑150 (MQ=255)
                                             agcagcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAACCGCGCCATTGGACTTCAAACCGTAATCGCGATGTGGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGt                                                                                                        <  2:1263829/150‑1 (MQ=255)
                                             agcagcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGt                                                                                                        >  1:1076940/1‑150 (MQ=255)
                                                 gcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTtg                                                                                                    <  1:190293/150‑1 (MQ=255)
                                                 gcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTtg                                                                                                    >  1:1331899/1‑150 (MQ=255)
                                                  cTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTtgt                                                                                                   <  1:236251/150‑1 (MQ=255)
                                                   tAAACTTGCGATGATGCCTAAACCACGAAGCGAGGAATTATATGTAGTGTTCTAAATGGCGCGATTTGACTTCAATCCGTATTGGCTATTTTTACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACTGGGTTCTGGTAGTTTtgtg                                                                                                  <  2:450428/150‑1 (MQ=255)
                                                   tAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTtgtg                                                                                                  <  1:445400/150‑1 (MQ=255)
                                                   tAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTtgtg                                                                                                  <  2:1302654/150‑1 (MQ=255)
                                                   tAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAAAATCAGGCACCGGTTGCCACTGGGTTCTGGTTGTTTtgtg                                                                                                  <  2:852801/150‑1 (MQ=255)
                                                   tAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGAAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACAGGTTGCCACGGGGTTCTGGTAGTTTtgtt                                                                                                  >  1:773385/1‑149 (MQ=255)
                                                    aaaCTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTgtgt                                                                                                 >  2:416139/1‑150 (MQ=255)
                                                    aaaCTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTgtgt                                                                                                 >  1:1119252/1‑150 (MQ=255)
                                                    aaaCTTGCGAAGATCATTAAATCACTAATCCAGGCAGTATATTAATTTTTCAAAAACTCGCCATTTGACTTGAAACCGTTATCGAGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTgtgt                                                                                                 <  2:957589/150‑1 (MQ=255)
                                                       cTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTAtt                                                                                              >  2:107795/1‑150 (MQ=255)
                                                        ttGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTAttt                                                                                             >  2:1633989/1‑150 (MQ=255)
                                                           cGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAg                                                                                          <  1:687832/150‑1 (MQ=255)
                                                             aaGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGtt                                                                                        >  2:367927/1‑150 (MQ=255)
                                                               gATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTg                                                                                      >  1:1425861/1‑150 (MQ=255)
                                                                   ccGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTaa                                                                                  >  1:1586624/1‑150 (MQ=255)
                                                                      aaaCCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACAt                                                                               <  2:259553/150‑1 (MQ=255)
                                                                         ccacaaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGTCACCGGTTGCCACGGGGTTATGGTACTTTTGTGTATTTTAGTTTGTTAACATTCt                                                                            >  2:815879/1‑150 (MQ=255)
                                                                           accaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGc                                                                          <  2:1180750/150‑1 (MQ=255)
                                                                            ccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCt                                                                         >  2:69214/1‑150 (MQ=255)
                                                                             caAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTa                                                                        >  2:1541498/1‑150 (MQ=255)
                                                                                gCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAAt                                                                     <  1:286838/150‑1 (MQ=255)
                                                                                cccAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAAt                                                                     <  1:772318/149‑1 (MQ=255)
                                                                                    ggCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATc                                                                 <  1:614240/150‑1 (MQ=255)
                                                                                           tatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGtt                                                                                    >  1:1610567/1‑124 (MQ=255)
                                                                                           tatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGtt                                                                                    <  2:1610567/124‑1 (MQ=255)
                                                                                           tatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGc                                                                          <  1:1056269/134‑1 (MQ=255)
                                                                                           tatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGc                                                                          >  2:1056269/1‑134 (MQ=255)
                                                                                             tGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCAt                                                        <  1:1358534/150‑1 (MQ=255)
                                                                                             tGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCAt                                                        >  2:109835/1‑150 (MQ=255)
                                                                                                   gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTcc                                                  >  2:1082735/1‑150 (MQ=255)
                                                                                                    ttCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCa                                                 <  2:1308091/150‑1 (MQ=255)
                                                                                                     tCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAg                                                >  1:1094398/1‑150 (MQ=255)
                                                                                                      cAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAAt                                                                     <  1:595221/128‑1 (MQ=255)
                                                                                                      cAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGa                                               >  2:589501/1‑150 (MQ=255)
                                                                                                      cAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGa                                               >  1:853369/1‑150 (MQ=255)
                                                                                                      cAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGa                                               <  2:234395/150‑1 (MQ=255)
                                                                                                            cgcgCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAggg                                         <  1:367927/150‑1 (MQ=255)
                                                                                                              cgcCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGgag                                       <  2:1119252/150‑1 (MQ=255)
                                                                                                                 cATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGata                                    <  1:690264/150‑1 (MQ=255)
                                                                                                                    ttgACTTGAAACCGTAATCGATATGTTTAATGTTCGGACTGTCAACATCAGGGACCGTTTTCCACGGGGTTATGTTAGTTTTGTGTATTTTAGTTTGTTAACATTCTTATAAATTATCGTGGGCCATCTTTCCAGATAAGGGAGATatga                                 <  2:1635481/148‑1 (MQ=255)
                                                                                                                    tGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATatga                                 <  1:819359/150‑1 (MQ=255)
                                                                                                                    tGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATatga                                 <  1:919613/150‑1 (MQ=255)
                                                                                                                     ggACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATatgat                                >  1:520120/1‑150 (MQ=255)
                                                                                                                      gACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATatgatg                               <  1:1633989/150‑1 (MQ=255)
                                                                                                                          tCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGc                           >  2:76824/1‑150 (MQ=255)
                                                                                                                           cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGga                                        <  1:1580712/136‑1 (MQ=255)
                                                                                                                           cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCa                          <  2:18578/150‑1 (MQ=255)
                                                                                                                            aaaCCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAg                         <  2:1143627/150‑1 (MQ=255)
                                                                                                                            aaaCCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACAGTCTGCTTAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAg                         >  2:305158/1‑150 (MQ=255)
                                                                                                                              aCCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGc                       <  1:1082735/150‑1 (MQ=255)
                                                                                                                              aCCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGc                       >  1:1173401/1‑150 (MQ=255)
                                                                                                                                 gTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGTGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCAtt                    <  1:815879/150‑1 (MQ=255)
                                                                                                                                 gTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCAtt                    >  2:1492186/1‑150 (MQ=255)
                                                                                                                                    aTCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACg                 >  1:510349/1‑150 (MQ=255)
                                                                                                                                     tCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGa                <  2:583937/150‑1 (MQ=255)
                                                                                                                                      cctgATGTTGACCGTCCGAACAGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTACTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAg               >  2:757565/4‑150 (MQ=255)
                                                                                                                                       gcgATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGt              <  1:1077150/150‑1 (MQ=255)
                                                                                                                                        tgATGTTGACCGTCCGAACAGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTa             >  2:1590885/2‑150 (MQ=255)
                                                                                                                                        tgATGTTGACCGTCCGAACAGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTa             >  2:1191942/2‑150 (MQ=255)
                                                                                                                                              tGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCAttt       <  1:1541498/150‑1 (MQ=255)
                                                                                                                                                aCTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTgg     <  2:807490/150‑1 (MQ=255)
                                                                                                                                                 cTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGt    >  1:1354322/1‑150 (MQ=255)
                                                                                                                                                  tGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGtt   <  2:1425861/150‑1 (MQ=255)
                                                                                                                                                   gTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCa                      >  1:491037/1‑130 (MQ=255)
                                                                                                                                                   gTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCa                      <  2:491037/130‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAACCTGAAGGTGGCGATCCCGCACTTTGTCAGGCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGTTT  >  NC_000913/777475‑777771

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: