Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 632,127 A→C 9.1% intergenic (+128/‑55) cstA → / → ybdD carbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913632,1270AC 9.1% 302.4 / 16.5 116intergenic (+128/‑55)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base A (54/51);  new base C (6/5);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 118 total aligned reads displayed.

CAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATAT  >  NC_000913/631979‑632276
                                                                                                                                                    |                                                                                                                                                     
ccggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                                                                                      <  1:499786/148‑1 (MQ=255)
 aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAgg                                                                                                                                                     <  1:721329/150‑1 (MQ=255)
   ccAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGAAGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCg                                                                                                                                                   <  2:275709/149‑1 (MQ=255)
    cAAAAGGCGCACACTAAAGGCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGt                                                                                                                                                  <  1:1892988/150‑1 (MQ=255)
      aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                <  1:1782253/150‑1 (MQ=255)
      aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                <  1:1571642/150‑1 (MQ=255)
         ggCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                             >  1:754976/1‑150 (MQ=255)
         ggCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                             >  1:1461315/1‑150 (MQ=255)
           cgcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                           <  2:1684107/150‑1 (MQ=255)
           cgcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                           >  1:1272802/1‑150 (MQ=255)
            gcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                          >  2:1062220/1‑150 (MQ=255)
                 cTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                     >  2:690750/1‑150 (MQ=255)
                 cTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGATCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  1:1380614/150‑1 (MQ=38)
                  aaaaGGCAGCGGGCGGGGGGGATGTTGGCGAATGCTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                    <  1:1994376/149‑1 (MQ=39)
                    acgTCAGAGGGAGGGTGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTTATAGCATTTGTAGGCGGGATAAGTCGTTCACGCCGCATCCGGc                                                                                                                                  <  2:1246037/148‑1 (MQ=39)
                    aaGTCAGCGGGCGGGGGGGATGTTGGCGACTGTTGGCGTAGGGCCCAGGGGTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                  <  1:369583/150‑1 (MQ=37)
                             ggAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCa                                                                                                                         <  2:75109/149‑1 (MQ=255)
                               aGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGc                                                                                                                       <  2:1127846/150‑1 (MQ=255)
                                  gggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAAc                                                                                                                    <  1:1781316/150‑1 (MQ=255)
                                        tGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGtt                                                                                                              >  1:1237475/1‑150 (MQ=255)
                                           tGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGc                                                                                                           >  2:1304997/1‑150 (MQ=255)
                                             gCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCgg                                                                                                         <  1:1109919/150‑1 (MQ=255)
                                                aaTGTTGGCTTAGTGCCGAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGGTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAgg                                                                                                      <  2:1807325/150‑1 (MQ=255)
                                                  tGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGcc                                                                                                    <  2:1212643/150‑1 (MQ=255)
                                                         ttAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGt                                                                                             <  1:1650464/150‑1 (MQ=255)
                                                          tAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGtt                                                                                            <  2:1461315/150‑1 (MQ=255)
                                                          tAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGtt                                                                                            <  2:754976/150‑1 (MQ=255)
                                                          tAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACGCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGtt                                                                                            >  1:1681726/1‑150 (MQ=255)
                                                              gCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                                                                        <  1:1543842/150‑1 (MQ=255)
                                                               cccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAtt                                                                                       <  1:1462999/150‑1 (MQ=255)
                                                               cccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAtt                                                                                       >  2:1611930/1‑150 (MQ=255)
                                                                ccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTc                                                                                      <  2:1390906/150‑1 (MQ=255)
                                                                ccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTc                                                                                      <  2:1950398/150‑1 (MQ=255)
                                                                  aGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCAc                                                                                    >  1:1607968/1‑150 (MQ=255)
                                                                   gggTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATATCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                                                                                   <  1:1025594/150‑1 (MQ=255)
                                                                   gggTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                                                                                   >  1:1642187/1‑150 (MQ=255)
                                                                      ttCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTTACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGc                                                                                <  1:689665/150‑1 (MQ=255)
                                                                        cccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaa                                                                              >  2:921490/1‑150 (MQ=255)
                                                                         ccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaaa                                                                             >  1:495162/1‑150 (MQ=255)
                                                                          ctctCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaaaa                                                                            >  2:350881/1‑150 (MQ=255)
                                                                              cACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCATCATGGTTAGCGGAGGCCAAGATGTTTGATTCTCTGGCAAAAGCCg                                                                        >  1:1576046/1‑150 (MQ=255)
                                                                              cACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCg                                                                        >  2:863155/1‑150 (MQ=255)
                                                                               aCCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGCCCGCGCGCGTTGATAGCATGTGTAGGCCGGATAAGGAGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                                                                       <  1:299333/150‑1 (MQ=255)
                                                                                cccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGa                                                                      >  1:560834/1‑150 (MQ=255)
                                                                                cccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGa                                                                      >  2:173564/1‑150 (MQ=255)
                                                                                cccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGa                                                                      >  2:872075/1‑150 (MQ=255)
                                                                                 ccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                                                                     <  1:1496087/150‑1 (MQ=255)
                                                                                   tAACCCTCTCCCCGGTGGGGCGAGGGGACTGCCCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCATGATGTTTGATTCACTGGCAAAAGCCGGaaaa                                                                   >  1:1877996/1‑150 (MQ=255)
                                                                                   tAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaaa                                                                   >  1:1109175/1‑150 (MQ=255)
                                                                                     ccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGCCAAAAGCCCGAAAata                                                                 >  2:1992027/2‑150 (MQ=38)
                                                                                      cccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGCTAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCCTCTTTCAGCAACATGGGTAGCGGAGGCCAAGATGTTTGAGTCACTGGCAAACGCCGGAAAattt                                                                >  1:1224136/1‑148 (MQ=38)
                                                                                       ccTCTCCCCGGTGGGGCGAAGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACAGGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt                                                               >  2:919712/1‑150 (MQ=255)
                                                                                         tctcCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTa                                                             <  2:473716/150‑1 (MQ=255)
                                                                                            ccccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACCGGCAAAAGCCGGAAAATATTTAGGt                                                          >  2:934375/1‑150 (MQ=255)
                                                                                             cccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTc                                                         >  1:1966186/1‑150 (MQ=255)
                                                                                                    gggCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCTCTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcggc                                                  >  1:142458/1‑150 (MQ=255)
                                                                                                    gggCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcggc                                                  <  2:1116763/150‑1 (MQ=255)
                                                                                                     ggCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATACGGCGGTCACGCCGCATCCGGCACTCTTTCAGCAACATGGGTAGCGGAGGCCAAGATGTTTGATTCTCTGGCAAAAGCCGGAAAATAGTTAGGTCAGgcggcg                                                 >  1:307929/1‑150 (MQ=255)
                                                                                                     ggCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGgcggcg                                                 >  1:220890/1‑150 (MQ=255)
                                                                                                             gACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgat                                         >  2:1405925/1‑150 (MQ=255)
                                                                                                               cTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                       >  1:611960/1‑150 (MQ=255)
                                                                                                               cTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                       >  1:1489958/1‑150 (MQ=255)
                                                                                                               cTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGAATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                       <  1:266442/150‑1 (MQ=255)
                                                                                                                 ttccGAGCGCGTTGATAGCATTTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGAtt                                     <  1:309344/148‑1 (MQ=255)
                                                                                                                 gACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGAtt                                     >  2:578698/1‑150 (MQ=255)
                                                                                                                  cccGAGCTCGTTGATATCATTTGTAGGCCGGATAAGGAGTTCACGCCTCATCCTGCACTCTTTCAGCAACATGGTTAGCGGAGGCCATGATGTTTTATTCACTGGCAAAAGCCGGAAAATATATAGGTCAGGCGGCGAAGCTGATGATTg                                    <  2:448587/149‑1 (MQ=255)
                                                                                                                  aCCGTGCGCGTTGATTGCATTTGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGAACTCTTTCAGCAACATGCTTTGCGGAGTCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAATTATTTAGGTCAGGCGGCGAAGCTGATGATTg                                    <  1:114066/150‑1 (MQ=255)
                                                                                                                   ccGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                   <  2:2003762/150‑1 (MQ=255)
                                                                                                                    cGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt                                  >  1:269095/1‑150 (MQ=255)
                                                                                                                    cGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt                                  >  1:1445754/1‑150 (MQ=255)
                                                                                                                        cgcgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGc                              >  2:1327976/1‑150 (MQ=255)
                                                                                                                        cgcgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGc                              >  1:224624/1‑150 (MQ=255)
                                                                                                                          cgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCt                            >  1:191686/1‑150 (MQ=255)
                                                                                                                          cgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCt                            >  1:826614/1‑150 (MQ=255)
                                                                                                                           gTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTg                           >  1:185041/1‑150 (MQ=255)
                                                                                                                           gTTGATAGCATTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCACTCTTTCAGCCACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTg                           <  1:1531288/150‑1 (MQ=255)
                                                                                                                            ttGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                          >  1:1243989/1‑150 (MQ=255)
                                                                                                                              gATAGCATTTGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAtt                        <  2:1208302/150‑1 (MQ=255)
                                                                                                                                tAGCATTTGTTGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGTCCTAGATGTTTGATTCACTGTCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                      <  1:1426290/150‑1 (MQ=255)
                                                                                                                                tAGCAGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTTGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                      <  2:1821405/150‑1 (MQ=255)
                                                                                                                                   cATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATGCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGAc                   >  1:833382/1‑150 (MQ=255)
                                                                                                                                   cATTTGTAGGCAGGATAAGGAGTTCACGCCGCATCCGGCACTGTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAACCAGGAAAATTTTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTTATTACGAc                   <  2:565532/150‑1 (MQ=255)
                                                                                                                                        gagtgCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTa              >  2:717018/5‑150 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTa              <  2:1839461/150‑1 (MQ=255)
                                                                                                                                          ttgCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGAAAACTATg            >  1:189419/3‑150 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGTGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATg            >  2:1345851/1‑150 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATg            >  2:1994975/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGt           >  2:1558443/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTc          >  2:1971551/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTc          >  1:1272905/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTc          <  1:1062220/150‑1 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTc          >  1:1755381/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTGCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCg         <  2:1216615/150‑1 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCg         <  2:1154324/150‑1 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCg         >  2:272605/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAATCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACAATGTCg         >  2:1027577/1‑150 (MQ=255)
                                                                                                                                                gATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTTGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTAGTATGCCTGATTACGACAACTATGTCGAAc      <  1:134534/150‑1 (MQ=255)
                                                                                                                                                 aTACGGCGTTCACGCCGCTTCCGGCACTCTTTCAGCAACATGGATAGCGGAGGCCAAGATGTTTGATTCACTGTCAAAAGGCGGAAAATATTTAGGTCAGGCGTCGAAGCTGATGATTTGTATGCCTGATTACGACAACTATGTCGAACa     <  1:984396/150‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACata   >  1:73816/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACatat  >  1:604694/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATAT  >  NC_000913/631979‑632276

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