Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,6360AC42.0% 131.6 / 16.4 81L52V (TTA→GTA) fliARNA polymerase, sigma 28 (sigma F) factor
Reads supporting (aligned to +/- strand):  ref base A (36/11);  new base C (0/34);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 2.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGA  >  NC_000913/2001488‑2001783
                                                                                                                                                    |                                                                                                                                                   
gCCACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAAcc                                                                                                                                                   <  1:187045/151‑1 (MQ=255)
gCCACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAAc                                                                                                                                                    <  1:714927/150‑1 (MQ=255)
 ccACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAAccc                                                                                                                                                  >  1:682375/1‑151 (MQ=255)
  cACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAAccc                                                                                                                                                  >  2:703526/1‑150 (MQ=255)
      tCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAAt                                                                                                                                               <  1:1335757/149‑1 (MQ=255)
      tCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATgc                                                                                                                                             >  1:304046/1‑151 (MQ=255)
       cACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATgcc                                                                                                                                            <  1:684551/151‑1 (MQ=255)
          gcgcgTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATgccgcc                                                                                                                                         >  1:592545/1‑151 (MQ=255)
             cgTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGccgcccgc                                                                                                                                      >  2:1138675/1‑151 (MQ=255)
             cgTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGccgcccgc                                                                                                                                      >  2:1130778/1‑151 (MQ=255)
                tGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTg                                                                                                                                   >  2:1074387/1‑151 (MQ=255)
                  cGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGt                                                                                                                                  >  2:350053/1‑150 (MQ=255)
                  cGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGt                                                                                                                                  >  2:1111273/1‑150 (MQ=255)
                   gTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGt                                                                                                                                  >  2:201387/1‑149 (MQ=255)
                         cGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGc                                                                                                                               >  2:1113132/1‑146 (MQ=255)
                            tGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTc                                                                                                                       >  1:1261133/1‑151 (MQ=255)
                            tGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTc                                                                                                                       >  1:1055110/1‑151 (MQ=255)
                                      ccAGTCACGGCTGCGAAGTTCATCCAGCATCGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTCCCTTGTGGGGCGTCATAGCGTTCGACGGCATTAAGTACCCCATTGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACg                                                                                                             <  2:516943/151‑1 (MQ=255)
                                       cAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGc                                                                                                            >  1:1300123/1‑151 (MQ=255)
                                        aGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCt                                                                                                           >  1:801840/1‑151 (MQ=255)
                                             cggcggCAAATTTCACCCAGAATAGCGCCACGGATACGCTCCCCTGCGTAAGTTGAAAAGCCCGTCCTTTGTGGGCCGTCATGGCTTTCGCCGGCTTAAAGTACCCAAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGc                                                                                                            <  1:550880/145‑1 (MQ=255)
                                             cGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCgg                                                                                                      >  2:1002674/1‑151 (MQ=255)
                                             cGGCTGCGAAGTCCATCCAGCATAGCCCCACGGATCCGCTGCACTGCGTAAGTTGAAAATGCCGTTCCTTGTGGGTCGTCATAGCGTCCGACGGCATTAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTcgcg                                                                                                       <  2:936546/150‑1 (MQ=255)
                                              ggCTGCGAATTTCATCCAGCATAGCGCCCCGGATACGCTGCCCTGGGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTACCCAAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCggg                                                                                                     <  2:819744/151‑1 (MQ=255)
                                              ggCTGCGAAGTTCATCCGGAATAGCGCCACGGATACGCTGCACTGCGTAAGTTGAAAATGCGGTCCCTTGTGGGGCGTCATAGCGTGCGGCGGCATAAAGTCCCCCAGTGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCgg                                                                                                      <  2:470166/150‑1 (MQ=255)
                                               gCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGc                                                                                                    >  1:1013289/1‑151 (MQ=255)
                                                  gCGAAGTTCATCCAGCATAGCGCCACGGATACGCTCCACTGCGTAAGTTGAAAATGCCTTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGAACCCAAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGt                                                                                                 <  2:1049828/151‑1 (MQ=255)
                                                    gAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCg                                                                                               >  2:963144/1‑151 (MQ=255)
                                                         tCATCCGGCATAGCCCCAGGGATCCGCTGCACGGCTTAGGTTGTAAATCCCTTTCCTTGTAGGGCGCCATAGCTTTCGCCGGCATTAAGTACCCCAGGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAAcc                                                                                           <  1:485933/150‑1 (MQ=255)
                                                           aTCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGAAAATGCCGTTCCTTGTAGGGCGTCATAGCTTCCGCCGGCATTAAGTACCCAAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCt                                                                                          <  2:162856/149‑1 (MQ=255)
                                                             ccAGCATGGCGCCACGGACCCGCTGCCCTGCGTAATTTGAAAAGCCCGTCCCTTGTAGGCCACCTTGGCGTGCCCCGGCATTAAGTAACCCACGCCCCCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCa                                                                                       <  1:325254/150‑1 (MQ=255)
                                                             ccAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGc                                                                                        >  2:789043/1‑149 (MQ=255)
                                                             ccAGCAAAGCGCCCGGGATCCCCTGCACTGCGTAATTTTAAAATCCCTTTCTTTGTGGGGCGTCATAGCGTCCGCCGGCATAAAGTACCCCAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCa                                                                                       <  1:226459/150‑1 (MQ=255)
                                                               aGCATAGCGCCACGGATCCGCTGCACTGCGTAATTTTTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAgg                                                                                     <  2:59757/150‑1 (MQ=255)
                                                                    aGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAAt                                                                               >  2:267975/1‑151 (MQ=255)
                                                                    aGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAAt                                                                               >  2:720226/1‑151 (MQ=255)
                                                                      ccccccGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGCCGGCATAAAGTAACCAAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGc                                                                             <  2:499486/146‑1 (MQ=255)
                                                                       gcCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCt                                                                            >  2:275450/1‑151 (MQ=255)
                                                                        ccACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATAAATTACCCAAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCtt                                                                           <  2:853928/151‑1 (MQ=255)
                                                                                 cGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGc                                                                   >  1:935697/1‑150 (MQ=255)
                                                                                  gCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCa                                                                  >  2:797140/1‑150 (MQ=255)
                                                                                  gCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCAAACCTGCAGGCGCAATGCTTCGTGACGCAc                                                                 >  2:748623/1‑151 (MQ=255)
                                                                                    tGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCAcc                                                                >  2:746867/1‑150 (MQ=255)
                                                                                      cccTGCGTAAGTTGAAAACGCCGTTCCTTGGGGGGCGTCATAGCGTCCGCCGGCATTAAGAACCCCAATCCCCCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGAc                                                                     <  2:1121067/141‑1 (MQ=255)
                                                                                          gCGTAATTTGAAAATGCCGTCCCTTGTAGGGCGTCATAGCGTCCGACGGCATAAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCg                                                            <  2:253823/148‑1 (MQ=255)
                                                                                          gCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCggg                                                          >  1:1006068/1‑150 (MQ=255)
                                                                                          gCGTAAGTTGAAAAGGCCGTTCCTGGTGGGGCGTCATAGCGTTCGACGGCATTAAGAACCCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGa                                                         <  1:253435/151‑1 (MQ=255)
                                                                                              aaGTTGTAAATGCCTTCCTTTGGAGGGCTTCTTGGGGTCCGACGGAATTAAGAACCCCATGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACa                                                       <  1:1130719/149‑1 (MQ=255)
                                                                                              aaGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATAAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGa                                                         <  2:1174655/147‑1 (MQ=255)
                                                                                              aaGTTGAAAATGCCGTTCTTTGTGGGGCGTCATAGCGTTCGACGGCATAAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACa                                                       <  2:363906/149‑1 (MQ=255)
                                                                                              aaGTTGAAAATGCCGTTCCTTGTGGGGCGTCATAGCGTCCGACGGCATTAAGTACCCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACa                                                       <  1:924184/149‑1 (MQ=255)
                                                                                                gTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAAc                                                   >  1:200069/1‑151 (MQ=255)
                                                                                                 ttGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAAc                                                   <  2:87629/150‑1 (MQ=255)
                                                                                                        tGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                           >  1:384589/1‑151 (MQ=255)
                                                                                                          ccttcctttGTGGGGCATCATAGCTTTCGACGGCATAAAGTACCCAAACCCCCCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcaca                                          <  2:385003/143‑1 (MQ=255)
                                                                                                          ccGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                           >  1:497586/1‑149 (MQ=255)
                                                                                                          ccGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                           >  2:1167597/1‑149 (MQ=255)
                                                                                                             ttCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTACCCCAAGCCCCCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCAACAGCg                                       <  2:1233674/150‑1 (MQ=255)
                                                                                                             ttCCTGGTAGGGCGTCATGGCGTTCGCCGGCATTAAGTAACCAAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCg                                       <  1:1313066/150‑1 (MQ=255)
                                                                                                             tcccTTGTAGGGCGTCATGGCGTTCGACGGCATTAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGc                                        <  2:1233929/147‑1 (MQ=255)
                                                                                                               ccTTGTAGGGCTTCATAGGGTCCGACGGCATTAAATACCCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgt                                    <  1:1002621/151‑1 (MQ=255)
                                                                                                                 ttGTGGGGCGTCATGGCGTTCGACGGCATTAATAACCCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAg                                     <  1:393535/148‑1 (MQ=255)
                                                                                                                 ttGTAGGGCGTCATCGCGTTCGACGGCATAAAGTACCCCAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgtg                                   <  2:766619/150‑1 (MQ=255)
                                                                                                                 ttGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgtgt                                  >  2:259112/1‑151 (MQ=255)
                                                                                                                 ttGGAGGGCGTCATGGCGTCCGCCGGCATTAATTACCCAAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgt                                    <  1:757824/149‑1 (MQ=255)
                                                                                                                  tGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGtt                                 >  1:1120583/1‑151 (MQ=255)
                                                                                                                   gtggggCGTCATAGCGTTCGACGCCATTAAGAACCCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGtt                                 <  1:1247760/147‑1 (MQ=255)
                                                                                                                   ggggggCGTCATGGCGTTCGACGGCATAAAGTAACCCAAGGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGttt                                <  2:389363/148‑1 (MQ=255)
                                                                                                                       ggCGTCATAGCGTTCGACGGAATTAAGTACCCAAAGCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATcc                            <  2:614252/151‑1 (MQ=255)
                                                                                                                             aTAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTac                       >  1:391780/1‑150 (MQ=255)
                                                                                                                             aTAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTac                       >  1:124450/1‑150 (MQ=255)
                                                                                                                                gCGTTCGCGGGCATAACGCACCCCACTCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCt                   <  1:411750/151‑1 (MQ=255)
                                                                                                                                   ttCGACGGCATAAACTACCCAAAGCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCtt                  <  1:1184083/149‑1 (MQ=255)
                                                                                                                                   tgcggcGGAATAAAGTACCCCAAGCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTc                 <  1:367983/145‑1 (MQ=255)
                                                                                                                                   tcggacgaaaTTAAGAAACCCAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCa                           <  2:312568/131‑1 (MQ=255)
                                                                                                                                      gACGGCATTAAGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCGCGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCg             <  2:730389/151‑1 (MQ=255)
                                                                                                                                       aCGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCgg            >  2:176450/1‑151 (MQ=255)
                                                                                                                                       aCGGCATAAAGTACCCAAATCCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCg             <  2:615583/150‑1 (MQ=255)
                                                                                                                                             ttAAGTACCCAAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATag       <  2:888035/150‑1 (MQ=255)
                                                                                                                                               aaaaaaCCCAATCCCCCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATaga      <  1:295488/145‑1 (MQ=255)
                                                                                                                                                aGTACCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTg   <  1:249631/151‑1 (MQ=255)
                                                                                                                                                 gTACCCCAAACCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGa  <  1:1060063/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCACTTCACGCGCGTTGCGTCGCACGCTGCGCGGCACCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGA  >  NC_000913/2001488‑2001783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: