Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,831 C→A 7.2% intergenic (+71/+14) zapB → / ← rraA FtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,8310CA 7.2% 386.9 / 13.5 124intergenic (+71/+14)zapB/rraAFtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein
Reads supporting (aligned to +/- strand):  ref base C (53/62);  new base A (6/3);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 125 total aligned reads displayed.

AACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118683‑4118977
                                                                                                                                                    |                                                                                                                                                  
aacaacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGc                                                                                                                                                    >  8:112225/1‑149 (MQ=255)
 acaacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCg                                                                                                                                                   <  7:204846/149‑1 (MQ=255)
   aacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                                 >  8:51225/1‑149 (MQ=255)
   aacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                                 >  5:262268/1‑149 (MQ=255)
   aacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                                 <  8:356806/149‑1 (MQ=255)
       aTCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCt                                                                                                                                             <  3:401167/149‑1 (MQ=255)
       aTCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCt                                                                                                                                             >  5:94280/1‑149 (MQ=255)
          tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          >  3:356521/1‑149 (MQ=255)
          tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTCCGGTGCCttta                                                                                                                                          >  5:372034/1‑148 (MQ=255)
           gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                         >  6:32833/1‑149 (MQ=255)
            aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          <  5:223339/147‑1 (MQ=255)
            aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          <  8:466966/147‑1 (MQ=255)
            aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          <  7:112225/147‑1 (MQ=255)
            aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                        <  1:130450/149‑1 (MQ=255)
            aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCtt                                                                                                                                            <  7:392696/145‑1 (MQ=255)
             aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTa                                                                                                                                       <  4:371993/149‑1 (MQ=255)
              aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttct                                                                                                                                      <  2:162212/149‑3 (MQ=255)
                aaCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCa                                                                                                                                    <  1:342350/149‑1 (MQ=255)
                 aCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAt                                                                                                                                   >  3:402132/1‑149 (MQ=255)
                 aCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAt                                                                                                                                   >  3:144311/1‑149 (MQ=255)
                   agcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  >  5:328708/1‑148 (MQ=255)
                    tcagAACGGCTGGCAGGAACGTCTGGAGGCCGTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGGCTACAGGTAACGCAGAAAAAAGTCACCTTGCTGTGCCTTTCTTATCATTCa                                                                                                                                <  8:448395/148‑1 (MQ=255)
                     cagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                               >  7:454221/1‑149 (MQ=255)
                      agAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAat                                                                                                                              <  5:219380/149‑1 (MQ=255)
                       gAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAata                                                                                                                             >  3:215331/1‑149 (MQ=255)
                        aaCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                               >  5:108105/1‑146 (MQ=255)
                         aCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATc                                                                                                                           <  6:176141/149‑1 (MQ=255)
                           ggCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATc                                                                                                                           >  3:429036/1‑147 (MQ=255)
                               ggCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                        <  3:100291/146‑1 (MQ=255)
                               ggCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                        >  4:100291/1‑146 (MQ=255)
                                  aGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTCCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtc                                                                                                                  <  1:186898/149‑1 (MQ=255)
                                  aGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtc                                                                                                                  >  3:9614/1‑149 (MQ=255)
                                           tGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGa                                                                                                                    <  7:304642/138‑1 (MQ=255)
                                            gCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATACAATATCCAGCGGATCTTCTGATAg                                                                                                         >  2:445598/1‑148 (MQ=255)
                                            gCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGAGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCTATATCCAGCGGATCTTCTGAAAg                                                                                                         >  2:453925/1‑148 (MQ=255)
                                             cAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  4:402132/146‑1 (MQ=255)
                                             cAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                         <  2:389287/147‑1 (MQ=255)
                                                gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  >  7:355309/1‑119 (MQ=255)
                                                gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                    >  5:351758/1‑149 (MQ=255)
                                                gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                    <  7:235421/149‑1 (MQ=255)
                                                gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                    >  8:365495/1‑149 (MQ=255)
                                                  ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATc                                                                                                  >  4:421641/1‑149 (MQ=255)
                                                   ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  >  8:306948/1‑116 (MQ=255)
                                                   ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  <  7:306948/116‑1 (MQ=255)
                                                      ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                              >  8:329639/1‑149 (MQ=255)
                                                      ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCg                                                                                               <  6:328708/148‑1 (MQ=255)
                                                      ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCg                                                                                               <  6:262268/148‑1 (MQ=255)
                                                      ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCg                                                                                               <  5:254548/148‑1 (MQ=255)
                                                       tgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                              >  5:397370/1‑148 (MQ=255)
                                                         ggTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGt                                                                                             <  4:199743/147‑1 (MQ=255)
                                                                tGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                    <  4:215331/149‑1 (MQ=255)
                                                                  gaatcggTCTGATTACGCGTCATCTTCCGCGTAATTTCTTTACTTTGTAGGCCAACAGGTAACGCATAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                   <  2:319677/143‑1 (MQ=255)
                                                                  gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCAt                                                                                  <  6:351758/149‑1 (MQ=255)
                                                                        gTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGAt                                                                            <  7:329639/149‑1 (MQ=255)
                                                                             aTTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                    >  1:49131/1‑120 (MQ=255)
                                                                             aTTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTc                                                                        <  8:454221/148‑1 (MQ=255)
                                                                             aTTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTc                                                                        <  8:349156/148‑1 (MQ=255)
                                                                                    ttCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGgag                                                                 >  5:595/1‑148 (MQ=255)
                                                                                         cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           >  6:230043/1‑149 (MQ=255)
                                                                                         cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           >  2:409348/1‑149 (MQ=255)
                                                                                         cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           >  2:84374/1‑149 (MQ=255)
                                                                                            ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           <  7:365495/146‑1 (MQ=255)
                                                                                            ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           <  6:94280/146‑1 (MQ=255)
                                                                                            ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  4:429036/148‑1 (MQ=255)
                                                                                             cgcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                       <  4:196029/149‑1 (MQ=255)
                                                                                             cgcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                       >  1:175990/1‑149 (MQ=255)
                                                                                              gcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTa                                                                                            >  2:300144/1‑111 (MQ=255)
                                                                                              gcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTa                                                                                            <  1:300144/111‑1 (MQ=255)
                                                                                              gcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                       >  8:61253/1‑148 (MQ=255)
                                                                                                gTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      <  6:429286/115‑1 (MQ=255)
                                                                                                gTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      >  5:429286/1‑115 (MQ=255)
                                                                                                     ttCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGcc                                                >  6:420438/1‑148 (MQ=255)
                                                                                                     ttCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGcc                                                >  2:98067/1‑148 (MQ=255)
                                                                                                      tCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAACGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCa                                               <  4:143946/148‑1 (MQ=255)
                                                                                                            cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                         <  6:397370/148‑1 (MQ=255)
                                                                                                                 tAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                    >  2:466668/1‑148 (MQ=255)
                                                                                                                    gCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                >  8:352478/1‑149 (MQ=255)
                                                                                                                      cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCACGGTATTGTCAGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                              >  6:321686/1‑149 (MQ=255)
                                                                                                                         ctggTAACGAAGAAAAAAGGCACCTTGCGGTGGCTTTCTTATCATTCAATATCCTGCGGATCTTCTGAAATAATAATCCCGGTATTGTCGGCTTAAAGATGGTCGCCGGAGAAGAAGGTGAGACCGCCAAAATTGa                                        <  6:210231/134‑1 (MQ=255)
                                                                                                                         cAGGTAACGCAGAAAAAAGGCACCTTGAGGTGCCTTTCTTATCATTCAATATACATCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCTGAGAAGAAGGTGACACCGACAAAATTGa                                        >  5:210231/1‑136 (MQ=255)
                                                                                                                            gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      <  7:86669/87‑1 (MQ=255)
                                                                                                                            gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                        <  1:409348/149‑1 (MQ=255)
                                                                                                                            gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                        >  3:27467/1‑149 (MQ=255)
                                                                                                                               aCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGc                      <  6:416091/148‑1 (MQ=255)
                                                                                                                               aCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGc                      >  8:231755/1‑148 (MQ=255)
                                                                                                                                 gCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  4:296352/1‑149 (MQ=255)
                                                                                                                                  gataagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAt                  >  8:207786/7‑149 (MQ=255)
                                                                                                                                   aGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg                 <  3:296352/149‑1 (MQ=255)
                                                                                                                                    taagaaaGGCACCGCAAGGTGCCTTTTTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGc                >  8:59534/5‑149 (MQ=255)
                                                                                                                                      agaaaGGCACCGCAAGGTGCCTTTTTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCTGCATAAAGATGGTCGCCGGAGAAGAAGGTGCCACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCt              >  6:32028/3‑149 (MQ=255)
                                                                                                                                       gaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt             >  1:11487/2‑149 (MQ=255)
                                                                                                                                       aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  3:406542/105‑1 (MQ=255)
                                                                                                                                       aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         >  4:406542/1‑105 (MQ=255)
                                                                                                                                        aaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTc            <  5:420438/149‑1 (MQ=255)
                                                                                                                                        aaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTc            <  3:420202/149‑1 (MQ=255)
                                                                                                                                               ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCg     <  1:396017/149‑1 (MQ=255)
                                                                                                                                               ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCg     <  4:393703/149‑1 (MQ=255)
                                                                                                                                                 ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   <  1:363633/133‑1 (MQ=255)
                                                                                                                                                 ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  2:363633/1‑133 (MQ=255)
                                                                                                                                                   gCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCa  <  1:466668/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118683‑4118977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: