Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,799,635 T→A 10.0% 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,6350TA10.0% 256.1 / 15.5 80intergenic (+22/+15)alaE/ygaCalanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base A (6/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAG  >  NC_000913/2799487‑2799783
                                                                                                                                                     |                                                                                                                                                    
tCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑At                                                                                                                                                      >  1:223456/1‑149 (MQ=255)
   gATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGc                                                                                                                                                   >  1:546693/1‑149 (MQ=255)
     ttATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCa                                                                                                                                                 >  3:838157/1‑149 (MQ=255)
        tGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTc                                                                                                                                              <  2:149449/149‑1 (MQ=255)
          gcggcgGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTcgc                                                                                                                                            >  3:181677/1‑149 (MQ=255)
            ggcggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGt                                                                                                                                          <  2:675716/149‑1 (MQ=255)
             gcggTCAGTTCAAACATCGGTGTGTCGATGTTGATGGGGGCGGTTTAGGGCTACTTCCTCGATTATTGCCTCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTACAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTc                                                                                                                                         <  3:16974/149‑1 (MQ=255)
               ggTCGGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       >  4:539583/1‑149 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       <  2:7783/149‑1 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       <  3:146508/149‑1 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       <  4:321387/149‑1 (MQ=255)
                 tCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑GATGCCAGTcg                                                                                                                                             <  3:591275/141‑1 (MQ=255)
                 tCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑GATGCCAGTcg                                                                                                                                             >  4:591275/1‑141 (MQ=255)
                  cAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTAAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTaa                                                                                                                                    <  4:63354/149‑1 (MQ=255)
                   aGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAAc                                                                                                                                   <  3:883232/149‑1 (MQ=255)
                     ttCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATTGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTTAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑GATGCCAGTctc                                                                                                                                            >  1:911068/1‑136 (MQ=255)
                     ttCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGAT‑GCCAGTcgc                                                                                                                                            <  2:911068/138‑1 (MQ=255)
                         aaCATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCtt                                                                                                                             >  1:320360/1‑149 (MQ=255)
                          aCAGCGGTGTTTCGATGTTGATGGGGGCGGTTTTTGGCTACTGCCTCGCTGCGGGGCGCCGACTGTTTTCAGGCTGTGGTTGCGAGCAGGCAACAGCATGGGGCTGTTTTTCGCGACGGGC‑GTCACCAGTTGCTCCTTTGACTGGCTTc                                                                                                                            >  4:676535/1‑149 (MQ=255)
                                ttgttTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGACGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCaaa                                                                                                                      >  2:567553/1‑149 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTc                                                                                                                                              >  1:814286/1‑121 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTc                                                                                                                                              <  2:814286/121‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                                 <  4:916471/149‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                                 <  2:47184/149‑1 (MQ=255)
                                      cGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGt                                                                                                                                               >  2:146123/1‑118 (MQ=255)
                                      cGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGt                                                                                                                                               <  1:146123/118‑1 (MQ=255)
                                       gATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAAGTTGCTTCGCCAAATAAACCa                                                                                                               >  3:760436/1‑149 (MQ=255)
                                       gATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCa                                                                                                               >  4:95268/1‑149 (MQ=255)
                                             gATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCaaa                                                                                                         <  2:21133/149‑1 (MQ=255)
                                                 ggggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAAc                                                                                                     <  2:543396/149‑1 (MQ=255)
                                                  gggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       >  1:391683/1‑114 (MQ=255)
                                                  gggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCtt                                                                                                                                       <  2:391683/114‑1 (MQ=255)
                                                    gCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGtt                                                                                                  >  3:354798/1‑149 (MQ=255)
                                                     cGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTc                                                                                                 <  4:762755/149‑1 (MQ=255)
                                                       gTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                               >  3:693303/1‑149 (MQ=255)
                                                           aTGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAgcgc                                                                                           <  1:567553/149‑1 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAacac                                                                                      >  1:807614/1‑149 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAacac                                                                                      >  4:453360/1‑149 (MQ=255)
                                                                    tCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGa                                                                                  <  1:848606/149‑1 (MQ=255)
                                                                      cTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAAc                                                                                >  4:217034/1‑149 (MQ=255)
                                                                       tCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACt                                                                               >  2:495867/1‑149 (MQ=255)
                                                                         gATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaa                                                                             <  4:181677/149‑1 (MQ=255)
                                                                           ttattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAg                                                                           <  2:320360/149‑1 (MQ=255)
                                                                                         tGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAcc                                                             >  4:35163/1‑149 (MQ=255)
                                                                                           tttAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGt                                                           >  4:84893/1‑149 (MQ=255)
                                                                                             tAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCgg                                                                   <  2:772653/139‑1 (MQ=255)
                                                                                             tAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCgg                                                                   >  1:772653/1‑139 (MQ=255)
                                                                                              aaaGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAAt                                                        >  2:250127/1‑149 (MQ=255)
                                                                                                  tCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTc                                                    <  1:495867/149‑1 (MQ=255)
                                                                                                              agcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGcc                                        <  3:217034/149‑1 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                       >  2:254778/1‑149 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                       >  2:563800/1‑149 (MQ=255)
                                                                                                                 agGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGc                                     <  3:539583/149‑1 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATaaa                                                                                                                  <  1:80591/71‑1 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATaaa                                                                                                                  >  2:80591/1‑71 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCt                                    <  4:354798/149‑1 (MQ=255)
                                                                                                                    tAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTaa                                  >  1:263649/1‑149 (MQ=255)
                                                                                                                     aaaaGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTaaa                                 <  2:546693/149‑1 (MQ=255)
                                                                                                                      aaaGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAAt                                >  4:847593/1‑149 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                               >  2:94901/1‑149 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                               <  4:693303/149‑1 (MQ=255)
                                                                                                                          ccTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATagag                            <  1:250127/149‑1 (MQ=255)
                                                                                                                             gACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCgg                         >  2:231175/1‑149 (MQ=255)
                                                                                                                              aCTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGa                        <  2:202787/149‑1 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt                    <  2:263649/149‑1 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt                    >  1:327051/1‑149 (MQ=255)
                                                                                                                                     ttttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt                 <  3:84893/149‑1 (MQ=255)
                                                                                                                                       ttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg               >  1:920901/1‑149 (MQ=255)
                                                                                                                                       ttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg               >  1:822453/1‑149 (MQ=255)
                                                                                                                                        acgcgACTGGC‑GATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAc                                                              >  4:238042/2‑101 (MQ=255)
                                                                                                                                        acgcgACTGGC‑GATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAc                                                              <  3:238042/100‑1 (MQ=255)
                                                                                                                                          gcgACTGGC‑GATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa            >  3:623013/1‑149 (MQ=255)
                                                                                                                                          gcgACTGGC‑GATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa            >  1:450140/1‑149 (MQ=255)
                                                                                                                                            gACTGGC‑GATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATg          >  1:519339/1‑149 (MQ=255)
                                                                                                                                            gACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATg          >  4:758273/1‑149 (MQ=255)
                                                                                                                                             aCTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGa         <  1:518857/149‑1 (MQ=255)
                                                                                                                                             aCTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGa         <  1:254778/149‑1 (MQ=255)
                                                                                                                                               tGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc       >  3:303739/1‑149 (MQ=255)
                                                                                                                                               tGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc       <  2:223456/149‑1 (MQ=255)
                                                                                                                                                 gC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc                  >  1:781203/1‑136 (MQ=255)
                                                                                                                                                 gC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc                  <  2:781203/136‑1 (MQ=255)
                                                                                                                                                     tCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAg  >  2:936532/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAG  >  NC_000913/2799487‑2799783

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