Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,799,635 T→A 9.3% intergenic (+22/+15) alaE → / ← ygaC alanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,799,6350TA 9.3% 275.0 / 15.7 87intergenic (+22/+15)alaE/ygaCalanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base A (5/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAA  >  NC_000913/2799487‑2799782
                                                                                                                                                     |                                                                                                                                                   
tCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCAt                                                                                                                                                     >  1:441957/1‑149 (MQ=255)
tCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCAt                                                                                                                                                     >  1:53626/1‑149 (MQ=255)
    aTTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGcc                                                                                                                                                 >  3:740926/1‑149 (MQ=255)
       aTGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGt                                                                                                                                              >  1:265980/1‑149 (MQ=255)
        tGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTc                                                                                                                                             <  1:800015/149‑1 (MQ=255)
        tGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                             <  2:657424/149‑1 (MQ=255)
        tGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                             <  2:743924/149‑1 (MQ=255)
                   aGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcgc                                                                                                                                           <  2:533558/140‑1 (MQ=255)
                   aGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcgc                                                                                                                                           >  1:533558/1‑140 (MQ=255)
                    gTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAGGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACt                                                                                                                                 <  4:213319/149‑1 (MQ=255)
                     ttCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTg                                                                                                                                >  4:200634/1‑149 (MQ=255)
                        aaaCATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCt                                                                                                                             >  4:875554/1‑149 (MQ=255)
                                ttgttTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCaaa                                                                                                                     <  2:776348/149‑1 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATaaa                                                                                                                 <  1:549861/149‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                                <  2:258394/149‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                                <  2:265980/149‑1 (MQ=255)
                                        aTGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAACAAACCAt                                                                                                             >  4:810715/1‑149 (MQ=255)
                                          gTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTc                                                                                                           >  1:264319/1‑149 (MQ=255)
                                                 ggggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGCCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACCGGCTTCGCCAAATACACCACTCAAATAAc                                                                                                    <  4:842907/149‑1 (MQ=255)
                                                 ggggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAAc                                                                                                    <  2:404003/149‑1 (MQ=255)
                                                   ggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGt                                                                                                  >  3:198078/1‑149 (MQ=255)
                                                    gCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGtt                                                                                                 <  3:563630/149‑1 (MQ=255)
                                                     cGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTc                                                                                                <  2:725640/149‑1 (MQ=255)
                                                      ggTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGCCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCa                                                                                               <  3:875554/149‑1 (MQ=255)
                                                       gTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                              <  1:221816/149‑1 (MQ=255)
                                                       gTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                              >  2:221816/1‑149 (MQ=255)
                                                             ggCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCaa                                                                                        >  4:719829/1‑149 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGt                                                                                                  <  2:831148/136‑1 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGt                                                                                                  >  1:831148/1‑136 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAacac                                                                                     >  3:83207/1‑149 (MQ=255)
                                                                 aCTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACAcg                                                                                    >  4:198209/1‑149 (MQ=255)
                                                                     ccTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGaa                                                                                >  3:132108/1‑149 (MQ=255)
                                                                      cTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAAc                                                                               >  4:200618/1‑149 (MQ=255)
                                                                       tCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACt                                                                              >  1:771287/1‑149 (MQ=255)
                                                                       tCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCAGCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACt                                                                              >  2:637891/1‑149 (MQ=255)
                                                                          attattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaaa                                                                           <  3:369223/149‑1 (MQ=255)
                                                                           ttattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAACGCAACACGCGAACTaaa                                                                           <  1:712679/148‑1 (MQ=255)
                                                                           ttattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCCCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaaa                                                                           >  2:712679/1‑148 (MQ=255)
                                                                                 ccgccgACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGc                                                                    <  3:121331/149‑1 (MQ=255)
                                                                                  cgccgACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCg                                                                   >  1:40294/1‑149 (MQ=255)
                                                                                       aCTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATa                                                              >  3:604709/1‑149 (MQ=255)
                                                                                           tttAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGt                                                          >  4:532926/1‑149 (MQ=255)
                                                                                              aaaGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAAt                                                       <  3:381881/149‑1 (MQ=255)
                                                                                               aaGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                              >  3:816448/1‑109 (MQ=255)
                                                                                               aaGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                              <  4:816448/109‑1 (MQ=255)
                                                                                                  tCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTaa                                                        >  4:780181/1‑144 (MQ=255)
                                                                                                  tCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTaa                                                        <  3:780181/144‑1 (MQ=255)
                                                                                                    aGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaaa                                                                           <  3:770635/123‑1 (MQ=255)
                                                                                                    aGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaaa                                                                           >  4:770635/1‑123 (MQ=255)
                                                                                                    aGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCgt                                                 <  2:441957/149‑1 (MQ=255)
                                                                                                    aGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCgt                                                 <  3:693040/149‑1 (MQ=255)
                                                                                                        gTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcg                                             >  2:1480/1‑149 (MQ=255)
                                                                                                         ttACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgc                                            <  4:746798/149‑1 (MQ=255)
                                                                                                         ttACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgc                                            <  1:317358/149‑1 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCAtt                                                                                                            <  3:276794/79‑1 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCAtt                                                                                                            >  4:276794/1‑79 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                      >  3:801546/1‑149 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                      >  2:217534/1‑149 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                      >  4:625600/1‑149 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCt                                   >  4:460600/1‑149 (MQ=255)
                                                                                                                   gTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAAt                                                       <  1:609646/128‑1 (MQ=255)
                                                                                                                   gTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAAt                                                       >  2:609646/1‑128 (MQ=255)
                                                                                                                   gTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTa                                  >  3:744610/1‑149 (MQ=255)
                                                                                                                     aaaaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTaaa                                <  4:694796/149‑1 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                              <  2:301877/149‑1 (MQ=255)
                                                                                                                          ccTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATagag                           <  4:29783/149‑1 (MQ=255)
                                                                                                                           cTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGc                          <  4:198078/149‑1 (MQ=255)
                                                                                                                              aCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGa                       >  1:637798/1‑149 (MQ=255)
                                                                                                                              aCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGa                       >  2:112586/1‑149 (MQ=255)
                                                                                                                                tcttctTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGaaa                     >  4:454541/1‑149 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt                   <  4:868713/149‑1 (MQ=255)
                                                                                                                                   tcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATACCGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTg                  <  2:638087/149‑1 (MQ=255)
                                                                                                                                   tcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTg                  <  4:515747/149‑1 (MQ=255)
                                                                                                                                    cttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc                 <  1:223890/149‑1 (MQ=255)
                                                                                                                                     ttttCGCGTCTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt                <  2:771287/149‑1 (MQ=255)
                                                                                                                                     ttttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt                >  2:159198/1‑149 (MQ=255)
                                                                                                                                     ttttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt                >  1:502512/1‑149 (MQ=255)
                                                                                                                                     aagacgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt                >  3:684538/5‑149 (MQ=255)
                                                                                                                                       ttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg              <  3:200618/149‑1 (MQ=255)
                                                                                                                                       gacgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCACAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc                 >  2:847781/3‑146 (MQ=255)
                                                                                                                                       gacgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCACAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc                 <  1:847781/144‑1 (MQ=255)
                                                                                                                                          gcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgct                                           >  2:793667/1‑117 (MQ=255)
                                                                                                                                          gcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgct                                           <  1:793667/117‑1 (MQ=255)
                                                                                                                                          gcgACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa           >  1:55463/1‑149 (MQ=255)
                                                                                                                                            gACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTTCTGGTAATg         >  1:294616/1‑149 (MQ=255)
                                                                                                                                               tGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc      >  4:227222/1‑149 (MQ=255)
                                                                                                                                               tGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc      <  3:719829/149‑1 (MQ=255)
                                                                                                                                                ggCATCGCCAGTGGCGTCTTTAACTGGCTTCGTCAAATAAACCATTTAAATAAGATTCAAGAGCAACACGCGAAATAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTGATGATCg     >  4:836505/1‑149 (MQ=255)
                                                                                                                                                    aTCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCaa  >  2:895943/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAA  >  NC_000913/2799487‑2799782

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