Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,042 Δ1 bp 9.2% coding (287/1323 nt) atoE → short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,325,0420A. 9.2% 275.5 / 20.6 87coding (287/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base A (28/51);  new base . (3/5);  total (31/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 = 4.70e-01

GTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTACATT  >  NC_000913/2324894‑2325190
                                                                                                                                                     |                                                                                                                                                    
gTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAg                                                                                                                                                     >  2:713495/1‑150 (MQ=255)
  gAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAggg                                                                                                                                                   <  1:626996/150‑1 (MQ=255)
    aaaTGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCg                                                                                                                                                 <  1:135689/150‑1 (MQ=255)
     aaTGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGt                                                                                                                                                <  2:1576789/150‑1 (MQ=255)
       tgttGGGTGACGGTTTCTTTAACTTGCTGGCGTTTGGTATGCAGATTGCTCTTATCATCGTTTCCGGTCTTGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTTCGTACTGCCGACTGCGCCGCAAAGACGCCCGTACAGGGCGTCa                                                                                                                                              <  2:806753/150‑1 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                          <  1:942395/150‑6 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                          <  1:1324713/150‑6 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                          <  2:644109/150‑6 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                          <  1:1126052/150‑6 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCt                                                                                                                                          <  2:310913/150‑1 (MQ=255)
           gggTGACGGTTTCTGGAACTTGCTGGCGTATGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCGCCGCCGCAAAGACGCCCTGTACGGGCGTctttgc                                                                                                                                          <  2:262601/150‑6 (MQ=255)
           gccttaCGGTTTCTCGAACTTGCTGGCGTTTGCTATCCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCt                                                                                                                                          >  2:1355426/6‑150 (MQ=255)
            ggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTg                                                                                                                                         <  2:482765/150‑1 (MQ=255)
            ggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTg                                                                                                                                         <  1:1088533/150‑1 (MQ=255)
            ggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTg                                                                                                                                         >  2:1304819/1‑150 (MQ=255)
             gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                        <  2:81631/150‑1 (MQ=255)
             gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                        <  1:115115/150‑1 (MQ=255)
             gTGACGGTGTCTTGAACTTGCTGGCGTTTGGTATGGAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                        <  2:838476/150‑1 (MQ=255)
                 cGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCAc                                                                                                                                    <  2:1190098/150‑1 (MQ=255)
                     ttCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCGGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTttc                                                                                                                                <  2:1290524/150‑1 (MQ=255)
                       cTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTttctt                                                                                                                              >  1:1518055/1‑150 (MQ=255)
                        tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                             <  1:62030/150‑1 (MQ=255)
                         ggAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCg                                                                                                                            <  1:677345/150‑1 (MQ=255)
                          gAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCgg                                                                                                                           <  1:487957/150‑1 (MQ=255)
                          gAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCgg                                                                                                                           >  1:1328915/1‑150 (MQ=255)
                           tagttGCTGGCGTTTGGTATGCAGATGGCGCTTTTCATCGTTACCGGTCATGCCCTTGCCTGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGt                                                                                                                          <  1:980769/147‑1 (MQ=255)
                                 cTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTc                                                                                                                    <  2:1476863/150‑1 (MQ=255)
                                  tGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCg                                                                                                                   >  1:220082/1‑150 (MQ=255)
                                   ggCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGc                                                                                                                  >  2:116628/1‑150 (MQ=255)
                                      atttGGTATGCAGATGGCGCTTATCATCGTTACCGGTTATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTtg                                                                                                               <  2:653040/149‑1 (MQ=255)
                                        ttGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTtgtg                                                                                                             <  2:1127173/150‑1 (MQ=255)
                                         tGGTAGGCAGATGGCGATTATCATGGTTAGCGGTCATGCCCTTGCCAGCTCTGCTGAGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTgtgt                                                                                                            <  2:1161924/150‑1 (MQ=255)
                                          ggTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGAGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtc                                                                                                           >  1:549246/1‑150 (MQ=255)
                                           gTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtca                                                                                                          <  1:21085/150‑1 (MQ=255)
                                             aTGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatc                                                                                                        <  1:270077/150‑1 (MQ=255)
                                                  gATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTg                                                                                                   >  2:1248947/1‑150 (MQ=255)
                                                   aTGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCTGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCTGTCGCTTGTGTCATCAACTgg                                                                                                  <  2:935783/150‑1 (MQ=255)
                                                    tGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTggg                                                                                                 <  2:1530073/150‑1 (MQ=255)
                                                    tGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTggg                                                                                                 >  1:785480/1‑150 (MQ=255)
                                                           tATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttgg                                                                                          >  1:804216/1‑150 (MQ=255)
                                                           tATCATCGTTACCGGGCATGCCCTTGCCAGCTCTGCTCCGGTGAAACGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGTAtttgg                                                                                          <  2:94926/150‑1 (MQ=255)
                                                            atcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggt                                                                                         >  2:1436476/1‑150 (MQ=255)
                                                              catcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                       >  2:1600264/1‑150 (MQ=255)
                                                                 cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATttggtttggt                                                                                    >  1:1128271/1‑150 (MQ=255)
                                                                     aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                                <  1:1324126/150‑1 (MQ=255)
                                                                     aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                                <  2:182427/150‑1 (MQ=255)
                                                                       cGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCgg                                                                              <  1:1248947/150‑1 (MQ=255)
                                                                          tCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgcgc                                                                           <  2:785480/150‑1 (MQ=255)
                                                                           cATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCa                                                                          <  1:936093/150‑1 (MQ=255)
                                                                              gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATg                                                                       <  1:1004324/150‑1 (MQ=255)
                                                                                   tGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCATACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGc                                                                  >  1:308874/1‑150 (MQ=255)
                                                                                    gCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGcc                                                                 <  2:220082/150‑1 (MQ=255)
                                                                                    gCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGcc                                                                 >  2:89325/1‑150 (MQ=255)
                                                                                        gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTg                                                             >  2:759029/1‑150 (MQ=255)
                                                                                        gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTg                                                             <  1:1355426/150‑1 (MQ=255)
                                                                                               gccGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGGCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGc                                                      <  1:315206/149‑1 (MQ=255)
                                                                                                   gTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCgg                                                  >  2:1208822/1‑150 (MQ=255)
                                                                                                    ggAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGc                                                 <  2:844520/149‑1 (MQ=255)
                                                                                                     gAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCg                                                <  1:1263996/150‑1 (MQ=255)
                                                                                                       aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                              <  2:550046/150‑1 (MQ=255)
                                                                                                          gTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTcc                                           <  1:1242539/150‑1 (MQ=255)
                                                                                                            ttGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTcccc                                         <  1:564652/150‑1 (MQ=255)
                                                                                                            ttGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGGTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTGTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTcccc                                         <  1:559174/150‑1 (MQ=255)
                                                                                                                 gcGTACTGCCGCCTCCGCCTCAAAGACGCCCGGACAGGGCGTCATGCTGGGCTCTTTCTTCGGTTCAGGCGCTTGTGGCATCAACAGGGGATTTTGTGTGGTTGGCTTCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTc                                    <  2:1119819/150‑1 (MQ=255)
                                                                                                                 gcGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTc                                    >  1:411921/1‑150 (MQ=255)
                                                                                                                     aCTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGAt                                <  1:495026/150‑1 (MQ=255)
                                                                                                                     aCTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGAt                                <  1:859386/150‑1 (MQ=255)
                                                                                                                     aCTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGAt                                >  2:502628/1‑150 (MQ=255)
                                                                                                                          cgccTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATcc                           <  1:1600264/150‑1 (MQ=255)
                                                                                                                             cTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt                        <  2:596568/150‑1 (MQ=255)
                                                                                                                             cTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt                        <  2:602146/150‑1 (MQ=255)
                                                                                                                               ccgccgCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTTAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGc                      >  1:1139529/1‑150 (MQ=255)
                                                                                                                                  ccgcAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCa                   >  2:87628/1‑150 (MQ=255)
                                                                                                                                   cagcatgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCAt                  >  1:1271663/7‑150 (MQ=255)
                                                                                                                                     gcatgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTg                >  1:994171/5‑150 (MQ=255)
                                                                                                                                      catgACGCCCTGTACGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc               >  1:670070/4‑150 (MQ=255)
                                                                                                                                      aaaGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGTTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc               <  2:1589614/150‑1 (MQ=255)
                                                                                                                                       aaGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc              <  2:1328915/150‑1 (MQ=255)
                                                                                                                                        aGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCt             <  1:782223/150‑1 (MQ=255)
                                                                                                                                           cGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg          >  2:511990/1‑150 (MQ=255)
                                                                                                                                           cGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg          >  2:128303/1‑150 (MQ=255)
                                                                                                                                            gCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcgc         >  1:2404/1‑150 (MQ=255)
                                                                                                                                             cccGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc        <  1:87628/150‑1 (MQ=255)
                                                                                                                                             cccGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc        >  2:1559712/1‑150 (MQ=255)
                                                                                                                                               cGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTa      <  2:411921/150‑1 (MQ=255)
                                                                                                                                                   aCAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCAGTTGCTCATTGCCTGCGCCTACAt   >  2:1233573/1‑150 (MQ=255)
                                                                                                                                                    cAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTACAtt  >  1:1510977/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTACATT  >  NC_000913/2324894‑2325190

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