Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,755,8800GA31.2% 1288.4 / inf 929intergenic (+502/‑279)smpB/intAtmRNA‑binding trans‑translation protein/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base G (143/495);  new base A (51/239);  total (194/735)
Fisher's exact test for biased strand distribution p-value = 9.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 946 total aligned reads displayed.

CAAAAGAGATCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTAC  >  NC_000913_3_bme_pgi/2755781‑2755976
                                                                                                   |                                                                                                
 naaaGAGATCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATg                                                                                                  >  1:2744050/2‑99 (MQ=255)
    agagATCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATGAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTaaa                                                                                              <  1:5271398/100‑1 (MQ=255)
    agagATCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTa                                                                                                >  2:1204445/1‑98 (MQ=255)
     gagaTCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAg                                                                                             <  1:1353836/100‑1 (MQ=255)
       gaTCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGt                                                                                                 >  2:861063/1‑94 (MQ=255)
        aTCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAgact                                                                                          >  2:1996872/1‑100 (MQ=255)
        aTCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAgact                                                                                          >  2:1664212/1‑100 (MQ=255)
          cgcgTCGAAGCCCTGCCTGGGGTTGAAGCGCTACAACTTTACCAGTCTAGTTTGTTAGTGGCGAGCCCGACGGCCGCCGGCCATCGCGTGTTAAgactca                                                                                        >  2:1845062/1‑98 (MQ=255)
           gcgTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAgactgac                                                                                       >  1:4409596/1‑100 (MQ=255)
             gTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAgac                                                                                           <  1:160400/94‑1 (MQ=255)
                gAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGc                                                                                  <  1:3737781/100‑1 (MQ=255)
                     cctgcctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATgta                                                                             <  1:5192227/100‑1 (MQ=255)
                        gcctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATgtagta                                                                          >  2:447413/1‑100 (MQ=255)
                         cctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAc                                                                         <  1:4329067/100‑1 (MQ=255)
                          ctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAcc                                                                        >  1:165876/1‑100 (MQ=255)
                          ctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTAcc                                                                        >  2:218399/1‑100 (MQ=255)
                           tgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTATCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAcc                                                                        >  1:4074737/1‑99 (MQ=255)
                            ggggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAcc                                                                        >  2:4657498/1‑98 (MQ=255)
                             gggTTGAAGCGTTGAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGACAAGCGAATATAAAGACTGTCTAAGCATGTAGTACCGAg                                                                     <  1:4890604/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCgg                                                                      <  1:45507/99‑2 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCg                                                                       <  1:1998777/98‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:1724129/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:1905181/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:662130/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:1214436/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:938735/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:4136284/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATATAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:5323289/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:3814639/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:4940381/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:449647/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:1288210/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:3186978/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:2146443/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:585489/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:2087560/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAg                                                                     <  1:219529/100‑1 (MQ=255)
                             gggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGCAc                                                                         <  1:2256193/96‑1 (MQ=255)
                              ggTTGTAGCGTTAAAACGTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:4458958/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:3239532/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:1111109/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:3409117/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:514077/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:2847693/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:732935/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:4672914/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:2349003/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     >  2:2787072/1‑99 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAg                                                                     >  2:3799782/1‑99 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGacn                                                                    <  1:2676636/100‑3 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGaag                                                                    <  1:2900819/100‑3 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:3061834/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:3294465/100‑1 (MQ=255)
                              ggTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCACAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAgg                                                                    <  1:3874080/100‑1 (MQ=255)
                               gTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGa                                                                      >  2:2441302/1‑97 (MQ=255)
                               gTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGa                                                                   >  2:1059470/1‑100 (MQ=255)
                                  gAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATg                                                                 <  1:2501221/99‑1 (MQ=255)
                                   aaGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCATCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGt                                                                <  1:4017051/99‑1 (MQ=255)
                                    aGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCAAGGATGTAg                                                              >  2:4205431/1‑100 (MQ=255)
                                    aGCGTTAAAACTTAATCAGGCTACTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAg                                                              <  1:5038173/100‑1 (MQ=255)
                                     gCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAgg                                                             <  1:4738518/100‑1 (MQ=255)
                                     gCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAgg                                                             <  1:2428972/100‑1 (MQ=255)
                                      cGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAgg                                                             <  1:3672579/99‑1 (MQ=255)
                                       gTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGag                                                           <  1:4993268/100‑2 (MQ=255)
                                       gTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGaa                                                           >  2:397432/1‑100 (MQ=255)
                                       gTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGaa                                                           <  1:470440/100‑1 (MQ=255)
                                        ttAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGaag                                                          <  1:2982867/100‑2 (MQ=255)
                                            aaCTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCTAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCg                                                      >  2:4897971/1‑100 (MQ=255)
                                            aaCTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCg                                                      <  1:2578358/100‑1 (MQ=255)
                                             aCTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGGAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCgg                                                     <  2:3156890/100‑1 (MQ=255)
                                              cTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGag                                                           <  1:1173111/93‑2 (MQ=255)
                                              cTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCgg                                                     <  1:3944863/99‑1 (MQ=255)
                                              cTTAATCAGGCTAGTTTGTTACTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGaa                                                           >  2:5250576/1‑93 (MQ=255)
                                                 aaTCGGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGAc                                                   <  1:14878/98‑1 (MQ=255)
                                                 aaTCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGAc                                                   <  1:1538767/98‑1 (MQ=255)
                                                        cTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCa                                          <  1:746513/100‑1 (MQ=255)
                                                               gTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTGCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTccc                                    <  2:5057833/99‑1 (MQ=255)
                                                                ttAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCg                                   <  2:3457718/99‑1 (MQ=255)
                                                                ttAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGc                                  <  1:5234130/100‑1 (MQ=255)
                                                                   gTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTa                                                               >  2:1513432/1‑68 (MQ=255)
                                                                    tGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGc                              <  1:3499084/100‑1 (MQ=255)
                                                                    tGGCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGc                              <  1:4632099/100‑1 (MQ=255)
                                                                     ggCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCt                             >  1:262563/1‑100 (MQ=255)
                                                                      gCGTGTCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTc                            <  1:1010222/100‑1 (MQ=255)
                                                                       cGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTcc                           <  1:1438418/100‑1 (MQ=255)
                                                                         tgtCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccac                         <  1:106855/100‑1 (MQ=255)
                                                                         tgtCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccac                         <  1:4797321/100‑1 (MQ=255)
                                                                         tgtCCGTCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccac                         <  1:2780757/100‑1 (MQ=255)
                                                                         tgtCCGGCCGCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATATAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccac                         <  1:4235229/100‑1 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCTAATTTAAAGACTGACTAAGCATGTAGTACCGAGGAAGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacc                        <  1:3644743/99‑1 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacct                       <  1:1837518/100‑2 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacc                        <  1:4853697/99‑1 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacc                        <  1:4186009/99‑1 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacc                        <  1:2301134/99‑1 (MQ=255)
                                                                           tccgtccgCAGCTGGCAAGCGAATATAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTccacc                        <  1:3565944/99‑1 (MQ=255)
                                                                                         gCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCgg         <  2:2320392/100‑1 (MQ=255)
                                                                                          cAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCgg         <  1:2468140/99‑1 (MQ=255)
                                                                                          cAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGt        >  1:828456/1‑100 (MQ=255)
                                                                                            aGCGAATGTAAAGACGGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGCCGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTg       <  2:360201/99‑1 (MQ=255)
                                                                                                aaTGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTAc  >  1:2527401/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CAAAAGAGATCGCGTGGAAGCCCTGCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTAC  >  NC_000913_3_bme_pgi/2755781‑2755976

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