Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000963,945,8190GA20.0% 110.3 / 10.1 50noncoding (2116/2904 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (15/25);  new base A (10/0);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 6.36e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCC  >  U00096/3945759‑3945967
                                                            |                                                                                                                                                    
ggAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg                                                              <  2:2892499/149‑1 (MQ=31)
ggAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg                                                              <  2:2892498/149‑1 (MQ=31)
 gAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg                                                              <  1:2892500/148‑1 (MQ=31)
 gAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAAATCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGa                                                             >  2:2892560/1‑149 (MQ=11)
  aaaGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGCCCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAc                                                            >  1:2892591/1‑149 (MQ=17)
   aaGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTTTAAAAACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAcc                                                           >  2:2892549/1‑149 (MQ=11)
    aGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAccc                                                          <  2:2892501/149‑1 (MQ=32)
    aGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAccc                                                          >  1:2892561/1‑149 (MQ=32)
      cccccGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAAGGTGGGAGGCATTGAAGTGATGACACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGt                                                        >  2:2892564/2‑149 (MQ=11)
       ccccGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTg                                                       <  1:2892502/149‑1 (MQ=32)
          cGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATc                                                    <  1:2892504/149‑1 (MQ=34)
          cGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATc                                                    <  2:2892503/149‑1 (MQ=34)
           gTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAAGTCTGCATGGAGACGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATcc                                                   >  2:2892580/1‑149 (MQ=11)
            tGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAACATTAAAAAAAAAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCg                                                  >  2:2892610/1‑149 (MQ=11)
              aaCCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATcc                                                   <  1:2892505/146‑1 (MQ=34)
               aCCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGt                                               <  2:2892506/149‑1 (MQ=34)
                ccTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGtt                                              <  1:2892508/149‑1 (MQ=34)
                ccTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGtt                                              <  2:2892507/149‑1 (MQ=34)
                   ttACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAACTCTGCATGGAGCCGACCTAGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCg                                           >  2:2892563/1‑149 (MQ=11)
                     aCTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGa                                         <  2:2892509/149‑1 (MQ=34)
                     aCTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGa                                         <  2:2892510/149‑1 (MQ=34)
                       tataGCTTGCCACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCATTTAATGTTTGATGTTCTAACGTGGACCCGAGATCCGGGTTGCGGACa                                       >  1:2892530/1‑149 (MQ=255)
                       tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGCCCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa                                       >  1:2892576/1‑149 (MQ=21)
                       tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa                                       >  1:2892567/1‑149 (MQ=34)
                       tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTCAATTGTGAACACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa                                       >  2:2892579/1‑149 (MQ=11)
                       tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa                                       >  2:2892577/1‑149 (MQ=11)
                        ataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTTAAATTTAAACTACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAg                                      >  2:2892574/1‑149 (MQ=11)
                         taGCTTGACACTGTACTTTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGAGGCCAGTCGGCATGGAGCCGACCTTGAAATACGTCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAgt                                     <  1:2892511/149‑1 (MQ=11)
                         taGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTATGCATGGAGCCGACCTTGAAATACCCCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAg                                      >  1:2892559/1‑148 (MQ=21)
                          aGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAgtg                                    <  2:2892513/149‑1 (MQ=34)
                            cTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTc                                  >  1:2892627/1‑149 (MQ=34)
                             ttGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCt                                 <  2:2892515/149‑1 (MQ=34)
                             ttGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCt                                 <  2:2892514/149‑1 (MQ=34)
                              tGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCtg                                >  2:2892566/1‑149 (MQ=11)
                                 cacTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCtggtg                             <  1:2892516/149‑1 (MQ=34)
                                   cTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTggg                           <  1:2892518/149‑1 (MQ=34)
                                    tGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGt                          <  2:2892519/149‑1 (MQ=34)
                                      aaCATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCCTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTa                         >  1:2892552/1‑148 (MQ=21)
                                       aCATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGt                       >  1:2892541/1‑149 (MQ=34)
                                        cATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGtt                      <  2:2892520/149‑1 (MQ=34)
                                         aTTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGttt                     <  1:2892521/149‑1 (MQ=34)
                                          ttGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTg                    >  1:2892638/1‑149 (MQ=34)
                                              ggcgtGATGGCTAGGATAGGTGGGAGGCTTTGAAGTCTTGACGCAGGTCTGCAGGGAGCCGACCTTGAAAAACCCCCCTTTCATGTTTGCAGATCGCACGTGGCCCCGTGATCCGGGTTGAGGACAGTGACTGCTGGGTAGATTGACTg                >  1:2892654/5‑149 (MQ=2)
                                               ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGTCGCCAGGCTGCTTGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACTTGGTCCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTgg               <  2:2892522/149‑1 (MQ=11)
                                               ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCCCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGAc                  >  1:2892592/1‑146 (MQ=21)
                                               ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTGGAAATACAACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTgg               >  1:2892635/1‑149 (MQ=21)
                                                cTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTggg              <  1:2892523/149‑1 (MQ=34)
                                                cTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTggg              >  1:2892584/1‑149 (MQ=34)
                                                     tgtgTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGt         >  1:2892588/1‑149 (MQ=34)
                                                          aGGATAGGGGGGAGGCTTTGAAGTGGGGACGCCAGGCGGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGAGCCGGGTTGCGGACAGTGGCTGGTGGGTAGTTTGACTGGGGCGGtctc      <  1:2892525/147‑1 (MQ=11)
                                                            gATAGGTGGGTGGCTTTGAAGTGTGGACGCCAGTCTGGATGGAGCCTACCTTGAAATACCACTCTTTAATGTTTGTTGTTCTAACGTGGACCCGTGATCCGGGTTGCGGTCAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCtcctcc  <  2:2892526/149‑1 (MQ=11)
                                                            |                                                                                                                                                    
GGAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCC  >  U00096/3945759‑3945967

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