Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,950,3060CG20.1% 419.1 / 106.4 181G73R (GGC→CGC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/81);  new base G (22/14);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 6.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913_3_bme_pgi/1950207‑1950405
                                                                                                    |                                                                                                   
cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                     >  1:726402/1‑100 (MQ=255)
 tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                     >  2:4852527/1‑99 (MQ=255)
  gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    >  2:2635772/1‑99 (MQ=255)
   cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    <  2:1844894/98‑1 (MQ=255)
   cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    <  1:1863547/98‑1 (MQ=255)
   cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                     <  2:1502526/97‑1 (MQ=255)
   cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                   <  2:2662691/99‑1 (MQ=255)
    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                    >  2:3745439/1‑97 (MQ=255)
     cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                 <  2:73124/99‑1 (MQ=255)
     cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                <  1:1811756/100‑1 (MQ=255)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:1012636/100‑1 (MQ=38)
      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:1501883/100‑1 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              <  1:1248725/100‑1 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                              <  2:3230919/100‑1 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                               <  2:4508574/99‑1 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGATGCCCGTGAcc                                                                                              >  2:3406263/1‑100 (MQ=255)
       gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTCAATATTTTTTTCGTCACGGTCACCGTAC‑GTGCGCGTGAc                                                                                               <  1:4892189/99‑1 (MQ=255)
        ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                                             <  1:3295122/100‑3 (MQ=37)
          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                            >  1:1728073/1‑99 (MQ=255)
          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                                            <  1:4062898/99‑4 (MQ=37)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                           <  1:374381/99‑5 (MQ=37)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                           <  1:374338/99‑5 (MQ=37)
           atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                           <  1:1547283/99‑5 (MQ=37)
            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                                         >  1:2133471/1‑100 (MQ=255)
            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                                         >  1:346290/1‑100 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                     >  1:4604269/1‑99 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                     <  2:1719825/99‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                     <  1:1654758/99‑1 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                     >  2:2293941/1‑99 (MQ=255)
                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                    <  1:2572464/100‑1 (MQ=255)
                      cGAGGGCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                           <  2:4297280/88‑5 (MQ=21)
                      cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                               >  2:4764104/1‑100 (MQ=255)
                      cGAGGACGCCAGCACTTCGAATTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCCCGCGCC‑CGTACGGTGCCCGTGACCGGAATGCAg                                                                                     >  2:249193/1‑94 (MQ=255)
                       gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                              >  1:520288/1‑100 (MQ=255)
                       gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATACTTTTTTCTTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                              >  2:3122754/1‑100 (MQ=255)
                        aGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                                 >  2:3519413/1‑96 (MQ=255)
                             gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                                        <  1:3990142/100‑1 (MQ=255)
                              cgaGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       <  2:2121914/98‑1 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                     >  2:4143134/1‑70 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                            <  2:1020673/95‑1 (MQ=255)
                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       >  2:1606319/1‑100 (MQ=255)
                                  cACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCACAACTCATTACGCAGtt                                                                     >  1:4095995/1‑98 (MQ=255)
                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCa                                                                   >  1:92403/1‑98 (MQ=255)
                                    cTTCGATTTCGCCGCTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACTGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                       <  2:2206365/94‑1 (MQ=255)
                                       cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                               <  1:3567126/99‑1 (MQ=255)
                                            tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                          <  1:864868/99‑1 (MQ=255)
                                            tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACGCATTACGCAGTTCAGAGGCCAg                                                           <  2:3867813/98‑1 (MQ=255)
                                            tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAAGCAGAACTCATTACGCAGTTCAGAGGCCAGCt                                                         >  2:3623682/1‑100 (MQ=255)
                                             cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCTCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                        >  1:3333753/1‑100 (MQ=255)
                                              gCCGGTCGCCATATCGCGGTTAATATTTTTGTCGGCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                 <  2:1439381/90‑1 (MQ=255)
                                               ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                       >  1:439497/1‑99 (MQ=255)
                                               ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                       >  1:1280381/1‑99 (MQ=255)
                                                 ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                     >  1:4023620/1‑99 (MQ=255)
                                                  gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   <  1:5000641/100‑1 (MQ=255)
                                                  gTCGCCATATCGCGGTTAATATTTTTTCCGTCACGCGCC‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                   >  1:4052050/1‑100 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTTTTTTTTTTTC‑CGTACGGAGCCCGTGACCCGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTCACGCGt                                                >  2:1145143/1‑100 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                 <  1:1004543/99‑1 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                >  1:4177789/1‑100 (MQ=255)
                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGAGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                >  1:5276951/1‑100 (MQ=255)
                                                      ccATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                >  2:4607545/1‑99 (MQ=255)
                                                       cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTcc                                              >  1:3884936/1‑100 (MQ=255)
                                                           tCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                          <  2:414275/100‑1 (MQ=255)
                                                              cgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                       >  1:5103622/1‑100 (MQ=255)
                                                                 ttAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCg                                    >  2:533338/1‑100 (MQ=255)
                                                                   aaTATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                                  >  2:832057/1‑100 (MQ=255)
                                                                    atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                 <  1:1257864/100‑1 (MQ=255)
                                                                       tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               <  1:2903591/99‑1 (MQ=255)
                                                                       tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               <  1:3205634/99‑1 (MQ=255)
                                                                       tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                               >  1:4663300/1‑99 (MQ=255)
                                                                          ttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                            >  1:4498798/1‑99 (MQ=255)
                                                                          ttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                            <  1:2068831/99‑1 (MQ=255)
                                                                           tttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                          >  1:1471436/1‑100 (MQ=255)
                                                                            caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                          >  1:655200/4‑99 (MQ=37)
                                                                            caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                          >  1:3846778/4‑99 (MQ=37)
                                                                            caggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                          >  2:683869/4‑99 (MQ=37)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  1:1307917/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  2:3315559/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  1:2603770/2‑100 (MQ=38)
                                                                              ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                       >  1:1526381/2‑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)
                                                                                tCACGCGTA‑CGTACGGTGCACGAGACCTGAATGCAGAACTCATTACGCAGTTCAGAGACCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      >  2:1775108/1‑99 (MQ=255)
                                                                                 cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                    >  1:3914409/1‑100 (MQ=255)
                                                                                  aCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc                   <  1:2402610/100‑1 (MQ=255)
                                                                                    gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                      <  2:2829142/95‑1 (MQ=255)
                                                                                    gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa                  <  2:4696059/99‑1 (MQ=255)
                                                                                      gcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa               >  2:180135/1‑100 (MQ=255)
                                                                                       cA‑CGTACGGTGCCCGTGACCTGGATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTCAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                 <  2:5022990/97‑1 (MQ=255)
                                                                                       cA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt              <  2:4521337/100‑1 (MQ=255)
                                                                                        a‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa             <  1:2330619/100‑1 (MQ=255)
                                                                                        a‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa             <  1:4283137/100‑1 (MQ=255)
                                                                                          cGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc            <  1:4540973/100‑1 (MQ=255)
                                                                                           gTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc           <  1:2173240/100‑1 (MQ=255)
                                                                                           gTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc            <  2:2772691/99‑1 (MQ=255)
                                                                                            tACGGTGCCCGTGACCTGACTGCAGAACTCATTGCGCAGTTCAGAGGCCAGc                                                          <  2:1129520/52‑1 (MQ=255)
                                                                                                gTGCCCGTGACCTGAATGCTGAACTCATTACGCTATTCAGAGGCCAGCTTTAACGAGTCCGCACTATCCGGATCGAAAAATACCTGCATGATACCTgcgc      <  2:4422706/100‑1 (MQ=255)
                                                                                                gTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgc      <  2:3078381/100‑1 (MQ=255)
                                                                                                   cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg    <  1:3778773/99‑1 (MQ=255)
                                                                                                   cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGt   <  2:3989355/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTC  >  NC_000913_3_bme_pgi/1950207‑1950405

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: ‑