Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,103 C→T 16.0% H112Y (CAT→TAT)  hyaD → hydrogenase 1 maturation protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,1030CT16.0% 293.7 / 39.9 128H112Y (CAT→TAT) hyaDhydrogenase 1 maturation protease
Reads supporting (aligned to +/- strand):  ref base C (59/48);  new base T (13/8);  total (72/56)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 129 total aligned reads displayed.

TGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTG  >  NC_000913/1036006‑1036202
                                                                                                 |                                                                                                   
tGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACAt                                                                                                   >  2:2146910/1‑100 (MQ=255)
  ctggAGAACGCATTCCGGCTTATCTCAGCGCGAAGAATATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGTGCTGCCGGATATACGTGGACATc                                                                                                  <  2:2697499/97‑1 (MQ=255)
  ccGGAGAACGCATTCCGGCTTATCTCAGCGGGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCt                                                                                                 >  2:78431/1‑100 (MQ=255)
  ccGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCt                                                                                                 >  1:2150444/1‑100 (MQ=255)
   cGGAGAACGCATTCCGGCTTATCTCATCGCGAAGAAAATGAGCTTGCATCAGAACAGTTTCTCCGAAGTGTTTGCGCTGGCGGATATCCGCGGACATCTg                                                                                                <  2:1311631/100‑1 (MQ=255)
   cGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                <  2:2934427/100‑1 (MQ=255)
   cGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                <  1:106842/100‑1 (MQ=255)
   cGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                >  1:2643566/1‑100 (MQ=255)
     gagaACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGcc                                                                                              >  1:2724884/1‑100 (MQ=255)
       gagcGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGcct                                                                                             <  2:1417335/96‑2 (MQ=255)
        aaCGCATTCCGGCTTATCTCTGCGCTAATAATATGATTCTTCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGc                                                                                           <  1:762401/100‑1 (MQ=255)
        aaCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGTATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAg                                                                                            <  1:808277/99‑1 (MQ=255)
        aaCGCATTCCGGCTTATATCATTGCGAAGAAATTGTGCGTGCATCAGAACAGTTTCGCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATc                                                                                                  <  1:142490/93‑1 (MQ=255)
         aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTACCAGca                                                                                          >  1:628479/1‑100 (MQ=255)
         aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                          <  1:2433246/100‑2 (MQ=38)
         aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGACGCTGGCGGATATCCGCGGACATCTGCCAGca                                                                                          >  1:1133784/1‑100 (MQ=255)
          cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcac                                                                                         <  2:1405051/100‑1 (MQ=255)
          cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgc                                                                                         <  2:2216921/100‑3 (MQ=37)
           gCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcaca                                                                                        <  1:1299754/100‑1 (MQ=255)
            cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                                       >  1:1753957/1‑100 (MQ=255)
            cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                                       >  1:2177803/1‑100 (MQ=255)
            cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgcca                                                                                       <  1:1762126/100‑5 (MQ=37)
             aTTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGGAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACata                                                                                      <  1:1079239/100‑1 (MQ=255)
              ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACatat                                                                                     >  2:1537337/1‑100 (MQ=255)
                ccGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                                    >  1:2574912/1‑99 (MQ=255)
                ccGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTg                                                                                   <  1:2910688/100‑1 (MQ=255)
                  ggCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTgt                                                                                  >  2:716780/1‑98 (MQ=255)
                  ggCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGcc                                                                                 <  2:2585984/100‑1 (MQ=255)
                    nttATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTACCCt                                                                               >  2:1359774/2‑100 (MQ=255)
                    cTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGccc                                                                                <  1:2992193/99‑1 (MQ=255)
                       aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgt                                                                            >  2:1191562/1‑100 (MQ=255)
                       aTCTCAGCCCTAATAAAATGAGCTTCCATCATAACAGTTTCTCCTAAGTGTTTGCTCTGGCTTTTATCCGCGGACATCTGCCAGcac                                                                                         <  2:185682/87‑1 (MQ=255)
                        tctcATCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtc                                                                           <  2:2602409/100‑1 (MQ=255)
                        tctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgt                                                                            >  1:2363032/1‑99 (MQ=255)
                         ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                          <  2:1124772/100‑1 (MQ=255)
                         ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGACGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                          >  1:109057/1‑100 (MQ=255)
                            aGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                                       >  2:879061/1‑100 (MQ=255)
                              cgcgAAGAAAATAAGCCTGCATCAGAACAGTTTCGACGAAGTGTTGGCGCTGGCGGAAAACCGCGGACATCTGCCAGCACATATTGCGCgcgtacgt                                                                        >  2:2136095/1‑93 (MQ=255)
                               gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                                       <  2:558907/97‑1 (MQ=255)
                                cgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCa                                                                   >  2:2383770/1‑100 (MQ=255)
                                 gaagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCaa                                                                  >  1:371236/1‑100 (MQ=255)
                                 gaagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCaa                                                                  >  1:1493817/1‑100 (MQ=255)
                                   agaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGTACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCa                                                                   <  1:2142153/97‑1 (MQ=255)
                                    gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAccc                                                               >  2:1186883/1‑100 (MQ=255)
                                       aaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCa                                                            <  1:1901039/100‑1 (MQ=255)
                                         tGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAAt                                                          <  1:1819067/100‑1 (MQ=255)
                                           ngCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATATGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATg                                                         >  2:1325470/2‑99 (MQ=255)
                                            gCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGTATATCCGCGGATATCTGCCAGCACATATTGCCCTCGTTGGTCTGCAACCCGCAATGCt                                                       <  2:477649/100‑1 (MQ=255)
                                            gCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCt                                                       >  1:607959/1‑100 (MQ=255)
                                             ccTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTc                                                      >  2:1676856/1‑100 (MQ=255)
                                             ccTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTc                                                      <  1:2383797/100‑1 (MQ=255)
                                              cTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                                     >  1:2754203/1‑100 (MQ=255)
                                               tGCATCAGAACTGTTTCTCCGAAGTGTTGGCGCTTGCGGATATCCGCGGACATCTGCCAGCACATATTACACTCATCGGTCTGCAACCCGCAATGcacaa                                                    >  2:2338003/1‑96 (MQ=255)
                                               tGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcga                                                    >  2:1701231/1‑100 (MQ=255)
                                                gCATCAGAACAGTTTCTCCTATGTGTTGGCGCTGGTGGATATCCGCGGACATCTGCCAGCACATATGGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                   <  1:538448/100‑1 (MQ=255)
                                                gCATCAGAACAGTTTCTCCGAAGTGTTGGCGTTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                   <  1:406306/100‑1 (MQ=255)
                                                gCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                   <  2:2099303/100‑1 (MQ=255)
                                                 cATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATCTCCGCGCACATCTCCCAGAACATATTACCCTCGTCAAACTACAACCCGCAATGCTcgacg                                                  >  1:1109401/1‑100 (MQ=255)
                                                   tCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgacg                                                  <  1:1185966/98‑1 (MQ=255)
                                                         cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggc                                          <  1:1512615/100‑1 (MQ=255)
                                                         cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggc                                          >  1:2239321/1‑100 (MQ=255)
                                                           gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                        <  2:2302920/100‑1 (MQ=255)
                                                           gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                        >  1:1525869/1‑100 (MQ=255)
                                                           gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                        <  1:1325482/100‑1 (MQ=255)
                                                           gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcag                                        >  2:1838415/1‑98 (MQ=255)
                                                            tttCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACAACTACGGCGGt                                       >  1:275792/1‑100 (MQ=255)
                                                             ttCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTa                                      >  2:928676/1‑100 (MQ=255)
                                                             ttCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAATGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTa                                      >  1:2508693/1‑100 (MQ=255)
                                                                 ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCACTGCTCGACGACTACGGCGGTAGCCt                                  >  2:2690042/1‑100 (MQ=255)
                                                                 ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGcc                                   >  1:647551/1‑99 (MQ=255)
                                                                 ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGATCTACAACCCGCAATGCTCGACGACTACGGCGGTAGCCt                                  >  1:2346383/1‑100 (MQ=255)
                                                                  cGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                 <  2:1570797/100‑1 (MQ=255)
                                                                  cGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                 <  2:2643534/100‑1 (MQ=255)
                                                                   gAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                 >  2:626451/1‑99 (MQ=255)
                                                                     aGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGTTCGACGACTACGGCGGTAGCCTGAGc                              <  1:1849240/100‑1 (MQ=255)
                                                                       aatatgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGa                            >  2:1694902/8‑100 (MQ=25)
                                                                         tatgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGaa                           >  2:1928994/6‑99 (MQ=37)
                                                                          atgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                         >  2:2894216/5‑100 (MQ=37)
                                                                          atgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                         >  1:2931226/5‑100 (MQ=37)
                                                                          atgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCACAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                         >  1:504373/5‑100 (MQ=25)
                                                                           ggCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCACAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTg                        >  1:320379/1‑100 (MQ=255)
                                                                            gtgcTGGGAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg                       >  1:1581671/3‑100 (MQ=37)
                                                                            gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGATAGACGACTACGGCGGTAGCCTGAGCGAACTgg                       >  1:2381660/3‑100 (MQ=25)
                                                                             tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                         >  1:889218/2‑97 (MQ=37)
                                                                             tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc                      >  2:620757/2‑100 (MQ=38)
                                                                                tGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACg                   >  1:219901/1‑100 (MQ=255)
                                                                                 ggCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACgg                  >  1:2821524/1‑100 (MQ=255)
                                                                                 ggCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACgg                  >  2:2019333/1‑100 (MQ=38)
                                                                                   cGGATTTCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGTTCTGCAACCCGCACTGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACgg                  <  2:623093/98‑1 (MQ=255)
                                                                                    agatgtCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg               <  1:2019356/95‑1 (MQ=255)
                                                                                        atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAc           >  1:282329/1‑100 (MQ=255)
                                                                                          ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACt          >  2:757911/1‑99 (MQ=255)
                                                                                          ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg         >  2:2673382/1‑100 (MQ=255)
                                                                                          ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg         >  1:2397755/1‑100 (MQ=255)
                                                                                           cgcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGc        <  1:84833/100‑1 (MQ=255)
                                                                                             cgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGccc      <  2:2724852/100‑1 (MQ=255)
                                                                                             cgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGccc      >  1:1315841/1‑100 (MQ=255)
                                                                                             cgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGcc       >  1:705246/1‑99 (MQ=255)
                                                                                              ggACATCTGCCAGCACATATTGCCCTTGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGccc      <  1:1163041/99‑1 (MQ=255)
                                                                                               gACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGc    <  2:2884462/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTG  >  NC_000913/1036006‑1036202

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