Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 234,308 T→C 22.1% T159A (ACA→GCA)  gloB ← hydroxyacylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913234,3080TC22.1% 298.3 / 85.8 164T159A (ACA→GCA) gloBhydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base T (65/62);  new base C (16/20);  total (82/82)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 164 total aligned reads displayed.

TTACGATAATAATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAGTGTGT  >  NC_000913/234211‑234403
                                                                                                 |                                                                                               
 tACGATAATAATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATc                                                                                              <  1:3176897/100‑1 (MQ=255)
  aCGATAATAATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTAt                                                                                               >  1:1490255/1‑98 (MQ=255)
    gATAATAATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTc                                                                                           <  1:3062767/100‑1 (MQ=255)
     ataataATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCa                                                                                          <  2:2990929/100‑1 (MQ=255)
     ataataATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCa                                                                                          <  1:1958620/100‑1 (MQ=255)
     ataataAGCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATGTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTc                                                                                           <  2:3191067/99‑1 (MQ=255)
      taataaTCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAg                                                                                         >  2:3749005/1‑100 (MQ=255)
        ataatCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGt                                                                                       <  1:791641/100‑1 (MQ=255)
         taatCATTTATGGACAAATCGTGCGGAAGATTACTCAAAGCAAAATTCATATTTTATAAGCTATATTTATGAGCACAACATACCAATGTTGGGTCAgg                                                                                        >  2:1179246/1‑98 (MQ=255)
         taatCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTa                                                                                      >  1:1114303/1‑100 (MQ=255)
         taatCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTa                                                                                      >  1:544569/1‑100 (MQ=255)
          aatCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTaa                                                                                     >  1:1409313/1‑100 (MQ=255)
          aatCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTaa                                                                                     >  2:2592993/1‑100 (MQ=255)
           atCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAACACAACATACCAATGTAACGTCAGGTAAc                                                                                    >  2:1863608/1‑100 (MQ=255)
              ntttATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGc                                                                                  >  2:249131/2‑99 (MQ=255)
               tttATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCAc                                                                                >  2:535458/1‑100 (MQ=255)
               tttATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCAc                                                                                <  1:1571768/100‑1 (MQ=255)
               tttATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCAc                                                                                <  2:1254641/100‑1 (MQ=255)
                ttATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACt                                                                               <  1:3622434/100‑1 (MQ=255)
                  aTGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAACAAAACATACCAATGTATCGTCAGGTAACGCACTTa                                                                             >  2:1575602/1‑100 (MQ=255)
                   tGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTaa                                                                            >  1:3894555/1‑100 (MQ=255)
                     gACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACt                                                                          >  1:3262455/1‑100 (MQ=255)
                        aaaTCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACtttt                                                                       <  1:1636103/100‑1 (MQ=255)
                         aaTCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACttttt                                                                      >  1:3704411/1‑100 (MQ=255)
                          aTCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTa                                                                     >  1:3565453/1‑100 (MQ=255)
                          aTCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTa                                                                     <  2:1283956/100‑1 (MQ=255)
                           tCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCAAAACATACCAATGTATCGTCAGGTAACGCAATTAACTTTTTaa                                                                    >  1:368780/1‑100 (MQ=255)
                                 ggAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgatt                                                              >  2:2188453/1‑100 (MQ=255)
                                 ggAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGtgatt                                                              >  1:2091718/1‑100 (MQ=255)
                                 ggAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAACACAACATACCAATGTATCGTCAGTTAATGAACTTTACTTTTTAAGtgc                                                                >  2:2913573/1‑97 (MQ=255)
                                 ggAAGAATACTCAAAGCAAACTTCATATTTGATAAGATATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgatt                                                              >  1:3587593/1‑100 (MQ=255)
                                  gaagaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgattg                                                             <  2:149943/100‑1 (MQ=255)
                                  gaagaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgattg                                                             >  2:3697721/1‑100 (MQ=255)
                                  gaagaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgatt                                                              >  2:2460412/1‑99 (MQ=255)
                                    agaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgattgat                                                           >  1:3835364/1‑100 (MQ=255)
                                    agaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGtgattgat                                                           >  2:3097815/1‑100 (MQ=255)
                                     gaaTACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATa                                                          <  2:909811/100‑1 (MQ=255)
                                         ncTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACAtt                                                      >  2:249156/2‑100 (MQ=255)
                                          cTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACAttt                                                     >  1:3356699/1‑100 (MQ=255)
                                          cTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACAttt                                                     >  1:3089944/1‑100 (MQ=255)
                                          cTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGAACTTAACTTTTTAAGTGATTGATACAttt                                                     >  2:365045/1‑100 (MQ=255)
                                           tCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTtg                                                    >  1:1381978/1‑100 (MQ=255)
                                            cAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTtgt                                                   <  1:3204659/100‑1 (MQ=255)
                                            cAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTtgt                                                   >  2:3708733/1‑100 (MQ=255)
                                              aaGCAAACGTCATATTTGATAAGGTATATTCATGAGCACAACATTCCTATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGa                                                 <  1:67250/100‑1 (MQ=255)
                                                 cAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGc                                              >  2:831419/1‑100 (MQ=255)
                                                 cAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGc                                              <  1:1129830/100‑1 (MQ=255)
                                                 cAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGc                                              <  2:1661668/100‑1 (MQ=255)
                                                     cTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGt                                           >  1:720616/1‑99 (MQ=255)
                                                     cTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAATGATTGATACATTTGTGATGCTGt                                           >  2:652356/1‑99 (MQ=255)
                                                     cTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGAAACGCAATTAACTTATTAAGCAATTGATACATTTGTGATGCTGTc                                          >  1:170175/1‑100 (MQ=255)
                                                     cTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTc                                          <  2:1043923/100‑1 (MQ=255)
                                                      ttCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTcc                                         >  2:2433730/1‑100 (MQ=255)
                                                       tCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTccc                                        <  1:1771340/100‑1 (MQ=255)
                                                        cATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCt                                       >  1:97346/1‑100 (MQ=255)
                                                        cATATTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCt                                       <  1:2971034/100‑1 (MQ=255)
                                                          tatTTGATATGGTATATTCATGAACACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCtt                                      >  2:2775425/1‑99 (MQ=255)
                                                          tatTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTc                                     <  1:411057/100‑1 (MQ=255)
                                                           atTTGATAAGGTATATTCATGAGCACAACATATCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCa                                    <  1:2459975/100‑1 (MQ=255)
                                                           atTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCa                                    >  1:190481/1‑100 (MQ=255)
                                                           atTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCa                                    <  2:3456443/100‑1 (MQ=255)
                                                           atTTGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCa                                    <  2:2924993/100‑1 (MQ=255)
                                                             ttGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCaaa                                  <  1:880454/100‑1 (MQ=255)
                                                             ttGATAAGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCaaa                                  <  2:2440751/100‑1 (MQ=255)
                                                               gATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAaaca                                <  1:3820145/100‑1 (MQ=255)
                                                                  aaGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCTCTTATCTTTTTAAGTGTTTGATATATTTGTGATGCTGTGCCTTCAaacaac                              <  2:527513/99‑1 (MQ=255)
                                                                  aaGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAAcc                             <  2:2355482/100‑1 (MQ=255)
                                                                  aaGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAAcc                             <  2:1503503/100‑1 (MQ=255)
                                                                   aGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCg                            >  1:2541221/1‑100 (MQ=255)
                                                                   aGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCg                            >  1:3096655/1‑100 (MQ=255)
                                                                   aGGTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCg                            >  2:1305154/1‑100 (MQ=255)
                                                                    ggTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGa                           >  2:2108277/1‑100 (MQ=255)
                                                                    ggTATATTCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGa                           >  1:2384031/1‑100 (MQ=255)
                                                                     gTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAc                          <  1:3817030/100‑1 (MQ=255)
                                                                     gTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAc                          >  1:1038250/1‑100 (MQ=255)
                                                                     cctatctCATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAc                          <  1:3121988/94‑1 (MQ=255)
                                                                      tataTTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAcc                         <  1:3546071/100‑1 (MQ=255)
                                                                      tataTTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAcc                         <  1:246589/100‑1 (MQ=255)
                                                                       atatTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACca                        <  1:2996222/100‑1 (MQ=255)
                                                                          ttCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCacac                     >  1:1485087/1‑100 (MQ=255)
                                                                            cATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCAcaccc                   >  1:1245680/1‑100 (MQ=255)
                                                                            cATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCAcaccc                   >  1:3911624/1‑100 (MQ=255)
                                                                            cATGAGCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGAcagatcgg                   >  2:3121958/1‑93 (MQ=255)
                                                                             aTGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCAcaccca                  <  2:3887186/100‑1 (MQ=255)
                                                                                aGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACccacca               <  1:626793/100‑1 (MQ=255)
                                                                                 gCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAg              >  1:1646583/1‑100 (MQ=255)
                                                                                 gCACAACATACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAg              >  1:2137167/1‑100 (MQ=255)
                                                                                  cacaACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGa             <  2:745578/100‑1 (MQ=255)
                                                                                   acaacaTACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGa             <  2:1037667/99‑1 (MQ=255)
                                                                                   acaacaTACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGaa            >  1:3359499/1‑100 (MQ=255)
                                                                                      actgaacAATGTATCGTCAGGTAACGCAATTAACTTGTTAAGTTATTGATACATTTGTGATGTTGTCCCTTCAATCTACCGACCACAACCACCAGAAAAc         <  2:2661216/94‑1 (MQ=255)
                                                                                         tACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgt      >  2:3818904/1‑100 (MQ=255)
                                                                                         tACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgt      >  1:1213142/1‑100 (MQ=255)
                                                                                         tACCAATGCATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAg       <  1:602678/99‑1 (MQ=255)
                                                                                           ccAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgt    >  1:2782796/1‑100 (MQ=255)
                                                                                           ccAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgt    >  1:2753479/1‑100 (MQ=255)
                                                                                           ccAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgt    >  1:1810588/1‑100 (MQ=255)
                                                                                            cAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgtg   <  1:1043403/100‑1 (MQ=255)
                                                                                            cAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgtg   <  1:2779478/100‑1 (MQ=255)
                                                                                            cAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAgtgtg   <  2:3092927/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TTACGATAATAATCATTTATGGACAAATCGTGCGGAAGAATACTCAAAGCAAACTTCATATTTGATAAGGTATATTCATGAGCACAACATACCAATGTATCGTCAGGTAACGCACTTAACTTTTTAAGTGATTGATACATTTGTGATGCTGTCCCTTCAAACAACCGACCACACCCACCAGAAAACAGTGTGT  >  NC_000913/234211‑234403

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