Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi892,4550GC20.0% 199.3 / 47.4 81V163V (GTG→GTC) rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base C (7/9);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAA  >  NC_000913_3_bme_pgi/892308‑892604
                                                                                                                                                    |                                                                                                                                                     
gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAt                                                                                                                                                     <  3:166840/150‑1 (MQ=255)
gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAt                                                                                                                                                     >  2:43952/1‑150 (MQ=255)
 cctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                    <  4:282561/150‑1 (MQ=255)
   tgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGa                                                                                                                                                  >  2:243686/1‑150 (MQ=255)
    gTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGAc                                                                                                                                                 <  3:55572/150‑1 (MQ=255)
       aCGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCt                                                                                                                                              >  3:206448/1‑150 (MQ=255)
         gggCATTGCGCATTCGCCCGATGATACCAGCGATTTAATCGACATGGACGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCttt                                                                                                                                            >  2:529914/1‑150 (MQ=255)
           gCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCTTTcc                                                                                                                                          <  1:741082/150‑1 (MQ=255)
           gCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCTTTcc                                                                                                                                          <  4:381562/150‑1 (MQ=255)
              ttGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgcg                                                                                                                                       >  3:334401/1‑150 (MQ=255)
              ttGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgcg                                                                                                                                       <  4:162194/150‑1 (MQ=255)
                gcgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGt                                                                                                                                     <  2:427505/150‑1 (MQ=255)
                 cgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                    <  2:745921/134‑1 (MQ=255)
                 cgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                    >  1:745921/1‑134 (MQ=255)
                 cgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTc                                                                                                                                    <  3:306607/150‑1 (MQ=255)
                     ttCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGaa                                                                                                                                >  2:367936/1‑150 (MQ=255)
                     ttCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGaa                                                                                                                                >  4:244428/1‑150 (MQ=255)
                         ccGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCTTTcc                                                                                                                                          >  2:75130/1‑136 (MQ=255)
                         ccGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCTTTcc                                                                                                                                          <  1:75130/136‑1 (MQ=255)
                               gatACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTTAACGCGCATAtt                                                                                                                      <  1:87922/150‑1 (MQ=255)
                                  aCCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGGTGAAGGCCCGCAGGGAACTGGCGTCGTGCTGGAGGAGACGCG‑TCAGGCGGCGGAAAGCGTGGTTGGCGCGTTCCCCGGTCTGCACGCGCATATTCTg                                                                                                                   >  2:61098/1‑150 (MQ=255)
                                      gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCCGCGGAAAGCGTCAATCACGCTTTCcgccgc                                                                                                                                     >  4:358130/1‑125 (MQ=38)
                                      gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCCGCGGAAAGCGTCAATCACGCTTTCcgccgc                                                                                                                                     <  3:358130/128‑4 (MQ=38)
                                        gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                             <  3:49128/150‑1 (MQ=255)
                                          tttAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                    >  2:236479/1‑122 (MQ=25)
                                          tttAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                    <  1:236479/126‑5 (MQ=25)
                                               tCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatata                                                                                                      >  2:274110/1‑150 (MQ=255)
                                                cGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                    >  1:522724/1‑116 (MQ=25)
                                                cGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                    <  2:522724/120‑5 (MQ=25)
                                                cGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatatat                                                                                                     <  1:43952/150‑1 (MQ=255)
                                                  aCATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcgcg                                                                                                                            <  1:521529/125‑1 (MQ=255)
                                                  aCATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcgcg                                                                                                                            >  2:521529/1‑125 (MQ=255)
                                                    aTGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaaa                                                                                                 <  1:458967/150‑1 (MQ=255)
                                                    aTGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTCTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaaa                                                                                                 <  2:584106/150‑1 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgcg                                                                                                                                       >  1:311553/1‑109 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgcg                                                                                                                                       <  2:311553/109‑1 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAagag                                                                                              <  4:367320/150‑1 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGACCGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGacgc                                                                                                                              <  2:809910/118‑4 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGACCGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGacgc                                                                                                                              >  1:809910/1‑115 (MQ=255)
                                                         cGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGgc                                                                                            <  4:334401/150‑1 (MQ=255)
                                                              gtccgccgcTGGTGGTCCCGTTGGTTGCAGGCACGCAGGGAATTGGCGTCGTCCTTGCGCAGACGCG‑TCAGGCGGCGGAACGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCCTATTCTGGTGCCGGAATATATCAAAGAGGCGCAAg                                                                                       <  2:641895/147‑1 (MQ=255)
                                                                  gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                               <  1:672633/122‑1 (MQ=255)
                                                                  gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                               >  2:672633/1‑122 (MQ=255)
                                                                  gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                               >  2:672603/1‑122 (MQ=255)
                                                                  gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                               <  1:672603/122‑1 (MQ=255)
                                                                         gtggtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATc                                                                            <  3:244428/150‑1 (MQ=255)
                                                                            gtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGc                                                                         >  3:182202/1‑150 (MQ=255)
                                                                               aaGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACgcgc                                                                                                                           >  1:104382/1‑97 (MQ=255)
                                                                               aaGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACgcgc                                                                                                                           <  2:104382/97‑1 (MQ=255)
                                                                                  ttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGAGTAAATGACTCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTg                                                                   >  3:403747/1‑150 (MQ=255)
                                                                                     gttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGTTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgtt                                                                <  1:529914/150‑1 (MQ=255)
                                                                                      ttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttg                                                               >  1:222742/1‑150 (MQ=255)
                                                                                      ttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttg                                                               >  2:627625/1‑150 (MQ=255)
                                                                                        gAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATAGCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgtt                                                             >  4:461344/1‑150 (MQ=255)
                                                                                          aGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTgg                                                           <  4:390177/150‑1 (MQ=255)
                                                                                                    ggAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAATGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGAAAGAGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTg                                                 <  1:243686/150‑1 (MQ=255)
                                                                                                           gCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTg                                          >  4:250656/1‑150 (MQ=255)
                                                                                                                  gCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGa                                   <  2:222742/150‑1 (MQ=255)
                                                                                                                    tGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAAc                                 >  2:644475/1‑150 (MQ=255)
                                                                                                                      gCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                 >  2:763703/1‑116 (MQ=255)
                                                                                                                      gCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                 <  1:763703/116‑1 (MQ=255)
                                                                                                                        ggAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAAcggcg                             >  2:283229/1‑150 (MQ=255)
                                                                                                                           gaCGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACggcgggc                          <  3:250656/150‑1 (MQ=255)
                                                                                                                              gcg‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGt                                               >  3:454558/1‑126 (MQ=255)
                                                                                                                              gcg‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGt                                               <  4:454558/126‑1 (MQ=255)
                                                                                                                                    aGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAgg                  <  3:234105/150‑1 (MQ=255)
                                                                                                                                    aGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAgg                  >  4:322989/1‑150 (MQ=255)
                                                                                                                                      gcggcgGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGc                >  1:678541/1‑150 (MQ=255)
                                                                                                                                         gcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGa                                                        <  3:331948/107‑1 (MQ=39)
                                                                                                                                         gcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGa                                                        >  4:331948/1‑107 (MQ=39)
                                                                                                                                          cggAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGaa                                  >  2:356320/1‑128 (MQ=255)
                                                                                                                                          cggAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGaa                                  <  1:356320/128‑1 (MQ=255)
                                                                                                                                          cggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                 <  1:786809/145‑1 (MQ=255)
                                                                                                                                          cggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                 >  2:786809/1‑145 (MQ=255)
                                                                                                                                          cggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAtt            >  3:366288/1‑150 (MQ=255)
                                                                                                                                             aaaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTc         <  1:283229/150‑1 (MQ=255)
                                                                                                                                               aGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGt       <  4:366288/150‑1 (MQ=255)
                                                                                                                                                 cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc     <  2:832222/150‑1 (MQ=255)
                                                                                                                                                 cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc     >  1:832222/1‑150 (MQ=255)
                                                                                                                                                   tGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTgg                                                           >  4:182891/1‑94 (MQ=255)
                                                                                                                                                   tGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTgg                                                           <  3:182891/94‑1 (MQ=255)
                                                                                                                                                    gATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCTCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCaa  >  1:836099/1‑150 (MQ=255)
                                                                                                                                                    gATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCaa  <  4:182202/150‑1 (MQ=255)
                                                                                                                                                    gATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCaa  <  1:367936/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAA  >  NC_000913_3_bme_pgi/892308‑892604

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