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

TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATC  >  NC_000913/2799487‑2799778
                                                                                                                                                    |                                                                                                                                                
tCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCAt                                                                                                                                                 >  1:223456/1‑149 (MQ=255)
   gATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGc                                                                                                                                              >  1:546693/1‑149 (MQ=255)
     ttATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCa                                                                                                                                            >  3:838157/1‑149 (MQ=255)
        tGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                         <  2:149449/149‑1 (MQ=255)
          gcggcgGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTcgc                                                                                                                                       >  3:181677/1‑149 (MQ=255)
            ggcggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGt                                                                                                                                     <  2:675716/149‑1 (MQ=255)
             gcggTCAGTTCAAACATCGGTGTGTCGATGTTGATGGGGGCGGTTTAGGGCTACTTCCTCGATTATTGCCTCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTACAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTc                                                                                                                                    <  3:16974/149‑1 (MQ=255)
               ggTCGGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  >  4:539583/1‑149 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  <  2:7783/149‑1 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  <  3:146508/149‑1 (MQ=255)
               ggTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  <  4:321387/149‑1 (MQ=255)
                 tCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcg                                                                                                                                        >  4:591275/1‑141 (MQ=255)
                 tCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcg                                                                                                                                        <  3:591275/141‑1 (MQ=255)
                  cAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTAAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTaa                                                                                                                               <  4:63354/149‑1 (MQ=255)
                   aGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAAc                                                                                                                              <  3:883232/149‑1 (MQ=255)
                     ttCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATTGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTTAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTctc                                                                                                                                       >  1:911068/1‑136 (MQ=255)
                     ttCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcgc                                                                                                                                       <  2:911068/138‑1 (MQ=255)
                         aaCATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCtt                                                                                                                        >  1:320360/1‑149 (MQ=255)
                          aCAGCGGTGTTTCGATGTTGATGGGGGCGGTTTTTGGCTACTGCCTCGCTGCGGGGCGCCGACTGTTTTCAGGCTGTGGTTGCGAGCAGGCAACAGCATGGGGCTGTTTTTCGCGACGGGCGTCACCAGTTGCTCCTTTGACTGGCTTc                                                                                                                       >  4:676535/1‑149 (MQ=255)
                                ttgttTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGACGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCaaa                                                                                                                 >  2:567553/1‑149 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                         <  2:814286/121‑1 (MQ=255)
                                    ttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                         >  1:814286/1‑121 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                            <  2:47184/149‑1 (MQ=255)
                                     tCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAAc                                                                                                            <  4:916471/149‑1 (MQ=255)
                                      cGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGt                                                                                                                                          >  2:146123/1‑118 (MQ=255)
                                      cGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGt                                                                                                                                          <  1:146123/118‑1 (MQ=255)
                                       gATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAAGTTGCTTCGCCAAATAAACCa                                                                                                          >  3:760436/1‑149 (MQ=255)
                                       gATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCa                                                                                                          >  4:95268/1‑149 (MQ=255)
                                             gATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCaaa                                                                                                    <  2:21133/149‑1 (MQ=255)
                                                 ggggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAAc                                                                                                <  2:543396/149‑1 (MQ=255)
                                                  gggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  >  1:391683/1‑114 (MQ=255)
                                                  gggCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCtt                                                                                                                                  <  2:391683/114‑1 (MQ=255)
                                                    gCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGtt                                                                                             >  3:354798/1‑149 (MQ=255)
                                                     cGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTc                                                                                            <  4:762755/149‑1 (MQ=255)
                                                       gTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCaa                                                                                          >  3:693303/1‑149 (MQ=255)
                                                           aTGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAgcgc                                                                                      <  1:567553/149‑1 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAacac                                                                                 >  1:807614/1‑149 (MQ=255)
                                                                tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAacac                                                                                 >  4:453360/1‑149 (MQ=255)
                                                                    tCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGa                                                                             <  1:848606/149‑1 (MQ=255)
                                                                      cTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAAc                                                                           >  4:217034/1‑149 (MQ=255)
                                                                       tCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACt                                                                          >  2:495867/1‑149 (MQ=255)
                                                                         gATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTaa                                                                        <  4:181677/149‑1 (MQ=255)
                                                                           ttattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAg                                                                      <  2:320360/149‑1 (MQ=255)
                                                                                         tGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAcc                                                        >  4:35163/1‑149 (MQ=255)
                                                                                           tttAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGt                                                      >  4:84893/1‑149 (MQ=255)
                                                                                             tAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCgg                                                              <  2:772653/139‑1 (MQ=255)
                                                                                             tAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCgg                                                              >  1:772653/1‑139 (MQ=255)
                                                                                              aaaGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAAt                                                   >  2:250127/1‑149 (MQ=255)
                                                                                                  tCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTc                                               <  1:495867/149‑1 (MQ=255)
                                                                                                              agcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGcc                                   <  3:217034/149‑1 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                  >  2:563800/1‑149 (MQ=255)
                                                                                                               gcagGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCt                                  >  2:254778/1‑149 (MQ=255)
                                                                                                                 agGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGc                                <  3:539583/149‑1 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATaaa                                                                                                             <  1:80591/71‑1 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATaaa                                                                                                             >  2:80591/1‑71 (MQ=255)
                                                                                                                  ggTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCt                               <  4:354798/149‑1 (MQ=255)
                                                                                                                    tAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTaa                             >  1:263649/1‑149 (MQ=255)
                                                                                                                     aaaaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTaaa                            <  2:546693/149‑1 (MQ=255)
                                                                                                                      aaaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAAt                           >  4:847593/1‑149 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                          <  4:693303/149‑1 (MQ=255)
                                                                                                                       aaGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATa                          >  2:94901/1‑149 (MQ=255)
                                                                                                                          ccTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATagag                       <  1:250127/149‑1 (MQ=255)
                                                                                                                             gACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCgg                    >  2:231175/1‑149 (MQ=255)
                                                                                                                              aCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGa                   <  2:202787/149‑1 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt               >  1:327051/1‑149 (MQ=255)
                                                                                                                                  ttcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt               <  2:263649/149‑1 (MQ=255)
                                                                                                                                     ttttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCt            <  3:84893/149‑1 (MQ=255)
                                                                                                                                       ttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg          >  1:822453/1‑149 (MQ=255)
                                                                                                                                       ttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg          >  1:920901/1‑149 (MQ=255)
                                                                                                                                        acgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAc                                                         >  4:238042/2‑101 (MQ=255)
                                                                                                                                        acgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATAc                                                         <  3:238042/100‑1 (MQ=255)
                                                                                                                                          gcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa       >  3:623013/1‑149 (MQ=255)
                                                                                                                                          gcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa       >  1:450140/1‑149 (MQ=255)
                                                                                                                                            gACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATg     >  1:519339/1‑149 (MQ=255)
                                                                                                                                            gACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATg     >  4:758273/1‑149 (MQ=255)
                                                                                                                                             aCTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGa    <  1:518857/149‑1 (MQ=255)
                                                                                                                                             aCTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGa    <  1:254778/149‑1 (MQ=255)
                                                                                                                                               tGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc  >  3:303739/1‑149 (MQ=255)
                                                                                                                                               tGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATc  <  2:223456/149‑1 (MQ=255)
                                                                                                                                                 gCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc             >  1:781203/1‑136 (MQ=255)
                                                                                                                                                 gCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGc             <  2:781203/136‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATC  >  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: