Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,094 A→G 16.8% I109V (ATC→GTC)  hyaD → hydrogenase 1 maturation protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,0940AG16.8% 368.8 / 65.7 168I109V (ATC→GTC) hyaDhydrogenase 1 maturation protease
Reads supporting (aligned to +/- strand):  ref base A (77/61);  new base G (10/18);  total (87/81)
Fisher's exact test for biased strand distribution p-value = 6.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

Only 100 of 174 total aligned reads displayed.

CTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAC  >  NC_000913/1035995‑1036193
                                                                                                   |                                                                                                   
cTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGata                                                                                                     <  1:1040488/100‑1 (MQ=255)
   cGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATcc                                                                                                  >  2:2643145/1‑100 (MQ=255)
       ccTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCcgcgn                                                                                              <  2:1098745/100‑2 (MQ=255)
       ccTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCgg                                                                                              >  1:612524/1‑100 (MQ=255)
        cTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGa                                                                                             <  1:757712/100‑1 (MQ=255)
         tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGAc                                                                                            >  2:636956/1‑100 (MQ=255)
         tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGAc                                                                                            <  1:1418187/100‑1 (MQ=255)
           tGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACAt                                                                                          >  1:1013736/1‑100 (MQ=255)
           tGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGGGTTGGCGCGCGCGGATATCCGCGGAc                                                                                            >  1:2696818/1‑98 (MQ=255)
            gCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATc                                                                                         >  1:3411975/1‑100 (MQ=255)
             ccGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCt                                                                                        >  2:3305874/1‑100 (MQ=255)
              cGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                       <  2:2549884/100‑1 (MQ=255)
                ttcgaCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGcc                                                                                     <  2:3309692/96‑1 (MQ=38)
                  gaACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGGTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAg                                                                                   <  1:1574920/100‑1 (MQ=255)
                   aaCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGc                                                                                  <  2:3269500/100‑1 (MQ=255)
                    aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGca                                                                                 >  2:3370141/1‑100 (MQ=255)
                    aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                 <  2:4294943/100‑2 (MQ=38)
                    aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                 <  2:1828099/100‑2 (MQ=38)
                    aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                 <  2:2768268/100‑2 (MQ=38)
                     cGCATTCCGGCTTATCTCAGCGCGACGAAAATGAGCCTGCATCATAACAGTGTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGc                                                                                  <  1:536499/98‑1 (MQ=37)
                     cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcac                                                                                <  2:3549550/100‑1 (MQ=255)
                     cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcac                                                                                >  2:2195191/1‑100 (MQ=255)
                     cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                 <  1:1450707/99‑2 (MQ=38)
                       cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgcca                                                                              <  2:2940768/100‑5 (MQ=37)
                       cATTCCGGCGTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                              >  1:3223630/1‑100 (MQ=255)
                        aTTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACata                                                                             >  1:3944249/1‑100 (MQ=255)
                        aTTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                              >  1:3355602/1‑99 (MQ=255)
                         ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgccaac                                                                            <  1:379437/100‑7 (MQ=37)
                          tCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                           >  2:2469552/1‑100 (MQ=255)
                            cGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAACTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                         <  1:2443314/100‑1 (MQ=255)
                             ggCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGcc                                                                        <  2:3262607/100‑1 (MQ=255)
                              gCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGccc                                                                       <  1:2798479/100‑1 (MQ=255)
                                 tATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcg                                                                    <  2:3341916/100‑1 (MQ=255)
                                  aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgt                                                                   <  2:1138688/100‑1 (MQ=255)
                                   tctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtc                                                                  <  2:3835184/100‑1 (MQ=255)
                                   tctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtc                                                                  >  1:3979417/1‑100 (MQ=255)
                                    ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                 >  1:2224347/1‑100 (MQ=255)
                                    ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                 >  1:2291464/1‑100 (MQ=255)
                                    ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                 >  2:3667261/1‑100 (MQ=255)
                                      cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGt                                                               >  2:626672/1‑100 (MQ=255)
                                      cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGt                                                               >  2:745226/1‑100 (MQ=255)
                                       aGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                              <  1:3022896/100‑1 (MQ=255)
                                       aGCGCGAAGAAAATGAGCCTGCATCACAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                              >  2:2680862/1‑100 (MQ=255)
                                         cgcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTg                                                            <  2:1219066/100‑1 (MQ=255)
                                         cgcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTg                                                            <  1:149834/100‑1 (MQ=255)
                                          gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                           >  1:951366/1‑100 (MQ=255)
                                          gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                           >  1:189521/1‑100 (MQ=255)
                                             aagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAc                                                        >  2:4283632/1‑100 (MQ=255)
                                               gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAccc                                                      <  2:1636157/100‑1 (MQ=255)
                                               gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAccc                                                      <  1:3574130/100‑1 (MQ=255)
                                                 aaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGc                                                    <  2:668922/100‑1 (MQ=255)
                                                 aaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGc                                                    <  2:2317882/100‑1 (MQ=255)
                                                  aaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCa                                                   >  1:529362/1‑100 (MQ=255)
                                                         cTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                            <  1:2594342/100‑1 (MQ=255)
                                                         cTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                            <  1:3549324/100‑1 (MQ=255)
                                                          tGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcga                                           <  2:854823/100‑1 (MQ=255)
                                                          tGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcga                                           >  2:1939520/1‑100 (MQ=255)
                                                           gCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                          <  2:3685682/100‑1 (MQ=255)
                                                              tCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCgacgac                                       >  2:1905570/1‑100 (MQ=255)
                                                                 gAACCGTTTCCGCGAAGTGTTCGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATCGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACtcc                                    <  1:2806703/100‑3 (MQ=255)
                                                                  aaCAGTCGCTCCGAAGCGTGGGCGCTGGCCGATCTCGGCGGACTGCCGCCAGATCATATTGCCCTCGTCGGACTGGAGCTCGCAATGCTcgtacagta                                     <  1:3403340/98‑8 (MQ=255)
                                                                   aCAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcgg                                  <  1:4270883/100‑1 (MQ=255)
                                                                    cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggc                                 <  1:2390000/100‑1 (MQ=255)
                                                                      gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                               <  1:2569190/100‑1 (MQ=255)
                                                                        ttCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTa                             >  1:1495812/1‑100 (MQ=255)
                                                                          ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGc                           >  2:504338/1‑100 (MQ=255)
                                                                          ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGc                           >  1:2942150/1‑100 (MQ=255)
                                                                          ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGACTGCAACCCGCAATGCTCGACGACTAcggcgg                               >  1:3688833/1‑96 (MQ=255)
                                                                            ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCt                         >  2:4071112/1‑100 (MQ=255)
                                                                             cGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                        <  2:1005090/100‑1 (MQ=255)
                                                                             cGAAGTGGTGGCGCTGACGAATATCAAAAGACAACTGCCAGCACATATTGACCTCGTAGATCTGCGACCAGCAAAGCCCGACGACTCCGGCGGTAGACt                         >  1:2313280/1‑99 (MQ=255)
                                                                                gcaatatgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                     >  2:2526760/10‑100 (MQ=21)
                                                                                gcaatatgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                     >  2:828845/10‑100 (MQ=21)
                                                                                aGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                     <  1:3496763/100‑1 (MQ=255)
                                                                                aGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                     <  2:1497116/100‑1 (MQ=255)
                                                                                 gtgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCg                    >  2:593814/1‑100 (MQ=255)
                                                                                  tgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGa                   >  1:1183300/1‑100 (MQ=255)
                                                                                  tgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGa                   <  2:100581/100‑1 (MQ=255)
                                                                                   gtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGaa                  >  1:2840313/1‑100 (MQ=255)
                                                                                   gtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGaa                  >  1:3683096/1‑100 (MQ=255)
                                                                                     tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                 <  1:530350/99‑1 (MQ=255)
                                                                                     tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                <  2:973511/100‑1 (MQ=255)
                                                                                      tgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTg               >  1:585750/4‑100 (MQ=37)
                                                                                      ggCGCTGGCGGATCTCCGCGGACACCTGCCAGCACATATCGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTCCGGCGGTAGCCTGAGCGAACTg               <  1:1967258/100‑1 (MQ=255)
                                                                                       gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg              >  2:4138016/3‑100 (MQ=37)
                                                                                       gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg              >  1:2909057/3‑100 (MQ=37)
                                                                                       gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg              >  2:831174/3‑100 (MQ=37)
                                                                                       gcgcTGGCCGATATACGCGGACATCTGCCAGGACCTAGTGCCATCGTCGGTCTGCAACCCGCAATGATCGACGACTACGGCGGTAGCCTGAGCGAACTgg              <  1:1540430/100‑1 (MQ=255)
                                                                                        tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc             >  1:1744388/2‑100 (MQ=38)
                                                                                        tgcTGGCAGATGTCCGCGGATATCCGCCAGAACGTATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGGGCGAACTGGc             >  1:2891110/2‑100 (MQ=25)
                                                                                        cgcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc             >  2:3028402/1‑100 (MQ=255)
                                                                                         gcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCa            <  1:2110672/100‑1 (MQ=255)
                                                                                         gcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCa            >  2:2008348/1‑100 (MQ=255)
                                                                                             gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACggg        <  2:215556/100‑1 (MQ=255)
                                                                                              cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGa       >  1:3463656/1‑100 (MQ=255)
                                                                                              cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGa       >  1:3906877/1‑100 (MQ=255)
                                                                                               ggATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg      <  1:3480168/100‑1 (MQ=255)
                                                                                               ggATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg      >  2:1179209/1‑100 (MQ=255)
                                                                                               ggATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg      <  1:116796/100‑1 (MQ=255)
                                                                                                   atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAc  <  1:2557552/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAC  >  NC_000913/1035995‑1036193

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: