Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,561:1 +C 17.2% coding (1120/1452 nt) uxaB ← altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,5611.C17.2% 306.6 / 73.1 99coding (1120/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base . (42/40);  new base C (7/10);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 102 total aligned reads displayed.

CCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGC  >  NC_000913/1609464‑1609657
                                                                                                  |                                                                                                
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                                >  1:2998823/1‑100 (MQ=255)
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                                >  2:3963305/1‑100 (MQ=255)
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                                <  2:2579412/100‑1 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAg                                                                                            >  2:1610858/1‑100 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAg                                                                                            <  2:742734/100‑1 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAAGGCGAAAGTACGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAg                                                                                            <  2:2748077/100‑1 (MQ=255)
     aCGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgg                                                                                           <  2:1551283/100‑1 (MQ=255)
      cgcgCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                          <  2:1439957/100‑1 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAgg                                                                                           >  2:1880124/1‑99 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                          <  1:2653402/100‑1 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                          <  2:1301002/100‑2 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                          <  1:2277647/100‑2 (MQ=255)
      cGCTCACCACGAGAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgg                                                                                           <  2:393550/99‑1 (MQ=255)
       gCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggc                                                                                         <  1:2318444/100‑3 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATg                                                                                        >  1:1692133/1‑100 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATg                                                                                        >  1:2351499/1‑100 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggca                                                                                        <  1:3106828/100‑4 (MQ=255)
          caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggcaga                                                                                      <  1:1607051/100‑6 (MQ=38)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                    <  1:300370/100‑1 (MQ=255)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                    >  1:1321742/1‑100 (MQ=255)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATCCggg                                                                                    <  1:2569958/100‑1 (MQ=255)
             cacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggcagaagg                                                                                   <  1:1493493/100‑9 (MQ=38)
              acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGt                                                                                   >  2:1779315/1‑99 (MQ=255)
               cGATAGAACGCAATAAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGTCGTTTGCCTTCTGCACTGCCAGCAGCTG‑TGGCAGGATGCGGGTg                                                                                  >  1:1076375/1‑99 (MQ=255)
                gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCg                                                                                <  2:550555/100‑1 (MQ=255)
                     aaCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACt                                                                           >  2:2312426/1‑100 (MQ=255)
                      aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACtt                                                                          >  1:4049756/1‑100 (MQ=255)
                       cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTg                                                                         >  2:1387996/1‑100 (MQ=255)
                       cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACt                                                                           <  2:3581200/98‑1 (MQ=255)
                        gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                        <  1:2412609/100‑1 (MQ=255)
                        gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                        <  1:3953642/100‑1 (MQ=255)
                         caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                       >  2:2551418/1‑100 (MQ=255)
                         caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                       <  1:3038847/100‑1 (MQ=255)
                           atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                     <  2:181821/100‑1 (MQ=255)
                           atcaatGCCGCTAATGCCAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                       <  2:1947953/98‑1 (MQ=255)
                            tcaatGCCGCTAATGCGAAAGTAAGGAGCGACGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCAt                                                                    >  2:160847/1‑100 (MQ=255)
                             caatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATa                                                                   >  2:2636439/1‑100 (MQ=255)
                              aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                  <  2:3517282/100‑1 (MQ=255)
                              aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                  >  1:3688516/1‑100 (MQ=255)
                                   ngCTAATGCGAAAGTAAGGCGCGCCGGACGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAGCGTTg                                                             >  2:1520437/2‑100 (MQ=255)
                                   cGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                             <  1:725087/100‑1 (MQ=255)
                                      tAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgc                                                          <  2:2337734/100‑1 (MQ=255)
                                         tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                       <  1:3349543/100‑1 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  2:818918/1‑100 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     <  1:1374515/100‑1 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  1:38329/1‑100 (MQ=255)
                                             aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                   <  2:135702/100‑1 (MQ=255)
                                              aaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTTTGCTCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  2:1619498/1‑97 (MQ=255)
                                              aaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGa                                                  >  2:3625595/1‑100 (MQ=255)
                                                gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                >  1:1938204/1‑100 (MQ=255)
                                                gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                <  2:3267247/100‑1 (MQ=255)
                                                  aaGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagc                                              <  1:451299/100‑1 (MQ=255)
                                                   aGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagca                                             >  1:2665775/1‑100 (MQ=255)
                                                    ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                            <  1:344749/100‑1 (MQ=255)
                                                     gcgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagc                                           <  2:3114285/100‑1 (MQ=255)
                                                      cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGCTAGACAGCAGCt                                          >  1:2890968/1‑100 (MQ=255)
                                                        cgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGa                                        >  1:1631799/1‑100 (MQ=255)
                                                         gcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCACCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                       <  2:2784636/100‑1 (MQ=255)
                                                         gcCGGAAGTGTGCCGTTTGCCGTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                       <  2:2784622/100‑1 (MQ=255)
                                                          ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                      <  1:1520449/100‑1 (MQ=255)
                                                             gAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCtt                                   >  2:2684855/1‑100 (MQ=255)
                                                               aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                 <  2:4128171/100‑1 (MQ=255)
                                                               aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                 <  2:38329/100‑1 (MQ=255)
                                                                gtgtGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAAt                                <  1:544869/100‑1 (MQ=255)
                                                                gtgtGCCGGTTGCCTCCTCCCCTGCCAGCAGCTG‑TGTCAGGATGCGGGTGCGGAATTTGGTCATACCGTTGAGCGCGCTAGACAGCAGCTGATGCtt                                   <  2:3356873/97‑1 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:2651432/1‑100 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:1191458/1‑100 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:3094686/1‑100 (MQ=255)
                                                                   tgCCGTTTGCCTTCTGCTCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTa                             >  2:1897652/1‑100 (MQ=255)
                                                                    gCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTAc                            >  1:2569322/1‑100 (MQ=255)
                                                                    gCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTAc                            >  2:108033/1‑100 (MQ=255)
                                                                     ccGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACg                           <  1:3396202/100‑1 (MQ=255)
                                                                       gtgtGCCGTCTGCCCTGCCAGCAGCTC‑TGGCAGGATGCGGGCGCGGAACTTTGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTa                             <  2:3714312/93‑1 (MQ=255)
                                                                        tttGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGt                        <  1:2551432/100‑1 (MQ=255)
                                                                            tctcgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:678257/92‑1 (MQ=38)
                                                                            ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:672556/100‑1 (MQ=255)
                                                                            ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:953231/100‑1 (MQ=255)
                                                                             cTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACgg                   >  1:618156/1‑100 (MQ=255)
                                                                             cTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  2:2514600/99‑1 (MQ=255)
                                                                               cgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACgaga                 >  2:678251/6‑97 (MQ=38)
                                                                               cgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaa                 >  2:1909692/6‑100 (MQ=38)
                                                                               cgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaa                 >  1:262925/6‑100 (MQ=38)
                                                                                 catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc               >  2:1595824/4‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGTAGCTGATGCTTAATGTACGGGTTACGGAAACg              <  1:3922202/100‑1 (MQ=255)
                                                                                  gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAACCg              >  2:3341606/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  2:1222552/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:1200710/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:525856/1‑100 (MQ=255)
                                                                                  atcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATCGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:1061714/3‑100 (MQ=255)
                                                                                   ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc             >  1:2506539/1‑100 (MQ=255)
                                                                                   atcctgcgacagcTGCTGGCAGGCTGGGGGTGAGGCACTTGGTCATACCGTTGAGCGCGATAGACAGCGGCGGATGCTTAATGTACGGGTTAAGGaa                 >  2:816604/10‑97 (MQ=255)
                                                                                    cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:1893485/1‑100 (MQ=255)
                                                                                    cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:2538578/1‑100 (MQ=255)
                                                                                    cctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:3934619/1‑100 (MQ=255)
                                                                                       gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGccc           >  2:493689/1‑98 (MQ=255)
                                                                                       gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg         <  1:2678556/100‑1 (MQ=255)
                                                                                       gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg         <  1:1865324/100‑1 (MQ=255)
                                                                                          agcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGa      >  2:1907172/1‑100 (MQ=255)
                                                                                           gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc     <  2:2616630/100‑1 (MQ=255)
                                                                                              gcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGc  >  1:2233831/1‑100 (MQ=255)
                                                                                                  |                                                                                                
CCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGC  >  NC_000913/1609464‑1609657

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