Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,214,658 A→G 16.3% E538E (GAA→GAG)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,6580AG16.3% 203.9 / 9.6 103E538E (GAA→GAG) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (50/33);  new base G (2/18);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.23e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGT  >  NC_000913/3214512‑3214802
                                                                                                                                                   |                                                                                                                                                
ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGCCCACCGAAAg                                                                                                                                                >  2:100253/1‑149 (MQ=255)
ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGAGTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAg                                                                                                                                                >  1:1022250/1‑149 (MQ=255)
   gAAGACAAGATCCGCAAAGTGCGGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCt                                                                                                                                             <  1:204049/149‑1 (MQ=255)
              ccGCAAAGTGCGGAAGGTCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTTCCGCTGGATTCTGCCACCACCGAAAGCCCTCGGGCGGCac                                                                                                                                  >  2:635763/1‑148 (MQ=255)
               cGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGAACACCGAAAGCCTGCGTGCGGCAAc                                                                                                                                 >  1:636752/1‑149 (MQ=255)
                gCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                                                                                                                                <  1:1009319/149‑1 (MQ=255)
                 cAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCTTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGc                                                                                                                               <  1:1287842/149‑1 (MQ=255)
                 cAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGc                                                                                                                               >  1:822807/1‑149 (MQ=255)
                  aaaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCa                                                                                                                              <  1:486783/149‑1 (MQ=255)
                   aaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATAAAGCA‑CTTCACCTGGGGGATCAGAGCAATGATACCACCCTCGAGCTGGCGCTGGATGCTGCGACCGCCGAAAGCCTGCGTGCGGCAACGCac                                                                                                                             >  1:629193/1‑149 (MQ=255)
                     gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCCGCGCGCCGCCACGCacga                                                                                                                           >  2:547464/1‑149 (MQ=255)
                     gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTCCGACCAGCGAACGCCTGCGTGCGGCAACGCccga                                                                                                                           >  1:1287657/1‑149 (MQ=255)
                        cTGAGGATCGCCAAGGAGCCAATCTCCATGGAAACGCCGGTCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCAGCGAGGATACCACCCTCTTGTTGCCTAAGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGt                                                                                                                        <  2:1143921/149‑1 (MQ=255)
                         tGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTg                                                                                                                       >  1:299754/1‑149 (MQ=255)
                             gATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCAAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAACCCTGCGTGCGGCAACGCACGACGTgctgg                                                                                                                   >  1:1171191/1‑149 (MQ=255)
                               tCGCCAAAGGGCCAACCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggct                                                                                                                 <  1:161824/149‑1 (MQ=255)
                                cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                                                                <  1:1194338/149‑1 (MQ=255)
                                  ccAACGAGCCTATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                              <  1:672673/149‑1 (MQ=255)
                                  ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGAGCACCGAAAGCCTGCGTGGGGGAACGCGCGACGTGCtggctggc                                                                                                              >  2:1047162/1‑149 (MQ=255)
                                   cAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc                                                                                                             <  1:1038546/149‑1 (MQ=255)
                                    aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGTTGCCGAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCt                                                                                                            <  2:464530/149‑1 (MQ=255)
                                     aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAACGCCTGCGTGCGGCACGGCACGCGGTCCTCGCTGGCCTg                                                                                                           >  2:824057/1‑149 (MQ=255)
                                            aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgc                                                                                                    <  1:625040/149‑1 (MQ=255)
                                            aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgc                                                                                                    >  1:502056/1‑149 (MQ=255)
                                            aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgc                                                                                                    >  1:1233853/1‑149 (MQ=255)
                                            aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑CCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgc                                                                                                    <  1:42312/149‑1 (MQ=255)
                                             aTCTCCATCGAAACGCCCATCGGCGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTTTTGTTGCCGCTTGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                                                                   <  2:184536/149‑1 (MQ=255)
                                               ctcCAGGGAAACGCCGACCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTg                                                                                                 <  1:500383/149‑1 (MQ=255)
                                               ctcCACGGAACCGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGTCCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTg                                                                                                 <  1:1219126/149‑1 (MQ=255)
                                                  cagggAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                                                              <  1:880001/146‑1 (MQ=255)
                                                    tGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGCGCCGGCCGGCCGGACCGCGCGCGAAGCa                                                                                            >  2:477050/1‑149 (MQ=255)
                                                      gAAACGCCGATCGGTGATGATGTAGGT‑TCGCATCTGGGGGATTTCATCGGGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCGGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCaaa                                                                                          <  1:116048/149‑1 (MQ=255)
                                                      gAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCCTCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCaaa                                                                                          <  1:1210183/149‑1 (MQ=255)
                                                          cgctcaTCGGTGATGATGAAGAT‑TCGCATCTGGGTGATGTCATCGAGGCTACCACCCTCGTGCTGTTGAAGGATCCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGtt                                                                                      <  2:954756/144‑1 (MQ=255)
                                                          cGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGCTACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGACTGACCGCGCGTGAAGCAAAAGtt                                                                                      <  1:50962/149‑1 (MQ=255)
                                                          cGCCGATCGGTGATGATGAAGAT‑CCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGtt                                                                                      <  1:703717/149‑1 (MQ=255)
                                                           tccGATCCGTTATGAGGAAGAT‑TCCTCTCTGGGGGATTTCATCGAGGATACCGCCCTTTTTTTTTTTAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTc                                                                                     <  2:518509/148‑1 (MQ=255)
                                                           gcggATCGGTAATGTGGACGTT‑TGGCTTCTGGGGGATTTCATCGAGGATACCGCCCTTTTTTTTTTGATGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTc                                                                                     <  2:196323/146‑1 (MQ=255)
                                                            ccGATCGGTGATGATTAAGCTGTTGC‑TCTTGGGGATTTCAGCGAGGATACCACCCTCTGGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                    <  2:268727/149‑1 (MQ=255)
                                                            ccGATCGGTGATGATGAAGTT‑TGACATCTTGGTGATTTCATCGAGGATACCACCCTCGATCTGCCTCTGGATACGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                    <  2:268891/149‑1 (MQ=255)
                                                            ccGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCCCCGAAAGCCTGCGTGCGGCAACGCACGACGTGCCGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                    <  1:949106/149‑1 (MQ=255)
                                                            ccGATCGGTGATGATGAAGAT‑TCGCATCAACCGGTTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGGAGCAAAAGTTCt                                                                                    >  1:956782/1‑149 (MQ=255)
                                                              taccGGTGCTGATGATCCT‑TCTCATCTGGGCGATTCCACCGCGGACCCCACCCTCGATTTTCTTAATGATCCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGc                                                                                  <  2:622771/146‑1 (MQ=255)
                                                               agcGGTGATGATGAAGAT‑TCCCATCTGGGGGATTTCACCGAGGATACCACCCTCGTGTTTTTGAAGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                                                                 <  2:1140543/147‑1 (MQ=255)
                                                                    tgatgatgaAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGc                                                                            >  2:570172/1‑149 (MQ=255)
                                                                          tgaAGCG‑TCGCATCTCGGGGATTTCACCGAGGATCCCCCCCCTGAGTTGCTGAATGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCg                                                                      <  2:112472/149‑1 (MQ=255)
                                                                          tgaAGAT‑TCGCATCTGGGGGATTTCATCGAGGGTACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAGAGTTCTGCGTATCCGTTTCg                                                                      >  1:545338/1‑149 (MQ=255)
                                                                          tgaAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCCCGCCGCGCTGGCTCGCCCGACCGCGCGTGAACCAAAAGAACTGCGTACGCGtcgcg                                                                      >  2:589786/1‑145 (MQ=255)
                                                                          tgaAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATGCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCCGACCGCCCGAGCAGCAAGAGGTCGGCGAAGGCGTTCCg                                                                      >  2:1235174/1‑149 (MQ=255)
                                                                           gaAGAT‑TCGCATCTGGGGGATTTGATCGAGGATACCGCCCTCGAGCTGCCGCTGGATTCTGCGAACCCCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCgg                                                                     >  1:684941/1‑149 (MQ=255)
                                                                           gaAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCgg                                                                     >  1:205068/1‑149 (MQ=255)
                                                                            aaGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCggg                                                                    >  1:1142819/1‑148 (MQ=255)
                                                                            aaGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGt                                                                    >  1:1087553/1‑149 (MQ=255)
                                                                             aGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCAGCGTGCGGCAACGCACCCCGTACGGACTGGCCTGACCGCCCGGACAGCACAAGTTCCGCGTATGCGGTTCggaa                                                                   >  2:301680/1‑147 (MQ=255)
                                                                             aGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTCTGCGTTTCggga                                                                   >  1:49990/1‑147 (MQ=255)
                                                                                       tCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGGAGCAAAAGTTCTGCGTATGCGTTTGGGTATccgtgtga                                                          >  1:187919/1‑142 (MQ=255)
                                                                                        cTGGGTGTTTTCATCGAGGATCCCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGaa                                                         <  1:997155/149‑1 (MQ=255)
                                                                                         tGGGGGATTTCATCGAGGATACCACCCTCGAGATGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGGGCGGCAACGCACGACGTGCTGGCTGGCCGGACCGCGCGGGAAGCAAAAGTCCTCCCACTGCGTTCTGGCATCGATATGAac                                                        >  2:930659/1‑149 (MQ=255)
                                                                                          gggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGAGATGAaca                                                       >  1:353104/1‑149 (MQ=255)
                                                                                            gggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc                                                     >  2:1037260/1‑149 (MQ=255)
                                                                                            gggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGGGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc                                                     >  1:131112/1‑149 (MQ=255)
                                                                                               ttttCATCGAGGATACCGCCCTTGTGTTGTCTATTGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGAc                                                  <  2:1123334/148‑1 (MQ=255)
                                                                                               tattCATAGAGGGTACAAGCCTCCAGTTGCGCAAGGATTCGGCGACCACCGAAAGCCGGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGAc                                                  <  2:893957/147‑1 (MQ=255)
                                                                                               aTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCAGGACGTGCTGGCTGGCCTGACCGCGCGTGGAGCAAAAATACTGCGTGTGCGTTTCGGTGTCGATATGAACACAGAc                                                  >  1:271024/1‑149 (MQ=255)
                                                                                                tttCATCGAGGAGACAACCCTCGAGCTGCTGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACt                                                 <  2:305962/149‑1 (MQ=255)
                                                                                                 ttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGCCCACCGAAAGCCTGCGTGCGGCAACGCACGCCGGGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCGCCGTATGCGTTCCGGGATCGAAAGCAACACAGAcaa                                                >  2:580609/1‑147 (MQ=255)
                                                                                                 ttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGACCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATTAGTTTCGGTATCGATATGAACACCGACTa                                                <  1:845424/149‑1 (MQ=255)
                                                                                                 ttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCGGGCCCGACCGCGCGTGAAGCAAAAGTTCTGCGGATGCGTTTCGGGATCGAGATGAACACCGAGTa                                                >  2:729423/1‑149 (MQ=255)
                                                                                                    aTCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCCATGCACGACGTGCTGGCGTTCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGGCTacac                                             >  1:627885/1‑149 (MQ=255)
                                                                                                       gAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGACTGACAGAGCGCGAAGAAAAACTTATGCGTATGCGTTTCGGCATAGATATCAAAACCGACTCAACGca                                          >  2:1077909/1‑148 (MQ=255)
                                                                                                       gAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCGAGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTAGGCGTTTCGGGATCGATATGAACCACGACTAAACGca                                          >  2:671127/1‑148 (MQ=255)
                                                                                                        acgATACCACCCTCGTGCTGCCGCAGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTg                                         <  2:299723/147‑1 (MQ=255)
                                                                                                        aGGCTACCACCCTTCGGCTGGCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTg                                         <  2:554924/149‑1 (MQ=255)
                                                                                                           aTACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCGAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACCCGCTGgaa                                      >  1:1058259/1‑149 (MQ=255)
                                                                                                              ccccccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGgaagaa                                   <  1:160846/146‑1 (MQ=255)
                                                                                                              ccccccGCGCGCTGCCGCTGGATTCTGCGACCACCCAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGGTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGgaagaa                                   <  1:164970/146‑1 (MQ=255)
                                                                                                                tctcTCGAGTTGCCGAGGGAGTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                 <  2:1293630/146‑1 (MQ=255)
                                                                                                                ccccTCGATCTGCCTAGTGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                 <  2:624882/148‑1 (MQ=255)
                                                                                                                accCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                 <  1:430126/149‑1 (MQ=255)
                                                                                                                accCTCGAGCTGCCGCTGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                 <  2:476057/149‑1 (MQ=255)
                                                                                                                  ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTTGCCTGACCGCGCGTGAAGCACAAGTATAGAGATGGGGGTTGGGGATCGATATGAACACCGACTACACGCTGGAAGAAGtag                               >  1:37243/1‑147 (MQ=255)
                                                                                                                  ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGTGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                               <  1:1004172/149‑1 (MQ=255)
                                                                                                                  ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                               <  1:1128660/149‑1 (MQ=255)
                                                                                                                   cTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTggg                              >  1:226726/1‑149 (MQ=255)
                                                                                                                   cTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGACCTGCGCAAGCGGGTCGGCAGCGATACGAAAACCGACAACACGCTGgaagaaagcgg                              >  2:459397/1‑144 (MQ=255)
                                                                                                                               catgATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGa                  <  2:1297543/146‑1 (MQ=255)
                                                                                                                                ttgATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGAc                 <  2:266050/147‑1 (MQ=255)
                                                                                                                                tGGATTCTCCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGGGCGGGCTGGCCTTGCCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGAc                 <  2:363904/149‑1 (MQ=255)
                                                                                                                                  gATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACCCTGGAAGAAGTGGGTAACAAGTTCGACGt               >  1:23723/1‑149 (MQ=255)
                                                                                                                                    ttCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGGTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTCTGCGTTTCGGTGTGGATATGGAGACCGCGTACACGCTGGAAGACGGGGGTGAACAGTTTGACggtg             >  2:979798/1‑146 (MQ=255)
                                                                                                                                    tccggCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTa             <  2:720187/145‑1 (MQ=255)
                                                                                                                                     tCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGGTCGACGTTAc            >  1:1069068/1‑149 (MQ=255)
                                                                                                                                      cTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAcc           <  2:629193/149‑1 (MQ=255)
                                                                                                                                       tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc          >  1:982391/1‑149 (MQ=255)
                                                                                                                                       ggCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc          <  2:354893/148‑1 (MQ=255)
                                                                                                                                             caccGAAAGCCTGCGTGCGGCGACGCCCGACGTGCTGGCTGGGCTGACCGCGCGTGAAGCAAAAGGTCTGCGTATGCGTTTCGGTATCACTATGAACCCAAGCGAAACGCTGGAAGAAGTCGGTGAACCGATCGCCGGTACCcgcggca    >  2:289948/1‑146 (MQ=255)
                                                                                                                                             caccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGGAGAAGTGGGTAAACAGTTCGACGTGACCCGCGAAc    >  1:983049/1‑149 (MQ=255)
                                                                                                                                             caccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAAc    >  1:1298574/1‑149 (MQ=255)
                                                                                                                                               ccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGt  >  2:665103/1‑149 (MQ=255)
                                                                                                                                               ccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGCTACCCGCGAACGt  >  1:263709/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAT‑TCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGT  >  NC_000913/3214512‑3214802

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