Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 569,862 G→A 5.8% V321I (GTA→ATA)  ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi569,8620GA 5.8% 176.1 / 2.9 71V321I (GTA→ATA) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base A (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

GATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAAC  >  NC_000913_3_bme_pgi/569721‑570008
                                                                                                                                             |                                                                                                                                                   
gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTgg                                                                                                                                              >  2:548895/1‑149 (MQ=255)
 atgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAACAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGa                                                                                                                                             >  1:796866/1‑149 (MQ=255)
  tgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGCAAATAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGaa                                                                                                                                            >  1:663010/1‑149 (MQ=255)
  tgagTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATACTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGaa                                                                                                                                            <  2:791107/149‑1 (MQ=255)
        gTTTTACGCTGTACAAGAAATTCGTTTGGCACCTTTTTGTAATAACAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGc                                                                                                                                     <  1:672151/149‑1 (MQ=255)
         ttttACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCa                                                                                                                                    <  1:529308/149‑1 (MQ=255)
             aCGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGc                                                                                                                                >  1:256064/1‑149 (MQ=255)
              cGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCATATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTATGGTAATACCATGATTGCTCATTCGGTATCTGGAGGATTGCATGACt                                                                                                                               >  2:465200/1‑149 (MQ=255)
               gCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCCTGGCta                                                                                                                              >  2:350203/1‑149 (MQ=255)
                   tACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACTATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtattat                                                                                                                          >  2:62438/1‑149 (MQ=255)
                   tACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtattat                                                                                                                          >  1:688189/1‑149 (MQ=255)
                       ggAAATTCGGTTGGACACTGTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTTATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                                      <  2:842149/149‑1 (MQ=255)
                       ggAAATTCGGTGTGCACCTTGTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAATTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                                      <  2:157522/149‑1 (MQ=255)
                           aTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTcc                                                                                                                  <  2:543301/149‑1 (MQ=255)
                                tgtGGTACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATga                                                                                                             <  1:610337/147‑1 (MQ=255)
                                gTTGGCAAGTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATga                                                                                                             <  1:163942/149‑1 (MQ=255)
                                    gCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGAGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaag                                                                                                         >  1:40673/1‑149 (MQ=255)
                                     cACTTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaaga                                                                                                        <  1:424408/149‑1 (MQ=255)
                                       ccTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTAAGATTGGTTGTGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGAtt                                                                                                      >  1:173262/1‑149 (MQ=255)
                                         ttttGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTAc                                                                                                    >  2:796934/1‑149 (MQ=255)
                                           ttGGTATTAACAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACAt                                                                                                  <  2:796866/149‑1 (MQ=255)
                                               tATTAGCAATAGTAGCATGAATCCTATTCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                                              <  1:240329/149‑1 (MQ=255)
                                               tATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                                              <  2:348683/149‑1 (MQ=255)
                                               tATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCAACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCTTGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                                              >  2:791094/1‑149 (MQ=255)
                                                 ttAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATAACATCGAtgt                                                                                            >  1:855422/1‑149 (MQ=255)
                                                 ttAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTCTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAtgt                                                                                            >  1:893431/1‑149 (MQ=255)
                                                     cAATAGTATCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACa                                                                                        <  2:350363/149‑1 (MQ=255)
                                                     cAATAGTAGCAAGAATCCTATTCTAATAAATCTACTTCGAACAGTTAAGAAGTGTGAGGCTTGTGGTAATACCGTGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGCGAAGATTACATCGATGTGACa                                                                                        >  2:386250/1‑149 (MQ=255)
                                                     cAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACa                                                                                        >  1:531600/1‑149 (MQ=255)
                                                     cAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACa                                                                                        >  1:779545/1‑149 (MQ=255)
                                                        tagtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt                                                                                                                       >  1:228117/1‑115 (MQ=255)
                                                        tagtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt                                                                                                                       <  2:228117/115‑1 (MQ=255)
                                                           tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGAACGGCCAt                                                                                  >  2:768099/1‑149 (MQ=255)
                                                             gCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCa                                                                                >  1:552775/1‑149 (MQ=255)
                                                             gCAAGAATCCTATGCTAATAAATCTACTTCGAAAAGTTATGAAGTGTGAGGCCTGTGGTAATACCATGATCGTTCATGCGGTATCTAGAAGTTTGCATGGCTATTATGTTTGTCCTATTAGAAGATTACATAGAAGTGGGAGGccctga                                                                                >  2:852289/1‑145 (MQ=255)
                                                               aaGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtatta                                                                                                                           >  1:303723/1‑104 (MQ=255)
                                                               aaGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtatta                                                                                                                           <  2:303723/104‑1 (MQ=255)
                                                                    tCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaaaaa                                                                         >  1:313714/1‑149 (MQ=255)
                                                                       tATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATAAATAAAAagag                                                                      >  1:259293/1‑149 (MQ=255)
                                                                            taataaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTg                                                                 >  2:723679/1‑149 (MQ=255)
                                                                              ataaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGAAATCAATAAAAAGAGATTTGGt                                                               >  1:553737/1‑149 (MQ=255)
                                                                               taaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGtt                                                              >  1:673736/1‑149 (MQ=255)
                                                                               taaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTAACAGGCCATCAATAAAAAGAGATTTGGtt                                                              >  1:873679/1‑149 (MQ=255)
                                                                                    cTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATta                                                         >  1:687411/1‑149 (MQ=255)
                                                                                      aCTGCGAACAGTTATGAAGTGGTAGGCTTGCGGTAATACCATGATTGTTCATGCGGTATGTGGAAGTTTGCATGGCTATTATGTTTGTCCTATGAGAAGATTACATAGATGTGACAGGCTATCAATAAAAAGAGATTTGGTTGATtata                                                       <  2:269242/149‑1 (MQ=255)
                                                                                            aataGTTATGAAGTGGTAGTATTGATGTAATACCATGATTGTTTATGCGGTATCTTGAAGTTTGCATGGCTATTATGTTTGTCCTATGAGAAGATTACATCGATGTGACAGGTCATCAATAAAAAGAGATTTGGTTGATTATAATATCa                                                 <  2:296734/146‑1 (MQ=255)
                                                                                             aCAGTTATGAAGTGTGAGGGTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCAt                                                <  2:590559/149‑1 (MQ=255)
                                                                                                 ttATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGACCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAAt                                            >  1:453442/1‑149 (MQ=255)
                                                                                                 taaggAAGTGTGAGTGTTGGGGTAATACCATGATTGTTGATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCAAACAATAAAAAGAGATTTGGTTGATTATAATATCATTAAt                                            <  2:332067/145‑1 (MQ=255)
                                                                                                    tGAAGTGTGAGTCTTGTGGTAATATCATGATTTTTCATGTGGTATATGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGaa                                         <  1:136252/149‑1 (MQ=255)
                                                                                                      aaGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTTTGTGGCGGTCTATTACCGGTGCCTGTTGGTCCGGTGCGAGGATTGCCGCGAGGGGGCTGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAAtt                                       <  2:639351/149‑1 (MQ=255)
                                                                                                         tgtgAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCACTAAAAAAAGATTTTGTTGATTATAATATCATTAATGAATTGCt                                    >  2:651831/1‑149 (MQ=255)
                                                                                                         tgtgAGGCGTGTGGTAATAATATGATTGTTCATGGGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCt                                    <  2:634335/149‑1 (MQ=255)
                                                                                                             aGGCTTGTGGTAATACCATGATTGTTCATGCGATATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGTCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                                >  2:484767/1‑149 (MQ=255)
                                                                                                             aGGCTTGTGGTAATACCATGATTGTTCATGCGATATCTGGAAGTTTGCATGGCTATTATGTTCGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                                >  2:840225/1‑149 (MQ=255)
                                                                                                                 ttGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAAAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAAtt                            >  1:366953/1‑149 (MQ=255)
                                                                                                                   ttgGTAATACCATTATTGTTTATGCCGTATCTGGAAGTCTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGATATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGt                          <  2:893431/148‑1 (MQ=255)
                                                                                                                    tgttAATCCAATGAGTGTTTATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCTATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTa                         <  1:229693/146‑1 (MQ=255)
                                                                                                                     ggTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTACGATGAGAAGATTGCATCGATGTGACAGCAAACGTAAAACAGGAAGTTAGGGGGATTATAATATCATTAATGAATTGCTTTTTAATTGTAg                        >  2:411556/1‑149 (MQ=255)
                                                                                                                       tAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACCTCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATAATTAATGAATTGCTTTTTAATTGTAGCa                      >  1:194066/1‑149 (MQ=255)
                                                                                                                              aTGATTGTTCATGCGGTGTGTGGAAGTTTGCATGGGGATTGTGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCa               >  2:666775/1‑149 (MQ=255)
                                                                                                                               tgcttGGTCAAGCGGGATCTGGAAGGTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATTGATGTGACAGGTCATCCATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCaa              <  1:350203/146‑1 (MQ=255)
                                                                                                                               tGATTTTTCATGCGTTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATCAATGAATTGCTTTTTAATTGTAGCAAAATTCaa              <  1:636854/149‑1 (MQ=255)
                                                                                                                               tGATTGTTCATGCGATATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGTCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCaa              <  1:484767/149‑1 (MQ=255)
                                                                                                                               tGATTGTTCATGCGATATCTGGAAGTTTGCATGGCTATTATGTTCGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCaa              <  1:840225/149‑1 (MQ=255)
                                                                                                                                gATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAAc             <  1:58169/149‑1 (MQ=255)
                                                                                                                                gATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAATATAATATCATTAATGAATTGCTATTTAACAGTAGCAAAATTCAAc             >  1:837289/1‑149 (MQ=255)
                                                                                                                                  tgtgtCATGGGCTATGTGGAATGTTCGATGGCTAGTATGTTTGCGCGATGAGAAGATTAGTTGGATGTGACAGGCCATCCATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATGCAACCa           <  2:167493/145‑1 (MQ=255)
                                                                                                                                    gTTCTTGCGGTATTTGGATGTTTGCTTGGCTTTTTTGTTTGTCGGGTGAGAGGGTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGt         <  1:116884/149‑1 (MQ=255)
                                                                                                                                           cGGTATCTGGAAGTTTGCATTGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAAc  <  1:314900/149‑1 (MQ=255)
                                                                                                                                           cGGTATCTGGAAGTTTGAATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCTATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAAc  <  1:233376/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAAC  >  NC_000913_3_bme_pgi/569721‑570008

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