Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi940,7540GA21.0% 207.5 / 49.3 100intergenic (+34/‑205)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base G (37/42);  new base A (11/10);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAA  >  NC_000913_3_bme_pgi/940613‑940902
                                                                                                                                             |                                                                                                                                                    
ggTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttt                                                                                                                                              >  3:415668/1‑150 (MQ=255)
 gTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttt                                                                                                                                             >  1:43282/1‑150 (MQ=255)
 gTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttt                                                                                                                                             <  3:395812/150‑1 (MQ=255)
 gTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttt                                                                                                                                             <  4:361412/150‑1 (MQ=255)
  tCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                                            >  2:662454/1‑150 (MQ=255)
  tCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttt                                                                                                                                            <  3:65088/150‑1 (MQ=255)
  tCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGCACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttt                                                                                                                                            <  1:362313/150‑1 (MQ=255)
   cGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttttt                                                                                                                                           <  1:947655/150‑1 (MQ=255)
    gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGCCATATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                          <  2:176701/150‑1 (MQ=255)
       cGCTGGTCGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtc                                                                                                                                       <  3:489161/150‑1 (MQ=255)
            gTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                                                                    <  1:293694/148‑1 (MQ=255)
            gTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                                                                    >  2:293694/1‑148 (MQ=255)
             ttGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCTCCTGCGGGCGCTTTTTTTGTCTCCCtt                                                                                                                                 <  1:771186/150‑1 (MQ=255)
              tGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                                                                                <  4:291897/150‑1 (MQ=255)
               gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                                                                               <  4:415668/150‑1 (MQ=255)
               gCAGTAATGGAAAACCATCCGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTCATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                                                                               <  4:309503/150‑1 (MQ=255)
                  gTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                                                                                <  4:528586/146‑1 (MQ=255)
                  gTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                                                                                >  3:528586/1‑146 (MQ=255)
                        gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                                                                      <  4:241050/150‑1 (MQ=255)
                         aaaaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                          >  2:935820/1‑129 (MQ=255)
                         aaaaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                          <  1:935820/129‑1 (MQ=255)
                                tCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                                                                              >  1:112956/1‑150 (MQ=255)
                                     agGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                         <  2:1010240/150‑1 (MQ=255)
                                       gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                                       >  3:472571/1‑150 (MQ=255)
                                       gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCAGTATCTGCACCCAATGCAACATAATCAATAAAAAATAATGTCTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                                       <  1:640290/150‑1 (MQ=255)
                                           aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                                            <  2:13864/109‑1 (MQ=255)
                                           aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                                            >  1:13864/1‑109 (MQ=255)
                                            tGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTac                                                                                                  <  3:193719/150‑1 (MQ=255)
                                               tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                                                                               <  2:43282/150‑1 (MQ=255)
                                               tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                                                                               >  4:156686/1‑150 (MQ=255)
                                                cGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGt                                                                                              >  4:427309/1‑150 (MQ=255)
                                                cGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGt                                                                                              >  3:116979/1‑150 (MQ=255)
                                                  tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                            <  1:720492/150‑1 (MQ=255)
                                                          tACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGCATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACtt                                                                                    >  4:641183/1‑150 (MQ=255)
                                                              aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                                                                                   <  3:175658/115‑1 (MQ=255)
                                                              aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                                                                                   >  4:175658/1‑115 (MQ=255)
                                                              aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                <  3:641183/150‑1 (MQ=255)
                                                                  gTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                            <  4:472571/150‑1 (MQ=255)
                                                                      tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAAt                                                                        <  3:205339/150‑1 (MQ=255)
                                                                         gTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                                                                     >  3:189206/1‑150 (MQ=255)
                                                                          tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                    <  2:112956/150‑1 (MQ=255)
                                                                             gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa                                                                 >  2:995505/1‑150 (MQ=255)
                                                                             gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa                                                                 >  3:439590/1‑150 (MQ=255)
                                                                                   gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCtt                                                           >  2:798410/1‑150 (MQ=255)
                                                                                     aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                                                                             >  3:611948/1‑98 (MQ=255)
                                                                                     aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                                                                             <  4:611948/98‑1 (MQ=255)
                                                                                         aCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                       <  3:115877/84‑1 (MQ=255)
                                                                                         aCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                       >  4:115877/1‑84 (MQ=255)
                                                                                          cTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtttt                                                    <  1:995505/150‑1 (MQ=255)
                                                                                           tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTc                                                   <  1:798410/150‑1 (MQ=255)
                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGTCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  >  3:335599/1‑117 (MQ=255)
                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGTCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  <  4:335599/117‑1 (MQ=255)
                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                                 >  2:923462/1‑150 (MQ=255)
                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                                 >  3:149233/1‑150 (MQ=255)
                                                                                                tatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                            <  2:40744/120‑1 (MQ=255)
                                                                                                tatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                            >  1:40744/1‑120 (MQ=255)
                                                                                                   atTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattatt                                           <  2:308738/150‑1 (MQ=255)
                                                                                                      ggCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGa                                        <  4:439590/150‑1 (MQ=255)
                                                                                                          aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                                    <  4:65365/150‑1 (MQ=255)
                                                                                                              cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGC‑TTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACaa                               >  3:462125/1‑150 (MQ=255)
                                                                                                               ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                <  4:452495/101‑1 (MQ=255)
                                                                                                               ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                >  3:452495/1‑101 (MQ=255)
                                                                                                               ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACaa                               >  4:533193/1‑150 (MQ=255)
                                                                                                               ccAATTGTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACaa                               >  2:412855/1‑150 (MQ=255)
                                                                                                                         tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                                       >  3:638718/1‑68 (MQ=255)
                                                                                                                         tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                                       <  4:638718/68‑1 (MQ=255)
                                                                                                                            acaaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                    <  1:453367/96‑1 (MQ=39)
                                                                                                                            acaaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                    >  2:453367/5‑100 (MQ=39)
                                                                                                                             tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa                 >  2:510936/1‑150 (MQ=255)
                                                                                                                              aaaaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                   <  4:546530/145‑1 (MQ=255)
                                                                                                                              aaaaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                   >  3:546530/3‑147 (MQ=255)
                                                                                                                              aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  >  4:385673/3‑84 (MQ=38)
                                                                                                                              aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  <  3:385673/82‑1 (MQ=38)
                                                                                                                              aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGGAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                >  4:325541/3‑150 (MQ=255)
                                                                                                                              aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                          >  4:69355/3‑124 (MQ=255)
                                                                                                                              aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                          <  3:69355/122‑1 (MQ=255)
                                                                                                                               tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                       >  4:133896/1‑46 (MQ=255)
                                                                                                                               tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                       <  3:133896/46‑1 (MQ=255)
                                                                                                                               aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                    <  4:281690/96‑1 (MQ=255)
                                                                                                                               aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                    >  3:281690/2‑97 (MQ=255)
                                                                                                                               aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt                                                       >  2:512300/2‑110 (MQ=255)
                                                                                                                               aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt                                                       <  1:512300/109‑1 (MQ=255)
                                                                                                                                aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                          <  2:505886/106‑1 (MQ=255)
                                                                                                                                aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                          >  1:505886/1‑106 (MQ=255)
                                                                                                                                aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                           <  1:658301/137‑1 (MQ=255)
                                                                                                                                aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                           >  2:658301/1‑137 (MQ=255)
                                                                                                                                aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTCTATAGCGAACTg              >  4:318794/1‑150 (MQ=255)
                                                                                                                                aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTg              >  1:1011315/1‑150 (MQ=255)
                                                                                                                                aaaaaGCGCCAGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTCCCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTg              >  2:47688/1‑150 (MQ=255)
                                                                                                                                 aaaaGCGCCCGCAGGCGCTTTTTTTTTTTTCTTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc             >  4:430259/1‑150 (MQ=255)
                                                                                                                                 aaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                                 >  3:173204/1‑114 (MQ=255)
                                                                                                                                 aaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                                 <  4:173204/114‑1 (MQ=255)
                                                                                                                                   aaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                          >  3:97066/1‑119 (MQ=255)
                                                                                                                                   aaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                          <  4:97066/119‑1 (MQ=255)
                                                                                                                                   aaGCGCCTGCGGGCGC‑TTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCtat          <  4:462125/150‑1 (MQ=255)
                                                                                                                                    aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                                                                                >  2:526705/1‑64 (MQ=255)
                                                                                                                                    aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                                                                                <  1:526705/64‑1 (MQ=255)
                                                                                                                                        ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  <  4:394291/74‑1 (MQ=255)
                                                                                                                                        ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                  >  3:394291/1‑74 (MQ=255)
                                                                                                                                            cGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAaataa  <  4:116979/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAA  >  NC_000913_3_bme_pgi/940613‑940902

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