Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,6580AG21.2% 172.5 / 13.3 86E538E (GAA→GAG) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (41/26);  new base G (0/18);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAT  >  NC_000913/3214513‑3214804
                                                                                                                                                 |                                                                                                                                                   
cGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAACGc                                                                                                                                                  >  1:632104/1‑149 (MQ=255)
cGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGACAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                                                                                                                                  >  1:893893/1‑149 (MQ=255)
    aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCAGCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGGAAGCCtgc                                                                                                                                              >  1:69298/1‑149 (MQ=255)
      aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgt                                                                                                                                            >  1:772551/1‑149 (MQ=255)
       cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtg                                                                                                                                           >  2:1239469/1‑149 (MQ=255)
          gATCCGCACAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCgg                                                                                                                                        <  1:1107942/149‑1 (MQ=255)
           aTCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCGCCCTCGGGCGGCCGCTGGATGCTGCGACCACCGAAAGCCTGCGTGCGGc                                                                                                                                       >  1:959999/1‑149 (MQ=255)
             ccGCAAAGTGCTGAAGATCGCCAAAAAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCaa                                                                                                                                     >  1:628276/1‑149 (MQ=255)
              cGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCTTGCGTGCGGCAAc                                                                                                                                    <  1:852428/149‑1 (MQ=255)
              cGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAAc                                                                                                                                    >  1:795934/1‑149 (MQ=255)
                     tGCTGAAGAACGCCAAAGCGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATGTCAACTAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgac                                                                                                                             <  1:174764/149‑1 (MQ=255)
                      gCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgacg                                                                                                                            >  1:766763/1‑149 (MQ=255)
                        tGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGGTTCTGCGACCACCGGAAGCCTGCGGGCGGCCCGGCGCGACGTg                                                                                                                          >  2:801478/1‑149 (MQ=255)
                          aaGAGCGCCAACGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGTTTCGATTCTGGGGGATTTCATCGAGGATACCACCCTTTTTCTTTCTAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgct                                                                                                                        <  2:530992/149‑1 (MQ=255)
                           aGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctg                                                                                                                       <  1:622562/149‑1 (MQ=255)
                            gATCGCCAAAGAGCCAATCGCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctgg                                                                                                                      <  1:842589/149‑1 (MQ=255)
                               cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGGGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCCCGACCTgcgggctg                                                                                                                   >  2:53311/1‑149 (MQ=255)
                                tccAAAGAGCCAATCTCGATGGAAGCGCCGATGGGTGGTGATGAAGATTCGCATCTGGGGGATTTCAGCGAGGATACCACCCTTTTTTTTTTTCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGctggctgg                                                                                                                  <  2:323767/148‑1 (MQ=255)
                                 ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                                                                   >  1:199548/1‑147 (MQ=255)
                                 ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                                                                   <  2:199548/147‑1 (MQ=255)
                                 ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                                 <  1:1145135/149‑1 (MQ=255)
                                 ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCAGCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctgac                                                                                                                 >  1:316010/1‑147 (MQ=255)
                                   aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGCCGTGCTGGCTGGCCt                                                                                                               >  1:763464/1‑149 (MQ=255)
                                       cgCCAATCTCCATGGAAACGGCGATCGGTGATGATGAAGACTCGCATCTGGGGGATTTCATCGAGGATACCACCCTTGTGTTTTTTAATGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAcc                                                                                                           <  2:1084660/148‑1 (MQ=255)
                                       agCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCCCGACGTCCTGGCTGGCCGGAcc                                                                                                           >  2:643131/1‑149 (MQ=255)
                                        gCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcg                                                                                                          >  1:529068/1‑149 (MQ=255)
                                         ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGCGCTGCCGCTGGGTTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgc                                                                                                         >  1:237574/1‑149 (MQ=255)
                                          cAATCTCCATGGAAACGCCGATCGGTGATGATGAGGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCCGGCCTGAGcggg                                                                                                        >  2:1161272/1‑147 (MQ=255)
                                           aacctcCAGGGAAGCGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTTCTGTCTATTGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgc                                                                                                       <  2:54392/146‑1 (MQ=255)
                                                  aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGc                                                                                                >  1:503670/1‑149 (MQ=255)
                                                           tcggtCGGTGAGTATGAAGATTCGCATCTGGGGGTGTTCATCGGGGATACCACCCTTGTGCTTTTGCAGGATACTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                       <  2:934831/145‑1 (MQ=255)
                                                           ccGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATTGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                       <  1:630028/149‑1 (MQ=255)
                                                             gATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTCCGTGCGGCAACGCGCGACGCGCGGGCTGGCCTGACCCAGCGTGAAGCATAAGATCAGc                                                                                     >  2:1020719/1‑149 (MQ=255)
                                                              ggcGGTTATGATGAAGTTTCGCATCTGGGGGATGTCGTCGAGGATACCACCCTCGTGTTGTCTAAGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                                                                    <  2:1177399/147‑1 (MQ=255)
                                                                cGGTGATGTGGAAGCTTCGCCTCTGGGGGATTTCATCGAGGATACCACCCTCGTTTTTCCGAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTa                                                                                  <  2:1087205/149‑1 (MQ=255)
                                                                cGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGAGCGGCAACGCACGACGTGCTGGCTCGCCTGACCGCGCGCGAAGACAAAGTTCTGCgga                                                                                  >  2:681423/1‑147 (MQ=255)
                                                                   tgatgatgaCGACTCGCATCTGGGTGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTGCTGCGTATGc                                                                               <  1:1015474/149‑1 (MQ=255)
                                                                   tgatgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGACAGGCCTGACCGCGCGTGAGGCAAAAGTTCAGCGCATGc                                                                               >  2:644911/1‑149 (MQ=255)
                                                                     atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGt                                                                             >  1:1187179/1‑149 (MQ=255)
                                                                      tgatgaAGATTCGTATCTGTGGGATTTCATCGCGGATCCCACCCTTGTTTTTTTTCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGtt                                                                            <  2:99343/149‑1 (MQ=255)
                                                                        atgaAGGTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                                                                          <  2:842581/149‑1 (MQ=255)
                                                                        atgaAGATTCGTATCTGGGGGATTTCATCGAGGTTACCACCTTTTTGCTTTTGATGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                                                                          <  2:230500/149‑1 (MQ=255)
                                                                        atgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                                                                          >  2:206796/1‑149 (MQ=255)
                                                                         tgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGGCCGCGCGTGAAGCAAAAGTTCGTCGGATGCGGTTCg                                                                         >  2:318392/1‑149 (MQ=255)
                                                                             gtgtCGCATTGGGGGGCTTTCATCGAGGATACCCCCCTCGTGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTAt                                                                     <  2:796929/146‑1 (MQ=255)
                                                                              aTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATc                                                                    <  2:503670/149‑1 (MQ=255)
                                                                                 cGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGat                                                                 <  1:1187607/149‑1 (MQ=255)
                                                                                 cGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGat                                                                 <  1:1189308/149‑1 (MQ=255)
                                                                                      cTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGaa                                                            <  2:956638/149‑1 (MQ=255)
                                                                                         ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac                                                         <  1:1239469/149‑1 (MQ=255)
                                                                                         ggaggtttCATCGACCATAACACCCCCGACCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac                                                         <  1:594107/144‑1 (MQ=255)
                                                                                            gATTTCATCGAGGATACCACCCTCGAGGTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGa                                                      <  2:217656/149‑1 (MQ=255)
                                                                                              tttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACCCCGACt                                                    >  1:722576/1‑149 (MQ=255)
                                                                                              tttCATCGAGGATACCACCCTCGAGCTGCCGCTAGATTCTGCGACCACCGAAAACCTGCGTTCGGCAACTCACGACCTGCTAGCTGCCCGGACCACGCGAGAAACAAAAGGTCTGCGGATGCGTTGCGGTATGGATCTGCACACCgcag                                                    >  2:412871/1‑146 (MQ=255)
                                                                                                tCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTac                                                  <  2:795934/149‑1 (MQ=255)
                                                                                                 cATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTaca                                                 <  1:943686/149‑1 (MQ=255)
                                                                                                   tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACtcaccg                                               >  1:587457/1‑144 (MQ=255)
                                                                                                   tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGCATGCGATGCGCTATCGGTATAGACACCGACGAGACg                                               >  2:933042/1‑149 (MQ=255)
                                                                                                   tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAACGCGCGTGAAGCAGAAGATCAGCGTATGCGGTTCAGTAGCGATCGGAACAGCCACAACACg                                               >  2:1093506/1‑149 (MQ=255)
                                                                                                    cGAGGACCCCACCCTTGTGCTGCTGCAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGc                                              <  2:993105/149‑1 (MQ=255)
                                                                                                         aTACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGgaa                                         <  1:206796/149‑1 (MQ=255)
                                                                                                             caccCTCGCGCTGTCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCCGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGaaga‑ag                                    <  2:632104/149‑1 (MQ=255)
                                                                                                               cccTCGATCTGCCGCTTGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTg                                  <  2:929274/149‑1 (MQ=255)
                                                                                                                ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTgg                                 >  2:1022589/1‑149 (MQ=255)
                                                                                                                   cgtgtttccGATGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTa                              <  2:158254/142‑1 (MQ=255)
                                                                                                                    gtgCTGCCTCATGAGTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTaa                             <  2:772551/147‑1 (MQ=255)
                                                                                                                     aGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGCCGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTaaa                            <  1:392519/149‑1 (MQ=255)
                                                                                                                          ccGCGGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGtt                       <  2:189971/149‑1 (MQ=255)
                                                                                                                           cgcgGGGTTCGGCAACCACCGAAAGCGTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTc                      <  2:286892/149‑1 (MQ=255)
                                                                                                                           cgatGGATTCGGCGACCACCGAAAGCCTGAGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTc                      <  2:915767/146‑1 (MQ=255)
                                                                                                                             gtGGGTTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTGCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGa                    <  1:437175/148‑1 (MQ=255)
                                                                                                                             cTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGa                    >  1:655014/1‑149 (MQ=255)
                                                                                                                               tgatacggCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACg                  <  2:882468/142‑1 (MQ=255)
                                                                                                                                  ttCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCGGGCCCGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCGGGAACA‑GGTGCGTAAACCGTTCGACGTTa               >  2:206159/1‑149 (MQ=255)
                                                                                                                                  tccTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGGGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTa               <  1:333903/147‑1 (MQ=255)
                                                                                                                                     tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTAccc            >  2:445832/1‑149 (MQ=255)
                                                                                                                                       cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGatat                                                               >  1:852969/1‑97 (MQ=255)
                                                                                                                                       cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGatat                                                               <  2:852969/97‑1 (MQ=255)
                                                                                                                                        gacccacGAAAGCCTGCGTGCGGGGACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCACAAGATCTGCGTGTGCGCTTCCGTAACGATAGGAACACCGACTACACACTCGGAGA‑AGTGGGGAGACAGTTCGCCGCTACCcccg         >  2:1137989/7‑149 (MQ=255)
                                                                                                                                        gACCACCGAAAGCCTGCGTGCGGCAACGCACGACGCGCTAACCGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGTGATCGGTAGGAACACTGACTAAACACTGGAAAA‑AGTGGCTAAACAGTACGCCCTTGACcgcg         >  2:578464/1‑149 (MQ=255)
                                                                                                                                        gACCACCGAAAGCCTGCGTGCGGCAACGCAAAAAGGGCTGGGAGGCATGCCCGCGTGTGAAGAAAAAGTGCTGCGTATGCGCTGCCGTCACGATATGAACACCGACTAAACGCTGGAACATACGGGGCAAAC‑GTTCAACGTTACCcggg         >  2:579556/1‑147 (MQ=255)
                                                                                                                                          ccaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGaa       <  1:534611/149‑1 (MQ=255)
                                                                                                                                           cgccGACAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAAc      <  1:283758/147‑1 (MQ=255)
                                                                                                                                           caccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGGGAGGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTGAACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAAc      >  2:629811/1‑149 (MQ=255)
                                                                                                                                              cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAACGTa   >  1:969150/1‑149 (MQ=255)
                                                                                                                                               gAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAACGttt  >  1:1064359/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGA‑AGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAT  >  NC_000913/3214513‑3214804

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