Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,308,355 G→T 22.1% L42M (CTG→ATG)  alsE ← allulose‑6‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,308,3550GT22.1% 139.5 / 34.5 59L42M (CTG→ATG) alsEallulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base T (8/5);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

CATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGC  >  NC_000913/4308074‑4308637
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                          
cagtaacGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTTGCGTCTCCGTGTTAAGGATCAGCCTCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTggg                                                                                                                                                                                                                                                                                      <  1:534349/282‑1 (MQ=255)
cATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTggg                                                                                                                                                                                                                                                                                      <  1:440503/288‑1 (MQ=255)
    aCCGTAATTTTTTCTGCCTTTTGGATATAGTATTTCCTGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACaa                                                                                                                                                                                                                                                                                 <  2:333798/289‑1 (MQ=255)
         aaTTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGa                                                                                                                                                                                                                                                                            <  1:301564/289‑1 (MQ=255)
           ttttACGGGCCCTATGACTATATTATGTCTTGGCCTCAACTGGCGTCTCCGGGTTAAAGATCGGCCCCACTATCATGTCCTGACGGCGGATTTGATCAATCAGGCGGACCGCCTGGCCGTTGATGGTACCCGGATGGAGAGTGATGAAATCTGCTCCCGCACGCGCCCGTTGGGCAATGTAAGCCTGCGGCCGCGTCACCATCAGATGGCAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGAAAAAGTGAcc                                                                                                                                                                                                                                                                          <  2:374235/289‑1 (MQ=255)
              tATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTc                                                                                                                                                                                                                                                                       <  1:248544/289‑1 (MQ=255)
                              ataGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTgg                                                                                                                                                                                                                                                       >  1:256602/1‑289 (MQ=255)
                              ataGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTgg                                                                                                                                                                                                                                                       >  2:405144/1‑289 (MQ=255)
                              ataGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTgg                                                                                                                                                                                                                                                       >  1:97635/1‑289 (MQ=255)
                                aGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGaa                                                                                                                                                                                                                                                     >  2:342419/1‑289 (MQ=255)
                                  tATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGACATCGATGTGGAAGt                                                                                                                                                                                                                                                   >  1:69876/1‑289 (MQ=255)
                                     ttCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTaa                                                                                                                                                                                                                                                 >  2:287469/1‑288 (MQ=255)
                                      tCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATc                                                                                                                                                                                                                                               >  1:525175/1‑289 (MQ=255)
                                                aaCTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATg                                                                                                                                                                                                                                         >  1:150045/1‑285 (MQ=255)
                                                                         aGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTaa                                                                                                                                                                                                            >  1:233020/1‑289 (MQ=255)
                                                                         aGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTCACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGAACGTCCATGATATCGATGTGGACGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTaa                                                                                                                                                                                                            >  1:155711/1‑289 (MQ=255)
                                                                          gCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATATGTTCTTTaaa                                                                                                                                                                                                           >  2:84675/1‑289 (MQ=255)
                                                                                          atgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATa                                                                                                                                                                                           <  2:150045/289‑1 (MQ=255)
                                                                                           tgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATac                                                                                                                                                                                          <  1:287469/289‑1 (MQ=255)
                                                                                               ggcggATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATAcaca                                                                                                                                                                                       <  2:256602/288‑1 (MQ=255)
                                                                                                 cggATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGCTAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTa                                                                                                                                                                                    >  1:293446/1‑289 (MQ=255)
                                                                                                    aTTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAAc                                                                                                                                                                                  <  1:342419/288‑1 (MQ=255)
                                                                                                     tttCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACg                                                                                                                                                                                 <  1:143190/288‑1 (MQ=255)
                                                                                                     ttatctcaagcaGGCGGAGCGCCTGGCGGTTGATGGTTTCCGGATGCAGAGTGATGAATTCTGCTCCCGCACGGGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGGGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACg                                                                                                                                                                                 <  2:69876/283‑1 (MQ=255)
                                                                                                          tcaatcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAggggg                                                                                                                                                                           <  2:233020/289‑1 (MQ=255)
                                                                                                                   cGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTc                                                                                                                                                                   <  1:405144/288‑1 (MQ=255)
                                                                                                                                   tggtGGTTTCCGGATGCATATTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCcttact                                                                                                                                                  <  2:320287/286‑1 (MQ=255)
                                                                                                                                   tGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCcttact                                                                                                                                                  <  1:196504/289‑1 (MQ=255)
                                                                                                                                   tGATGCATTCCGGATGCAGAGTGATAAAAGCTTCACCCGAACGCGCCATTTGGGCAAAGTAATACTGCGGCCGCGTCACCATCAGATGACAGTCGGGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGAGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTTTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCcttact                                                                                                                                                  <  2:293446/289‑1 (MQ=255)
                                                                                                                                                          atgaAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCg                                                                                                                           >  2:359725/1‑289 (MQ=255)
                                                                                                                                                                tCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGcc                                                                                                                     >  1:152886/1‑289 (MQ=255)
                                                                                                                                                                tCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGcc                                                                                                                     >  2:361218/1‑289 (MQ=255)
                                                                                                                                                                  tGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTAATGTTCAATTCATTACTTACTGATAAGACGGTCAAGGGCGACAGcag                                                                                                                    >  2:332033/1‑286 (MQ=255)
                                                                                                                                                                              cgcCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATc                                                                                                        >  2:218424/1‑288 (MQ=255)
                                                                                                                                                                                     tGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGccc                                                                                                   <  1:359725/286‑1 (MQ=255)
                                                                                                                                                                                        gCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATcacc                                                                                             >  2:282344/1‑289 (MQ=255)
                                                                                                                                                                                                      cGGCCGCGTCACCATCGAAGGCAAGTGGGGAGGTTTAGTTGCCAGTTTTTTAACCGGACTTACGAGGAACAGTGAGAGTGGCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAAt                                                                                <  2:152886/288‑1 (MQ=255)
                                                                                                                                                                                                      cGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTcagca                                                                                                                                                                                                   >  2:289857/1‑173 (MQ=255)
                                                                                                                                                                                                      cGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTcagca                                                                                                                                                                                                   <  1:289857/173‑1 (MQ=255)
                                                                                                                                                                                                      cGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATa                                                                               >  1:234375/1‑289 (MQ=255)
                                                                                                                                                                                                                             aGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCcttactt                                                                                                                                                 >  2:70197/1‑200 (MQ=255)
                                                                                                                                                                                                                             aGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCcttactt                                                                                                                                                 <  1:70197/200‑1 (MQ=255)
                                                                                                                                                                                                                               tCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAg                                                       >  1:238736/1‑288 (MQ=255)
                                                                                                                                                                                                                                  aGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCg                                                   >  1:3480/1‑289 (MQ=255)
                                                                                                                                                                                                                                       tttAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCgttgtt                                              >  2:259084/1‑289 (MQ=255)
                                                                                                                                                                                                                                         tAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTcagca                                                                                                                                                                                                   >  2:131992/1‑138 (MQ=255)
                                                                                                                                                                                                                                         tAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTcagca                                                                                                                                                                                                   <  1:131992/138‑1 (MQ=255)
                                                                                                                                                                                                                                         gaGTTGCCAGGTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGa                                            <  2:3480/288‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCCAGTTTTGTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAATGGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATg                                          <  2:487023/287‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACgg                                                                                                                                    >  1:527656/1‑193 (MQ=255)
                                                                                                                                                                                                                                                 aGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACgg                                                                                                                                    <  2:527656/193‑1 (MQ=255)
                                                                                                                                                                                                                                                    tttttAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCg                                    <  1:282344/286‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAAcc                          <  1:361218/289‑1 (MQ=255)
                                                                                                                                                                                                                                                               aCTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAAt                                                                                <  2:162723/231‑1 (MQ=255)
                                                                                                                                                                                                                                                               aCTTACGAAGAACGGTGAGAGTGTCATATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAAt                                                                                >  1:162723/1‑231 (MQ=255)
                                                                                                                                                                                                                                                                cTTACGAGGAACGGGGAGAGTGTAAGATTGGGGACAAAGTGACCGTCCATGAGTTCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGGTCGATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGcc                      <  2:297225/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaCGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACcaca            >  2:379557/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                           cGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACaga          >  2:547895/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   gtgtCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATgc  >  1:106778/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                          
CATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACCGTCCATGATATCGATGTGGAAGTAATCGGCATGGCTGTCGATAAATTCGATCTGTTCTTTAAATTTCAGCAGATCCATACACATTAACGAGGGGGAGATTTTCATGTTCAATTCCTTACTTACTGATAAGACGGTCAAGGGCGACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGC  >  NC_000913/4308074‑4308637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: