Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,555 Δ1 bp 15.5% coding (1126/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,5550G.15.5% 322.2 / 57.3 110coding (1126/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (49/44);  new base . (11/6);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-02

Only 100 of 128 total aligned reads displayed.

TGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  NC_000913/1609457‑1609654
                                                                                                  |                                                                                                    
tGTTTCGCCGTTACGTTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagc                                                                                                     <  1:1633437/100‑1 (MQ=255)
tGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagc                                                                                                     >  1:3402681/1‑100 (MQ=255)
tGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCACGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcag                                                                                                      >  1:239972/1‑99 (MQ=255)
 gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagca                                                                                                    <  1:3438807/100‑1 (MQ=255)
 gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagca                                                                                                    <  1:1681257/100‑1 (MQ=255)
  tttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagcag                                                                                                   >  2:3420949/1‑100 (MQ=255)
   ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagatcg                                                                                                  >  2:2464868/1‑96 (MQ=255)
   ttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCagcagc                                                                                                  >  1:2160197/1‑100 (MQ=255)
      gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑t                                                                                              >  1:3197802/1‑100 (MQ=255)
       ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                             >  1:1456982/1‑100 (MQ=255)
       ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                             <  1:3413010/100‑1 (MQ=255)
         gTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGc                                                                                           <  2:554130/100‑1 (MQ=255)
         gTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGc                                                                                           >  1:1997836/1‑100 (MQ=255)
          ttACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCTCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCa                                                                                          <  1:185443/100‑1 (MQ=255)
           tACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAg                                                                                         <  1:3773570/100‑1 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:1502557/100‑2 (MQ=255)
               cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGATTGTTAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggc                                                                                      <  1:1949148/99‑3 (MQ=255)
               cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATg                                                                                     >  1:2787461/1‑100 (MQ=255)
               cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATg                                                                                     >  1:537892/1‑100 (MQ=255)
               cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggc                                                                                      <  1:3609767/99‑3 (MQ=255)
                 caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCg                                                                                   <  2:242692/100‑1 (MQ=255)
                 caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCg                                                                                   >  1:863359/1‑100 (MQ=255)
                   ccacGATATAACGCAATCAATGCCGCTAATGGGAAAGTAAAGCGCGCAGGAAGTGTGCCGTTTGCCTACTACCCTGCCAGCAGCTT‑TAGCa                                                                                          >  2:1098069/1‑91 (MQ=255)
                   ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 <  2:3140329/100‑1 (MQ=255)
                   ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 >  1:1065120/1‑100 (MQ=255)
                   ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 >  2:24415/1‑100 (MQ=255)
                     acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTg                                                                               >  2:3787789/1‑100 (MQ=255)
                      cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGc                                                                              <  1:1530346/100‑1 (MQ=255)
                        aTAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCgg                                                                            <  1:1522235/100‑1 (MQ=255)
                         tAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGa                                                                           <  2:2689747/100‑1 (MQ=255)
                           gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAAc                                                                         <  1:1519429/100‑1 (MQ=255)
                           gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAAc                                                                         <  2:4014/100‑1 (MQ=255)
                           gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAAc                                                                         >  1:2203171/1‑100 (MQ=255)
                           gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAAc                                                                         >  1:399276/1‑100 (MQ=255)
                             aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACtt                                                                       >  2:1527076/1‑100 (MQ=255)
                               gCAACCAAAGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCGGTTTGCCTTCTGCCCTGCCAGTAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgt                                                                     <  2:2798708/100‑2 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                     <  1:2850750/99‑1 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGACGGAAATGTGCCGTTTAACTTCTGCCCTGACAGAAGCTG‑TGGCAAAATGCGGGTGcaaaac                                                                         >  2:2988125/1‑90 (MQ=255)
                                caagcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    <  2:843908/100‑1 (MQ=255)
                                 aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTc                                                                   <  2:1390252/100‑1 (MQ=255)
                                 aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTc                                                                   <  2:40317/100‑1 (MQ=255)
                                 aatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTc                                                                   >  1:1859161/1‑100 (MQ=255)
                                  atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                  <  1:2867779/100‑1 (MQ=255)
                                    caatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATa                                                                >  1:2075796/1‑100 (MQ=255)
                                    caatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATa                                                                >  1:6801/1‑100 (MQ=255)
                                       tGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCg                                                             >  2:1899231/1‑100 (MQ=255)
                                          cGCTAATGCGAAAGTAAGGCGCACCGGAAGTGTGCCGCTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                          <  2:2648760/100‑1 (MQ=255)
                                              aaTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAATTTGGTCATACCGTTGAgcg                                                      <  2:1937994/100‑1 (MQ=255)
                                              aaTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                      >  1:3215566/1‑100 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:203479/1‑100 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  1:163424/1‑100 (MQ=255)
                                                 gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                   >  2:3657211/1‑100 (MQ=255)
                                                  cGAAAGTAAGGCGCTCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                  <  1:1843321/100‑1 (MQ=255)
                                                  cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGTGCTAt                                                  <  2:275117/100‑1 (MQ=255)
                                                   gAAAGTAAGGCGCGGCGGAAGTGTGCCGTTTGCCTTCTGCTCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                 <  1:323008/100‑1 (MQ=255)
                                                   gAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                 <  2:494795/100‑1 (MQ=255)
                                                    aaaGTAAGGCGCTCCGGAAGTGTGCCGTTTGCCTTTTTCCGTGCCAGCAGGTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                <  1:1297240/100‑1 (MQ=255)
                                                    aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                >  2:335087/1‑100 (MQ=255)
                                                    aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGCAACTTGGTCATACCGTTGAGCGCGATAg                                                <  2:3730104/100‑1 (MQ=255)
                                                      aGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAc                                              >  2:776148/1‑100 (MQ=255)
                                                       gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                             >  1:1133757/1‑100 (MQ=255)
                                                       ctaacacccgcCCGAAGTGTACCGTTTGCCTTCTGCCCTGCCAGCAGCTA‑TGACAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCACGATAGAca                                             >  2:717979/9‑100 (MQ=255)
                                                         aaGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagc                                           <  1:2761272/100‑1 (MQ=255)
                                                          aGGCACGCCGGAAGTGTGCCGTTTGCCTTATGCCCTGCCAGCAGCTG‑TGGCAAGATGCGGGTGCGGATCTTGGTCATACCGTTGAGACAGATAGCcagca                                          <  2:573349/100‑1 (MQ=255)
                                                           ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                         >  1:2583043/1‑100 (MQ=255)
                                                             cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACagcagcc                                       >  1:2562381/1‑99 (MQ=255)
                                                             cgcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCt                                       >  2:2181374/1‑100 (MQ=255)
                                                             cgcgcCGAAAGTGTGCCGTTTTTCTTCTCATCTGCCAGTAGCTG‑TGGCAGGATGAGGGTGCGGACCTTGGTCATACCGTTGAGCGCGATAGACagcagc                                        <  1:660191/99‑1 (MQ=255)
                                                                 ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGAt                                    >  1:487626/1‑99 (MQ=255)
                                                                   ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCt                                 <  1:460498/100‑1 (MQ=255)
                                                                   ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCt                                 <  1:824295/100‑1 (MQ=255)
                                                                       gtgtGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAAt                             >  1:1105906/1‑100 (MQ=255)
                                                                         gtgCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                           >  1:1759250/1‑100 (MQ=255)
                                                                             cGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAACGCGATAGACAGCAGCTGATGCTTAATATACgg                       >  1:3759193/1‑100 (MQ=255)
                                                                              gTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACggg                      <  2:1732278/100‑1 (MQ=255)
                                                                              gTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACggg                      <  1:885548/100‑1 (MQ=255)
                                                                                 tGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTa                   >  1:1090331/1‑100 (MQ=255)
                                                                                  cacccgcatcctgccA‑CAGCTGCTGGCAGGATGCGGCTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTAc                  >  1:2763515/10‑100 (MQ=37)
                                                                                   ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                 >  1:1249833/1‑100 (MQ=255)
                                                                                   acccgcatcctgccA‑CAGCTGCTGGCAGGAGGCAGGTGCGGAACTTCGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGGACGGGTTACg                 >  2:2580425/9‑100 (MQ=25)
                                                                                    nttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                 >  1:834206/2‑99 (MQ=255)
                                                                                     ttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa               <  2:3890290/100‑1 (MQ=255)
                                                                                     ttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa               >  2:2140918/1‑100 (MQ=255)
                                                                                       gcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa             >  1:99055/5‑100 (MQ=255)
                                                                                        tgccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            <  2:3801907/100‑1 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGCTAGACAGCATCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:2999161/4‑100 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:3615373/4‑100 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:2905517/4‑100 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGATGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  2:2571179/4‑100 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGACAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:3809379/4‑100 (MQ=255)
                                                                                           cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc         >  1:1481840/1‑100 (MQ=255)
                                                                                           cctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc         >  2:109930/1‑100 (MQ=255)
                                                                                              gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg      <  2:813172/100‑1 (MQ=255)
                                                                                               ccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGt     <  1:2140937/100‑1 (MQ=255)
                                                                                                cagcagCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTg    >  1:3280893/1‑100 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  >  2:2932460/1‑100 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  >  1:2603848/1‑100 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  <  1:2080746/100‑1 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  >  2:2872778/1‑100 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  >  2:307313/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
TGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  NC_000913/1609457‑1609654

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