Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,893 G→A 5.3% 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.3% 238.2 / 2.5 76intergenic (+28/+12)insQ/ydcOIS609 transposase B/BenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base A (2/2);  total (35/41)
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.06e-01

AATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGC  >  NC_000913/1504749‑1505039
                                                                                                                                                |                                                                                                                                                  
aaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGt                                                                                                                                                 <  1:123549/148‑1 (MQ=255)
aaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                                                                                                                >  6:469365/1‑149 (MQ=255)
aaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                                                                                                                >  3:369410/1‑149 (MQ=255)
  tGGCGCTCGCAACATTTTAGCGGTGGGGGACGCCTTTATTGCATGTGGAGAGATGGTGCAGTCTGGCCGCCCGTTGAAGCAGGAACCCTCCGAAATGTTTCAGGCGACATCCTGTACGTAGCAGGGATCCACGTCCTTCAGGTCGTgga                                                                                                                                              <  6:30537/149‑1 (MQ=255)
    gcgcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggat                                                                                                                                             <  6:315534/148‑2 (MQ=255)
     cgcCCGCAACATTTTCGCGGCGGGGCACGCCGGTCTTGCCTGTGGAGCGCTGGCGCAGTCAGGCCGCCCGTTGACGCAGGACCCCCCCGAACTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggagga                                                                                                                                           <  4:271955/149‑1 (MQ=255)
       cTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATg                                                                                                                                         >  1:12709/1‑149 (MQ=255)
       cTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATg                                                                                                                                         >  5:343686/1‑149 (MQ=255)
       cTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATg                                                                                                                                         >  2:175979/1‑149 (MQ=255)
          gCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTc                                                                                                                                       >  7:320397/1‑148 (MQ=255)
          gCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTc                                                                                                                                       <  1:396820/148‑1 (MQ=255)
            aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGt                                                                                                                                                 >  8:242627/1‑136 (MQ=255)
            aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGt                                                                                                                                                 <  7:242627/136‑1 (MQ=255)
               aTTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCACCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATc                                                                                                                                 >  1:251723/1‑149 (MQ=255)
                ttttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCg                                                                                                                                <  8:68018/149‑1 (MQ=255)
                      cggcggGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCt                                                                                                                          <  6:21592/149‑1 (MQ=255)
                          ggggCACGCCGTTCTTGCCTGTGGAGAGAGGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                                                                                                                      <  7:454528/149‑1 (MQ=255)
                               aCGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGAt                                                                                                                                          <  5:25292/124‑1 (MQ=255)
                               aCGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGAt                                                                                                                                          >  6:25292/1‑124 (MQ=255)
                                cGCCGTTCTTGCCTGTGGAGAGATGGTGCTGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGGCTGAACGTAGCAGGGATCCTCGTCCTTCAGGTCGTGGTGGATGGCAATATCGGTTTCTGTCAGCGAtt                                                                                                                <  3:400663/149‑1 (MQ=255)
                                  ccGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTaa                                                                                                              <  8:320397/149‑1 (MQ=255)
                                      tCTTGCCTGTGGAGAGATTGTGCAGTCAGGCCGCCCGTTGAAGCAGGCACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                          <  2:223847/149‑1 (MQ=255)
                                        ttGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCACGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa                                                                                                        <  4:369410/149‑1 (MQ=255)
                                          gCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca                                                                                                      >  2:26190/1‑149 (MQ=255)
                                               tgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcaccac                                                                                                  <  7:326861/148‑1 (MQ=255)
                                               tgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAGTAACTTCAACACCACg                                                                                                 >  6:86681/1‑149 (MQ=255)
                                                    agaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAAACCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAc                                                                                             >  6:224639/1‑148 (MQ=255)
                                                       tGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAAc                                                                                         <  7:368838/149‑1 (MQ=255)
                                                       tGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAAc                                                                                         <  5:469365/149‑1 (MQ=255)
                                                          tGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                       >  8:223570/1‑148 (MQ=255)
                                                            cAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                                                                                    <  6:343686/149‑1 (MQ=255)
                                                            cAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                                                                                    <  1:175979/149‑1 (MQ=255)
                                                              gTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCg                                                                                  >  4:151197/1‑149 (MQ=255)
                                                               tCAGGCCGCCCGTCGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCgg                                                                                 >  1:498512/1‑149 (MQ=255)
                                                                 aGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa                                                                                                        <  7:309991/124‑1 (MQ=255)
                                                                 aGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa                                                                                                        >  8:309991/1‑124 (MQ=255)
                                                                    ccgccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAAtt                                                                            <  2:12709/149‑1 (MQ=255)
                                                                        ccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATc                                                                        <  7:477027/149‑1 (MQ=255)
                                                                         cGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc                                                                       >  3:146159/1‑149 (MQ=255)
                                                                           ttGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                       <  2:353846/131‑1 (MQ=255)
                                                                           ttGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                       >  1:353846/1‑131 (MQ=255)
                                                                                           ccGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGTGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAAt                                                     <  1:26133/149‑1 (MQ=255)
                                                                                            cGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATc                                                    <  2:498512/149‑1 (MQ=255)
                                                                                             gAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCACCACCACGTAACAAACGCCTCCGGCAAGTAATCCCCAAAACGCAGAACCAATcc                                                   >  6:417289/1‑149 (MQ=255)
                                                                                              aaaTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATccc                                                  >  1:107095/1‑149 (MQ=255)
                                                                                              aaaTGATTCAGGCGACAGCCTGAAAGTAGCAGGGATACCCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTTAATTCAACACCACGTAACAAACGCCTCAGGCAATTAATCCACAAAACGCAGAACCAATccc                                                  >  7:385431/1‑149 (MQ=255)
                                                                                                  gATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCgcc                                               >  2:486979/1‑146 (MQ=255)
                                                                                                    ttCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAACCGCAGACCCAATCCCGCCCAg                                            >  6:71640/1‑149 (MQ=255)
                                                                                                       aGGCGACAGCCTGAACGTAGGAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAGACCACGTAACTAACGCCTCCGGCAATTAATCCCCATAACGGAGAACCTATCCCGACCAGCGt                                         <  6:307658/149‑1 (MQ=255)
                                                                                                       aGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGa                                                           >  2:479517/1‑131 (MQ=255)
                                                                                                       aGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGa                                                           <  1:479517/131‑1 (MQ=255)
                                                                                                           gACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTTACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAAt                                     <  8:6064/149‑1 (MQ=255)
                                                                                                           gACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAAt                                     <  7:223570/149‑1 (MQ=255)
                                                                                                           gACAGCCTGAACGTAGCAGGGATCCACGTCATTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGACTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAAt                                     <  4:13361/149‑1 (MQ=255)
                                                                                                                     aCGTAGCAGGGATCCACGTCATTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCg                           <  2:123408/149‑1 (MQ=255)
                                                                                                                      cGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGTTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGt                          <  1:190329/149‑1 (MQ=255)
                                                                                                                        tAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTa                        >  4:454974/1‑149 (MQ=255)
                                                                                                                          gacatcctCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAc                       >  6:302986/8‑148 (MQ=255)
                                                                                                                          gCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAc                       <  2:107095/148‑1 (MQ=255)
                                                                                                                           cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCa                     <  3:151197/149‑1 (MQ=255)
                                                                                                                            catcctCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAg                    <  5:302986/144‑1 (MQ=255)
                                                                                                                            aGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAc                       >  6:415831/1‑146 (MQ=255)
                                                                                                                               gATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGaaa                 >  1:58805/1‑149 (MQ=255)
                                                                                                                                aTCCACGTCCTTCATGGCGTGGAGGATGTCAATATCGGTTGCTGTCAGCGATTAAAATCAACCCCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGACGTTACCAGAAAt                <  8:250144/149‑1 (MQ=255)
                                                                                                                                 tCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATg               >  6:356609/1‑149 (MQ=255)
                                                                                                                                   cACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGACCCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGc              >  6:262941/1‑148 (MQ=255)
                                                                                                                                   cACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGc              <  1:26190/148‑1 (MQ=255)
                                                                                                                                    aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTg                              <  3:61763/131‑1 (MQ=255)
                                                                                                                                    aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTg                              >  4:61763/1‑131 (MQ=255)
                                                                                                                                    aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca            <  3:454974/149‑1 (MQ=255)
                                                                                                                                    aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca            >  1:240142/1‑149 (MQ=255)
                                                                                                                                         cTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAACCGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCaccacc       <  4:146159/149‑1 (MQ=255)
                                                                                                                                          ttCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                          <  1:191447/49‑1 (MQ=255)
                                                                                                                                          ttCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                          >  2:191447/1‑49 (MQ=255)
                                                                                                                                          ttCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCAccg      >  3:442762/1‑149 (MQ=255)
                                                                                                                                          ttCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCAccg      >  4:300429/1‑149 (MQ=255)
                                                                                                                                              gggCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcc             >  6:89664/1‑138 (MQ=255)
                                                                                                                                              gggCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcc             <  5:89664/138‑1 (MQ=255)
                                                                                                                                              gggCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGGCAATCCACTTGCCGTTACCAGAAATGCCACCACcgccgc  <  5:262941/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGC  >  NC_000913/1504749‑1505039

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: