Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,037:1 +T 9.8% coding (282/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,0371.T 9.8% 290.0 / 29.9 82coding (282/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base . (25/49);  new base T (3/5);  total (28/54)
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.17e-01

GCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTA  >  NC_000913/2324889‑2325186
                                                                                                                                                     |                                                                                                                                                     
gCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCg                                                                                                                                                      <  1:518415/150‑1 (MQ=255)
gCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCg                                                                                                                                                      >  2:677345/1‑150 (MQ=255)
 cATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGccctg                                                                                                                                                      >  2:1324713/1‑148 (MQ=255)
     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:1304819/1‑150 (MQ=255)
                 ggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTg                                                                                                                                     <  1:1088533/150‑1 (MQ=255)
                 ggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTg                                                                                                                                     <  2:482765/150‑1 (MQ=255)
                  gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                    <  1:115115/150‑1 (MQ=255)
                  gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTgg                                                                                                                                    <  2:81631/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                                                                                             >  1:785480/1‑150 (MQ=255)
                                                         tGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTggg                                                                                             <  2:1530073/150‑1 (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:89325/1‑150 (MQ=255)
                                                                                         gCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGcc                                                             <  2:220082/150‑1 (MQ=255)
                                                                                             gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTg                                                         <  1:1355426/150‑1 (MQ=255)
                                                                                             gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTg                                                         >  2:759029/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                     
GCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCGTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTA  >  NC_000913/2324889‑2325186

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: