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.

CCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGC  >  NC_000913/1609464‑1609657
                                                                                                  |                                                                                                
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgtg                                                                                                >  1:2998823/1‑100 (MQ=255)
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgtg                                                                                                >  2:3963305/1‑100 (MQ=255)
ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgtg                                                                                                <  2:2579412/100‑1 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAg                                                                                            >  2:1610858/1‑100 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAg                                                                                            <  2:742734/100‑1 (MQ=255)
    tACGCTCACCACGATAGAACGCAATCAATGCCGCTAAGGCGAAAGTACGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAg                                                                                            <  2:2748077/100‑1 (MQ=255)
     aCGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgg                                                                                           <  2:1551283/100‑1 (MQ=255)
      cgcgCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGa                                                                                          <  2:1439957/100‑1 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAgg                                                                                           >  2:1880124/1‑99 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGa                                                                                          <  1:2653402/100‑1 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAggg                                                                                          <  2:1301002/100‑2 (MQ=255)
      cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAggg                                                                                          <  1:2277647/100‑2 (MQ=255)
      cGCTCACCACGAGAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgg                                                                                           <  2:393550/99‑1 (MQ=255)
       gCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgggc                                                                                         <  1:2318444/100‑3 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATg                                                                                        >  1:1692133/1‑100 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATg                                                                                        >  1:2351499/1‑100 (MQ=255)
        cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgggca                                                                                        <  1:3106828/100‑4 (MQ=255)
          caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgggcaga                                                                                      <  1:1607051/100‑6 (MQ=38)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                    <  1:300370/100‑1 (MQ=255)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                    >  1:1321742/1‑100 (MQ=255)
            ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATCCggg                                                                                    <  1:2569958/100‑1 (MQ=255)
             cacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAgggcagaagg                                                                                   <  1:1493493/100‑9 (MQ=38)
              acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGt                                                                                   >  2:1779315/1‑99 (MQ=255)
               cGATAGAACGCAATAAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGTCGTTTGCCTTCTGCACTGCCAGCAGCTGTGGCAGGATGCGGGTg                                                                                  >  1:1076375/1‑99 (MQ=255)
                gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCg                                                                                <  2:550555/100‑1 (MQ=255)
                     aaCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACt                                                                           >  2:2312426/1‑100 (MQ=255)
                      aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACtt                                                                          >  1:4049756/1‑100 (MQ=255)
                       cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                         >  2:1387996/1‑100 (MQ=255)
                       cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACt                                                                           <  2:3581200/98‑1 (MQ=255)
                        gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                        <  1:2412609/100‑1 (MQ=255)
                        gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                        <  1:3953642/100‑1 (MQ=255)
                         caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGt                                                                       >  2:2551418/1‑100 (MQ=255)
                         caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGt                                                                       <  1:3038847/100‑1 (MQ=255)
                           atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                     <  2:181821/100‑1 (MQ=255)
                           atcaatGCCGCTAATGCCAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGt                                                                       <  2:1947953/98‑1 (MQ=255)
                            tcaatGCCGCTAATGCGAAAGTAAGGAGCGACGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCAt                                                                    >  2:160847/1‑100 (MQ=255)
                             caatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATa                                                                   >  2:2636439/1‑100 (MQ=255)
                              aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                  <  2:3517282/100‑1 (MQ=255)
                              aatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATAc                                                                  >  1:3688516/1‑100 (MQ=255)
                                   ngCTAATGCGAAAGTAAGGCGCGCCGGACGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATAGCGTTg                                                             >  2:1520437/2‑100 (MQ=255)
                                   cGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                             <  1:725087/100‑1 (MQ=255)
                                      tAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgc                                                          <  2:2337734/100‑1 (MQ=255)
                                         tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                       <  1:3349543/100‑1 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  2:818918/1‑100 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     <  1:1374515/100‑1 (MQ=255)
                                           cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  1:38329/1‑100 (MQ=255)
                                             aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                   <  2:135702/100‑1 (MQ=255)
                                              aaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTTTGCTCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                     >  2:1619498/1‑97 (MQ=255)
                                              aaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGa                                                  >  2:3625595/1‑100 (MQ=255)
                                                gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                >  1:1938204/1‑100 (MQ=255)
                                                gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                                <  2:3267247/100‑1 (MQ=255)
                                                  aaGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagc                                              <  1:451299/100‑1 (MQ=255)
                                                   aGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagca                                             >  1:2665775/1‑100 (MQ=255)
                                                    ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                            <  1:344749/100‑1 (MQ=255)
                                                     gcgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagc                                           <  2:3114285/100‑1 (MQ=255)
                                                      cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGCTAGACAGCAGCt                                          >  1:2890968/1‑100 (MQ=255)
                                                        cgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGa                                        >  1:1631799/1‑100 (MQ=255)
                                                         gcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCACCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                       <  2:2784636/100‑1 (MQ=255)
                                                         gcCGGAAGTGTGCCGTTTGCCGTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                       <  2:2784622/100‑1 (MQ=255)
                                                          ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                      <  1:1520449/100‑1 (MQ=255)
                                                             gAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCtt                                   >  2:2684855/1‑100 (MQ=255)
                                                               aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                 <  2:4128171/100‑1 (MQ=255)
                                                               aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                                 <  2:38329/100‑1 (MQ=255)
                                                                gtgtGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAAt                                <  1:544869/100‑1 (MQ=255)
                                                                gtgtGCCGGTTGCCTCCTCCCCTGCCAGCAGCTGTGTCAGGATGCGGGTGCGGAATTTGGTCATACCGTTGAGCGCGCTAGACAGCAGCTGATGCtt                                   <  2:3356873/97‑1 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:2651432/1‑100 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:1191458/1‑100 (MQ=255)
                                                                  gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                              >  1:3094686/1‑100 (MQ=255)
                                                                   tgCCGTTTGCCTTCTGCTCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTa                             >  2:1897652/1‑100 (MQ=255)
                                                                    gCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTAc                            >  1:2569322/1‑100 (MQ=255)
                                                                    gCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTAc                            >  2:108033/1‑100 (MQ=255)
                                                                     ccGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACg                           <  1:3396202/100‑1 (MQ=255)
                                                                       gtgtGCCGTCTGCCCTGCCAGCAGCTCTGGCAGGATGCGGGCGCGGAACTTTGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTa                             <  2:3714312/93‑1 (MQ=255)
                                                                        tttGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGt                        <  1:2551432/100‑1 (MQ=255)
                                                                            tctcgcatcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:678257/92‑1 (MQ=38)
                                                                            ccTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:672556/100‑1 (MQ=255)
                                                                            ccTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  1:953231/100‑1 (MQ=255)
                                                                             cTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACgg                   >  1:618156/1‑100 (MQ=255)
                                                                             cTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                    <  2:2514600/99‑1 (MQ=255)
                                                                               cgcatcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACgaga                 >  2:678251/6‑97 (MQ=38)
                                                                               cgcatcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaa                 >  2:1909692/6‑100 (MQ=38)
                                                                               cgcatcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaa                 >  1:262925/6‑100 (MQ=38)
                                                                                 catcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc               >  2:1595824/4‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGTAGCTGATGCTTAATGTACGGGTTACGGAAACg              <  1:3922202/100‑1 (MQ=255)
                                                                                  gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAACCg              >  2:3341606/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  2:1222552/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:1200710/1‑100 (MQ=255)
                                                                                  gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:525856/1‑100 (MQ=255)
                                                                                  atcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATCGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg              >  1:1061714/3‑100 (MQ=255)
                                                                                   ccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc             >  1:2506539/1‑100 (MQ=255)
                                                                                   atcctgcgacagcTGCTGGCAGGCTGGGGGTGAGGCACTTGGTCATACCGTTGAGCGCGATAGACAGCGGCGGATGCTTAATGTACGGGTTAAGGaa                 >  2:816604/10‑97 (MQ=255)
                                                                                    cctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:1893485/1‑100 (MQ=255)
                                                                                    cctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:2538578/1‑100 (MQ=255)
                                                                                    cctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc            >  1:3934619/1‑100 (MQ=255)
                                                                                       gccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGccc           >  2:493689/1‑98 (MQ=255)
                                                                                       gccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg         <  1:2678556/100‑1 (MQ=255)
                                                                                       gccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg         <  1:1865324/100‑1 (MQ=255)
                                                                                          agcagcTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGa      >  2:1907172/1‑100 (MQ=255)
                                                                                           gcagcTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc     <  2:2616630/100‑1 (MQ=255)
                                                                                              gcTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGc  >  1:2233831/1‑100 (MQ=255)
                                                                                                  |                                                                                                
CCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGC  >  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: