Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,0670AG37.6% 94.0 / 88.1 85noncoding (309/2904 nt)rrlH23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (20/33);  new base G (8/24);  total (28/57)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTG  >  NC_000913/225937‑226154
                                                                                                                                  |                                                                                       
cGAACCGGGGGAACTGAAATATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGACTGGgaa                                                                                       >  2:300830/1‑133 (MQ=12)
cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGaaa                                                                                       >  2:493083/1‑133 (MQ=32)
cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGaaa                                                                                       >  1:375656/1‑133 (MQ=32)
cGAACCGGGGGAACTGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGgaa                                                                                       >  2:260627/1‑133 (MQ=17)
 gAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAg                                                                                      >  2:478163/1‑133 (MQ=17)
 gAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAg                                                                                      <  1:38507/133‑1 (MQ=17)
 gAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATGCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAg                                                                                      <  1:402496/133‑1 (MQ=12)
   aCCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGgc                                                                                    <  1:262931/133‑1 (MQ=17)
     cGGGGGAATTGAAACTTCCAAGTACCCCGACGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                                                                  <  1:130324/133‑1 (MQ=255)
     cGGGGGAACTGAAACTTCTAAGTACCCCGACGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                                                                  <  1:232807/133‑1 (MQ=17)
     cGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                                                                  <  1:455375/133‑1 (MQ=32)
      gggggAACTTCAACATCTAAGTACCCGGAGGAAAAGAAATCCACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGgcgcg                                                                                 <  2:241774/133‑1 (MQ=11)
      gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGATCGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTGGTGGAAGCGTCTGGGAAGgcgcg                                                                                 <  1:419354/133‑1 (MQ=12)
      gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                                                                 <  2:424169/133‑1 (MQ=32)
      gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                                                                 <  1:298374/133‑1 (MQ=32)
      gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                                                                 <  1:291009/133‑1 (MQ=32)
      gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGAAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                                                                 <  2:254261/133‑1 (MQ=17)
       ggggAACTGAGACATCTAAGTACCGCGAGGAAGAGAAATCAGCCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTG‑ATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGcgcgcg                                                                               <  2:221606/133‑1 (MQ=11)
       ggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGgcgcgc                                                                                >  1:107696/1‑133 (MQ=17)
       ggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTG‑ATCAGTGTGTCTGTTAGTGGAATCGTCTGGGACGGcgcgcg                                                                               >  1:221606/1‑133 (MQ=11)
         gtaaCTGAAAAATCTTAGTACCCCGAGGAAAAGAAATCAACTGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGa                                                                              <  2:244061/131‑1 (MQ=11)
         ggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGa                                                                              <  1:285922/133‑1 (MQ=17)
           aaCTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATa                                                                            <  1:461777/133‑1 (MQ=32)
            aCTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATAc                                                                           <  2:370150/133‑1 (MQ=32)
            aCTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATAc                                                                           <  1:377591/133‑1 (MQ=32)
              tGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAg                                                                         <  1:293555/133‑1 (MQ=32)
              tGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAg                                                                         <  1:424101/133‑1 (MQ=32)
                   cATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGa                                                                    <  1:464016/133‑1 (MQ=32)
                   cATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGa                                                                    <  2:474961/133‑1 (MQ=32)
                     tCTAAGTACCCCGAGGAAAAGAAATCAACCTAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACa                                                                  >  2:460375/1‑133 (MQ=12)
                        aaGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcc                                                               >  2:300091/1‑133 (MQ=32)
                        aaGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcc                                                               >  2:285697/1‑133 (MQ=32)
                          gTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGcccc                                                             <  2:77702/133‑1 (MQ=17)
                          gTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcccc                                                             >  1:80383/1‑133 (MQ=32)
                           tACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGGGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACCGGGTGACAGCCCCg                                                            <  2:255805/133‑1 (MQ=12)
                             ccccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTa                                                          <  1:91278/133‑1 (MQ=17)
                             ccccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTa                                                          >  1:117348/1‑133 (MQ=32)
                              cccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                         >  2:320861/1‑133 (MQ=17)
                              cccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                         <  2:461115/133‑1 (MQ=32)
                              cccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                         >  2:20755/1‑133 (MQ=32)
                              cccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                         >  2:103686/1‑133 (MQ=32)
                                cGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACCGGGTGACAGCCCCGTacac                                                       >  1:416560/1‑133 (MQ=17)
                                cGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTacac                                                       <  2:105682/133‑1 (MQ=17)
                                cGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTacac                                                       >  2:435220/1‑133 (MQ=32)
                                cGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTacac                                                       >  2:75573/1‑133 (MQ=32)
                                cGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTacac                                                       >  2:341438/1‑133 (MQ=32)
                                 gAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTAcaca                                                      >  1:33726/1‑133 (MQ=17)
                                  aGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaa                                                     >  1:428032/1‑133 (MQ=32)
                                   ggAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaa                                                    >  1:360910/1‑133 (MQ=17)
                                   ggAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGGGCGATACAGGGGGAGAGCCCCGTACACaaa                                                    >  1:198533/1‑133 (MQ=12)
                                   ggAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaa                                                    >  2:344418/1‑133 (MQ=32)
                                   gcaaaaGAAGTCATCTGAGATTCCACTACTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaa                                                    <  2:252178/131‑1 (MQ=11)
                                     aaaaGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaaaa                                                  <  1:266756/133‑1 (MQ=32)
                                       aaGAAATCAACCTAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATg                                                <  1:460375/133‑1 (MQ=17)
                                          aaaTCAACCGAGATTTCCCCAGTATCGGGGAGCGAACGGGGTGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGcac                                             <  2:182197/133‑1 (MQ=11)
                                          aaaTCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACA‑GGTGACAGCCCCGTACACAAAAATGcaca                                            <  2:289225/133‑1 (MQ=12)
                                           aaTCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGcaca                                            <  2:375656/133‑1 (MQ=32)
                                           aaTCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGcaca                                            <  2:487335/133‑1 (MQ=32)
                                            aTCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGGGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACAt                                           >  2:390015/1‑133 (MQ=17)
                                               aaCCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACAAAAAAATCCACATGCt                                        >  1:248883/1‑133 (MQ=17)
                                                aCCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCtg                                       <  1:303701/133‑1 (MQ=17)
                                                 ccGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACTGTCACGTACACATAAATGCACATGCtgt                                      >  2:226141/1‑133 (MQ=255)
                                                 ccGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCtgt                                      >  1:327343/1‑133 (MQ=32)
                                                   gagaTTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGa                                    <  1:320861/133‑1 (MQ=17)
                                                     gaTTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                  <  1:259940/133‑1 (MQ=32)
                                                     gaTTCCCCCAGTAGCCGCAAGCCAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                  <  1:137281/133‑1 (MQ=255)
                                                        tCCCCCAGTAACGGCAAGCGAACGCGGAGGAGGCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCg                               <  2:132579/133‑1 (MQ=255)
                                                               gTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                        <  1:265926/133‑1 (MQ=17)
                                                               gTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                        >  1:264739/1‑133 (MQ=32)
                                                                tAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGGTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGTGTAg                       >  2:369342/1‑133 (MQ=12)
                                                                  gcggcgAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAggg                     <  1:300091/133‑1 (MQ=32)
                                                                   cggcgAGCGAACGGGGAGGAGCGCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                    <  1:226141/133‑1 (MQ=17)
                                                                   cggcgAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                    <  2:327343/133‑1 (MQ=32)
                                                                     gcgagcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                  >  1:258076/1‑133 (MQ=32)
                                                                     gcgagcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                  >  1:883/1‑133 (MQ=32)
                                                                       gagcgaACGGGGAGGAGCCCAGAGGCTGGATCAGTGTGTGTATTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCTATGAGTAGGGCGGGa                <  2:210683/133‑1 (MQ=255)
                                                                        agcgaACGTGGAGGATCCTAGAACCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  1:243037/133‑1 (MQ=255)
                                                                        agcgaACGGGGAGGATCCTAGATCCTGAATTAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGACGAGTAGGGCGGGac               <  1:158322/133‑1 (MQ=11)
                                                                        agcgaACGGGGAGGAGCCCAGAGGCTGCATCAGTGGGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGAACAAAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               >  2:282312/1‑133 (MQ=11)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               >  1:258749/1‑133 (MQ=17)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  1:478163/133‑1 (MQ=17)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  2:89698/133‑1 (MQ=32)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  2:99386/133‑1 (MQ=32)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  2:356396/133‑1 (MQ=32)
                                                                        agcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac               <  1:414049/133‑1 (MQ=32)
                                                                         gcgaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGaca              >  1:40402/1‑133 (MQ=32)
                                                                           gaACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACg            <  2:33726/133‑1 (MQ=17)
                                                                                 ggaggaGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt      <  1:344418/133‑1 (MQ=32)
                                                                                  gaggagCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc     <  1:410983/133‑1 (MQ=17)
                                                                                  gaggagCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc     <  2:414977/133‑1 (MQ=32)
                                                                                  gaggagCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc     <  2:80383/133‑1 (MQ=32)
                                                                                    ggagCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCt   <  1:75573/133‑1 (MQ=31)
                                                                                     gagCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGGAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg  <  2:500487/133‑1 (MQ=11)
                                                                                                                                  |                                                                                       
CGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTG  >  NC_000913/225937‑226154

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