Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,893 G→A 5.1% intergenic (+28/+12) insQ → / ← ydcO IS609 transposase B/BenE family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,8930GA 5.1% 376.6 / 5.5 117intergenic (+28/+12)insQ/ydcOIS609 transposase B/BenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (57/54);  new base A (2/4);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 118 total aligned reads displayed.

AATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTC  >  NC_000913/1504749‑1505042
                                                                                                                                                |                                                                                                                                                     
aaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGTCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgg                                                                                                                                                  >  1:81642/1‑150 (MQ=255)
aaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGAAGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgg                                                                                                                                                  >  2:614725/1‑150 (MQ=255)
  tGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggag                                                                                                                                                <  2:1036307/150‑1 (MQ=255)
  tGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggag                                                                                                                                                <  2:1707622/150‑1 (MQ=255)
  tGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTCTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggat                                                                                                                                                <  1:140036/150‑2 (MQ=255)
   ggCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCCGGGCGTggagg                                                                                                                                               >  1:1279414/1‑150 (MQ=255)
   ggCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggatc                                                                                                                                               <  2:711077/150‑3 (MQ=255)
   ggCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGGGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggatc                                                                                                                                               <  2:1123968/150‑3 (MQ=255)
      gcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATg                                                                                                                                            <  2:1042052/150‑1 (MQ=255)
      gcTCGCAAAATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGCTGCAGTCAGGCCGCCCGTTGCAGCAGGCACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGACCTTCAGGGCGTGGAGGATg                                                                                                                                            >  1:1700646/1‑150 (MQ=255)
       cTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGt                                                                                                                                           >  2:1667874/1‑150 (MQ=255)
       cTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggatccctg                                                                                                                                           <  2:1246375/150‑7 (MQ=255)
         cGCAACATTTTAGCGGCGGGGCCCGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCAGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCa                                                                                                                                         <  2:1729518/150‑1 (MQ=255)
          gCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCaa                                                                                                                                        >  1:1523398/1‑150 (MQ=255)
           cAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                                                                                                                       >  1:375645/1‑150 (MQ=255)
           cAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                                                                                                                       >  2:206332/1‑150 (MQ=255)
           cAACATTTTAGCGGCGGGGCACGCCGTTATTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                                                                                                                       >  1:1578036/1‑150 (MQ=255)
                ttttAGCGGCGGGGCAGGCCGTTCTTGCCTGTGGAGCGATGGGGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGTTCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCgg                                                                                                                                  <  2:135189/150‑1 (MQ=255)
                ttttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCgg                                                                                                                                  <  1:804129/150‑1 (MQ=255)
                ttttAGCGGCGGGGCACGCCGTTCTTGCCGTTGGAGAGCTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGACATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGGGGAGGATGTCAATATCgg                                                                                                                                  <  2:148302/150‑1 (MQ=255)
                tattAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGTAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATggg                                                                                                                                  <  2:800433/148‑1 (MQ=255)
                  ttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                                                                                                                <  1:390522/150‑1 (MQ=255)
                   tAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGttt                                                                                                                               >  1:163474/1‑150 (MQ=255)
                    aGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTc                                                                                                                              <  2:1194234/150‑1 (MQ=255)
                      cggcggGGCACGCCGTTCTTGCCTGGGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCTCCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTg                                                                                                                            <  2:348072/150‑1 (MQ=255)
                       ggcggGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGt                                                                                                                           <  1:1594410/150‑1 (MQ=255)
                        gcggGGCACGCCGTGCTTGCCTCTGGCGTGATGGGGCAGTCCGGCCGCCCGGTGCAGCAGGAACCCACCGAAATGAATCAGGCGACAGCATGAACGTAGCAGGGATCCCCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTc                                                                                                                          <  2:340521/150‑1 (MQ=255)
                          ggggCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAg                                                                                                                        <  2:1521267/150‑1 (MQ=255)
                           gggCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGc                                                                                                                       >  1:971910/1‑150 (MQ=255)
                            ggCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCg                                                                                                                      >  1:504916/1‑150 (MQ=255)
                            ggCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCg                                                                                                                      >  2:621508/1‑150 (MQ=255)
                              cACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAt                                                                                                                    <  1:1387171/150‑1 (MQ=255)
                                cGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATTGGTTTCTGTCAGCGATTa                                                                                                                  <  2:160733/150‑1 (MQ=255)
                                cGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTa                                                                                                                  <  2:791491/150‑1 (MQ=255)
                                   cGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAAt                                                                                                               <  2:1099808/150‑1 (MQ=255)
                                   cGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAAt                                                                                                               <  2:1206164/150‑1 (MQ=255)
                                     ttCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                             <  2:713646/150‑1 (MQ=255)
                                        ttGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGCACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac                                                                                                          <  1:1111756/150‑1 (MQ=255)
                                             tgtgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcaccac                                                                                                     <  1:176365/150‑1 (MQ=255)
                                                  agagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCATCGATTAAATTCAACACCACTTAAc                                                                                                >  2:393581/1‑150 (MQ=255)
                                                  agagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATACAACACCACGTAAc                                                                                                >  2:386895/1‑150 (MQ=255)
                                                   gagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAGTAAATTCAACACCACGTAACa                                                                                               >  2:1293950/1‑150 (MQ=255)
                                                    agaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACaa                                                                                              <  2:1564888/150‑1 (MQ=255)
                                                       tGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACg                                                                                           >  1:1433785/1‑150 (MQ=255)
                                                        ggTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                          >  2:1131555/1‑150 (MQ=255)
                                                        ggTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                          >  2:194679/1‑150 (MQ=255)
                                                          tGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCt                                                                                        >  2:1718938/1‑150 (MQ=255)
                                                           gCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGCCAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                                                                                       >  2:1535530/1‑150 (MQ=255)
                                                               tCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGc                                                                                   <  2:1467014/150‑1 (MQ=255)
                                                                 aGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCaa                                                                                 <  1:399987/150‑1 (MQ=255)
                                                                 aGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAAAATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCaa                                                                                 >  1:1675414/1‑150 (MQ=255)
                                                                  ggCCGCCCGTTGAAGCAGGACCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACGCCACGTAACAAACGCCTCCGGCAAt                                                                                >  2:531880/1‑150 (MQ=255)
                                                                  ggCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAGCGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAAt                                                                                >  2:1150263/1‑150 (MQ=255)
                                                                  ggCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAAt                                                                                <  2:1518897/150‑1 (MQ=255)
                                                                  ggCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAAt                                                                                <  1:1710800/150‑1 (MQ=255)
                                                                   gccgccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACATCCTGAACGTACCAGGGATCCACGTCCATCAGAGAGTGGAGGAGGTCAATATCGTATTCTTAAAGAGAATAAATTAAAAACCACGTAACAAACGCCTCCCGCAAtt                                                                               >  2:1211609/1‑150 (MQ=255)
                                                                     cgccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTaa                                                                             >  2:1226297/1‑150 (MQ=255)
                                                                     cgccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTaa                                                                             >  2:1421974/1‑150 (MQ=255)
                                                                        ccGTTGAAGCATGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAATAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc                                                                          <  2:881624/150‑1 (MQ=255)
                                                                        ccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc                                                                          <  2:81642/150‑1 (MQ=255)
                                                                         cGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATccc                                                                         <  2:185255/150‑1 (MQ=255)
                                                                         cGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATccc                                                                         <  1:1588562/150‑1 (MQ=255)
                                                                           ttGAAGCAGGAACCCACCGATATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCa                                                                       <  1:281432/150‑1 (MQ=255)
                                                                             gAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCaaa                                                                     >  1:203952/1‑150 (MQ=255)
                                                                               aGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCTATTAATCCCCAAAAc                                                                   >  1:162719/1‑150 (MQ=255)
                                                                                   ggAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAg                                                               <  2:1543124/150‑1 (MQ=255)
                                                                                      aCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAAc                                                            <  2:293614/150‑1 (MQ=255)
                                                                                       cccACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACAACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAAcc                                                           >  2:1243687/1‑150 (MQ=255)
                                                                                               aaTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCgc                                                   >  1:1222809/1‑149 (MQ=255)
                                                                                                aTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACGCCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAc                                                  >  2:1004171/1‑150 (MQ=255)
                                                                                                 tGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCACCACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAcc                                                 >  1:888217/1‑150 (MQ=255)
                                                                                                 tGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAcc                                                 >  1:1178608/1‑150 (MQ=255)
                                                                                                   aTTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAg                                               >  1:1583170/1‑150 (MQ=255)
                                                                                                   aTTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAg                                               >  2:652040/1‑150 (MQ=255)
                                                                                                     tCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCg                                             >  2:1531130/1‑150 (MQ=255)
                                                                                                      cAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGATGTTGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGt                                            <  1:1337311/150‑1 (MQ=255)
                                                                                                            aCAGCCTGAACGAATCAGGCATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGAAACTCTCAGCGATAAAATTCAACACCACCTAACAAACGCCTCCTGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAAAcc                                      <  1:1091020/150‑1 (MQ=255)
                                                                                                               gCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACATACGCCTCCGGCAATTAATCCCCAAACCGCAGAACCAATCCCGCCCAGCGTCAATCCACt                                   >  2:1821289/1‑150 (MQ=255)
                                                                                                               gCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACt                                   <  1:1473144/150‑1 (MQ=255)
                                                                                                                       gTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTa                           <  1:206332/150‑1 (MQ=255)
                                                                                                                       gTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTACTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGCCCAGCGTCAATCCACTTGCCGTTa                           >  1:812924/1‑150 (MQ=255)
                                                                                                                          gCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCa                        <  1:1718938/150‑1 (MQ=255)
                                                                                                                           cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCACCACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGCCCAGCGTCAATCCACTTGCCGTTACCAg                       >  1:1609367/1‑150 (MQ=255)
                                                                                                                           cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAg                       >  1:1602353/1‑150 (MQ=255)
                                                                                                                           cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAg                       <  2:1725793/150‑1 (MQ=255)
                                                                                                                            aGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGa                      <  2:341702/150‑1 (MQ=255)
                                                                                                                               gATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGACCCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAAt                   >  1:992848/1‑150 (MQ=255)
                                                                                                                               gATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCAAGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAAt                   >  1:1528969/1‑150 (MQ=255)
                                                                                                                                 tCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAAACAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGc                 >  1:1423792/1‑150 (MQ=255)
                                                                                                                                 tCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTATGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGc                 <  2:1461274/150‑1 (MQ=255)
                                                                                                                                 tCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGc                 >  1:1515351/1‑150 (MQ=255)
                                                                                                                                  ccACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGCACCAATCCCGACCAGCGTCAATCCACGTGCCGTTCCCAGAAATGcc                >  2:1611861/1‑150 (MQ=255)
                                                                                                                                   aaCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGTGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca               <  1:54523/149‑1 (MQ=255)
                                                                                                                                    aCGTCCTTGAGGGCGTGGAGGATGTCAATATCGGTTTGTGTCAGCTATTAAATTCAACAACAAGTAACATACGCCTCCGGCAATTTATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccac              <  2:852901/150‑1 (MQ=255)
                                                                                                                                    aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccac              <  1:181353/150‑1 (MQ=255)
                                                                                                                                     cGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCATAACCAATCCCTACCAGAGTCAATCCACTTGCCGTTACCAGAAAAGccacc             >  1:1097012/1‑150 (MQ=255)
                                                                                                                                        ccTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCaccacc          >  1:346458/1‑150 (MQ=255)
                                                                                                                                          ttCAGGGCGTGCAGGATGTCAATATCGGTTTCTGTCATCGATTAAATTCAAACCCACGTACCAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAAGACCGACCCCCGACAATCCACTTGCCGTTCCCAGAAATGCCACCAccgc        >  1:863134/1‑150 (MQ=255)
                                                                                                                                             aGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACcgccgc     >  2:1769466/1‑150 (MQ=255)
                                                                                                                                                gCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCATCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGtc  <  1:1221752/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
AATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTC  >  NC_000913/1504749‑1505042

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