Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,799,634 A→G 10.7% intergenic (+21/+16) 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,6340AG10.7% 238.4 / 11.6 81intergenic (+21/+16)alaE/ygaCalanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (7/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATC  >  NC_000913/2799487‑2799778
                                                                                                                                                    |                                                                                                                                                
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                                                                                                                                        >  4:591275/1‑141 (MQ=255)
                 tCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑GATGCCAGTcg                                                                                                                                        <  3:591275/141‑1 (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                                                                                                                                         <  2:814286/121‑1 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTc                                                                                                                                         >  1:814286/1‑121 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                            <  2:47184/149‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                            <  4:916471/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:563800/1‑149 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                  >  2:254778/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                          <  4:693303/149‑1 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                          >  2:94901/1‑149 (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               >  1:327051/1‑149 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt               <  2:263649/149‑1 (MQ=255)
                                                                                                                                     ttttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt            <  3:84893/149‑1 (MQ=255)
                                                                                                                                       ttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg          >  1:822453/1‑149 (MQ=255)
                                                                                                                                       ttCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg          >  1:920901/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)
                                                                                                                                                    |                                                                                                                                                
TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGC‑ATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATC  >  NC_000913/2799487‑2799778

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: ‑