Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,950,2920CG20.4% 404.8 / 106.4 178A77A (GCG→GCC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (67/73);  new base G (22/14);  total (91/87)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 184 total aligned reads displayed.

CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCA  >  NC_000913_3_bme_pgi/1950196‑1950389
                                                                                                 |                                                                                                  
cGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCa                                                                                                 <  1:767972/100‑1 (MQ=255)
 ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCa                                                                                                 <  2:1461365/99‑1 (MQ=255)
  gCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑Cg                                                                                              <  1:4799393/100‑1 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTa                                                                                            <  1:5066653/100‑1 (MQ=255)
     gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTAc                                                                                           >  1:3807122/1‑100 (MQ=255)
       aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACgg                                                                                         <  2:1718679/100‑1 (MQ=255)
        cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGt                                                                                        >  1:1368095/1‑100 (MQ=255)
           cTGAGCGGTTGATGATAGTCAGCGAGGACGCCAGCTCTTCGATTACGCCGGTCGCCATATCGCGGCGAATA‑TTTTTTTCGCCACGCGCA‑CGTCCGGTgac                                                                                     >  2:4430605/1‑98 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                    <  1:1863547/98‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                    <  2:1844894/98‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                     <  2:1502526/97‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                   <  2:2662691/99‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTGCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                  <  2:3595270/100‑1 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                  >  2:3922894/1‑99 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                  >  2:3743143/1‑99 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                 <  1:4278579/99‑1 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                 <  2:73124/99‑1 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                <  1:1811756/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                               <  1:3751942/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                               <  1:1012636/100‑1 (MQ=38)
                  gTTGATGATAGTCAGCGGGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                              <  2:4990087/100‑1 (MQ=255)
                  gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                               <  2:4508574/99‑1 (MQ=255)
                  gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGATGCCCGTGAcc                                                                              >  2:3406263/1‑100 (MQ=255)
                  gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTCAATA‑TTTTTTTCGTCACGGTCACCGTAC‑GTGCGCGTGAc                                                                               <  1:4892189/99‑1 (MQ=255)
                   ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                             <  1:3295370/100‑3 (MQ=37)
                    tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCt                                                                             >  2:1799604/1‑99 (MQ=255)
                     gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                            >  1:1728073/1‑99 (MQ=255)
                     gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                            <  1:4062898/99‑4 (MQ=37)
                      atgatAGTCAGCGAGGGCGCCAGCACGTCGATTTCGCCGGGCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTACGG‑GCGCGTGAcgaaaa                                                                          <  2:2817349/100‑6 (MQ=25)
                      atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                           <  1:374338/99‑5 (MQ=37)
                      atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                           <  1:1547283/99‑5 (MQ=37)
                      atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTGTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                               <  2:1502118/95‑1 (MQ=255)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                         >  1:2133471/1‑100 (MQ=255)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                         >  1:2404116/1‑100 (MQ=255)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                         >  1:346290/1‑100 (MQ=255)
                         atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                       <  2:2577005/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCa                                                                      <  1:4129102/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCa                                                                      <  1:262410/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                     >  2:2309956/1‑99 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                     >  1:1268785/1‑99 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                     >  1:1526134/1‑99 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                     >  1:4604269/1‑99 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                     <  1:1654758/99‑1 (MQ=255)
                             tCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                    <  1:2603663/99‑1 (MQ=255)
                               aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                  <  1:3277551/99‑1 (MQ=255)
                                 cGAGGGCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                           <  2:4297280/88‑5 (MQ=21)
                                 cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                               >  2:4764104/1‑100 (MQ=255)
                                   aGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                 >  2:3519413/1‑96 (MQ=255)
                                    ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                              >  2:4799371/1‑98 (MQ=255)
                                        gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                        <  1:3990142/100‑1 (MQ=255)
                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                       >  2:1606319/1‑100 (MQ=255)
                                           aGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGtt                                                     >  1:3855710/1‑100 (MQ=255)
                                              aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                  <  2:3121772/100‑1 (MQ=255)
                                               cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCa                                                   >  1:92403/1‑98 (MQ=255)
                                                  cGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                               <  1:3567126/99‑1 (MQ=255)
                                                  cGATTTCGCCGGTCGCCATAACGCGGTTAATATTTTTTTTCGTCACGAGCA‑CGTAAGGTGCCCGTGACCTGAACGCAGGACTCATTACGCAGTTCACAgg                                               >  2:3253560/1‑100 (MQ=255)
                                                     tttCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                           <  1:3345076/100‑1 (MQ=255)
                                                       tCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                          <  1:864868/99‑1 (MQ=255)
                                                         gCCGGTCGCCATATCGCGGTTAATA‑TTTTTGTCGGCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                 <  2:1439381/90‑1 (MQ=255)
                                                          ccGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                       >  1:439497/1‑99 (MQ=255)
                                                          ccGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                      <  1:1251433/100‑1 (MQ=255)
                                                            ggTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                     >  1:4023620/1‑99 (MQ=255)
                                                             gTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                   <  1:5000641/100‑1 (MQ=255)
                                                             gTCGCCATATCGCGGTTAATA‑TTTTTTCCGTCACGCGCC‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                   >  1:4052050/1‑100 (MQ=255)
                                                                gCCATATCGCGGTTAATA‑TTTTTTTTTTTTTTTTTC‑CGTACGGAGCCCGTGACCCGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTCACGCGt                                >  2:1145143/1‑100 (MQ=255)
                                                                gCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                 <  1:1004543/99‑1 (MQ=255)
                                                                gCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGAGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                >  1:5276951/1‑100 (MQ=255)
                                                                 ccATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                >  2:4607545/1‑99 (MQ=255)
                                                                 ccATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACCGTGCCCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCCTTAACGCGtt                               >  2:3218259/1‑99 (MQ=255)
                                                                  cATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                              >  1:3884936/1‑100 (MQ=255)
                                                                  cATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                              >  1:2358096/1‑100 (MQ=255)
                                                                    tatCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                             <  2:1722966/99‑1 (MQ=255)
                                                                      tCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                          <  2:414275/100‑1 (MQ=255)
                                                                        gcgGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                        >  1:5234157/1‑100 (MQ=255)
                                                                           gTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                     >  1:4073840/1‑100 (MQ=255)
                                                                            ttAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCg                    >  2:533338/1‑100 (MQ=255)
                                                                              aaTA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                  >  2:832057/1‑100 (MQ=255)
                                                                                ta‑tTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                <  1:4619548/100‑1 (MQ=255)
                                                                                   tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg               <  1:1298831/99‑1 (MQ=255)
                                                                                   tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg               <  1:2903591/99‑1 (MQ=255)
                                                                                    ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa              >  2:2197299/1‑99 (MQ=255)
                                                                                      ttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGGATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg               <  2:1087263/96‑1 (MQ=255)
                                                                                      ttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa            >  1:4498798/1‑99 (MQ=255)
                                                                                      ttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa            <  1:2068831/99‑1 (MQ=255)
                                                                                        caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa          >  1:655200/4‑99 (MQ=37)
                                                                                        caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa          >  2:683869/4‑99 (MQ=37)
                                                                                          ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc       >  2:3856293/2‑100 (MQ=38)
                                                                                          ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc       >  1:956099/2‑100 (MQ=38)
                                                                                          ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc       >  1:1526381/2‑100 (MQ=38)
                                                                                          ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATCCAGCACTCACTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGCTCGAAAAATAc       >  2:1896040/2‑100 (MQ=25)
                                                                                           gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc      >  2:2254110/1‑100 (MQ=38)
                                                                                           gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc      >  1:2667867/1‑100 (MQ=38)
                                                                                           gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc      >  1:3077202/1‑100 (MQ=255)
                                                                                            tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt     >  1:544844/1‑100 (MQ=38)
                                                                                             cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCATAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg    >  1:5141943/1‑100 (MQ=255)
                                                                                             cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg    >  1:3183533/1‑100 (MQ=255)
                                                                                             cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg    >  1:3914409/1‑100 (MQ=255)
                                                                                                gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGGACGAGCCGGATc                <  2:4196976/85‑1 (MQ=255)
                                                                                                gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc      <  2:2829142/95‑1 (MQ=255)
                                                                                                gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa  <  2:4696059/99‑1 (MQ=255)
                                                                                                 |                                                                                                  
CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATA‑TTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCA  >  NC_000913_3_bme_pgi/1950196‑1950389

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