Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,555 Δ1 bp 18.4% 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.18.4% 325.8 / 76.3 114coding (1126/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (35/57);  new base . (9/12);  total (44/70)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

TGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  NC_000913/1609457‑1609654
                                                                                                  |                                                                                                    
tGTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagc                                                                                                     >  1:680628/1‑100 (MQ=255)
 gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagca                                                                                                    >  1:505850/1‑100 (MQ=255)
  tttCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagatc                                                                                                   >  1:583977/1‑97 (MQ=255)
      gCCGTTACGCTCACCACGATAGAACGCAATCAATGCTGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑t                                                                                              >  2:1902032/1‑100 (MQ=255)
       ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                             >  1:3999404/1‑100 (MQ=255)
       ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                             <  1:4038341/100‑1 (MQ=255)
            aCGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgg                                                                                        <  1:1302375/100‑1 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                       >  1:1499916/1‑100 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                       <  1:3410823/100‑1 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:3859239/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  1:717952/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:148897/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  1:2790478/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:55817/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:3879811/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  2:2570025/100‑2 (MQ=255)
             cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCATTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                       >  2:414234/1‑100 (MQ=255)
              gCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggc                                                                                      <  1:3088506/100‑3 (MQ=255)
               cTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATg                                                                                     <  1:230252/100‑1 (MQ=255)
               cTCACCACGATAGAACGCAATCAAGGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCAGCCA‑CAGCTGCTGGCAgggta                                                                                     <  2:3525735/100‑4 (MQ=255)
                tCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGc                                                                                    <  2:3504228/100‑1 (MQ=255)
                 caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCg                                                                                   <  1:2254366/100‑1 (MQ=255)
                 caccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCg                                                                                   <  2:2121966/100‑1 (MQ=255)
                  accacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCgg                                                                                  <  2:1914436/100‑1 (MQ=255)
                  accacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAgggcagaa                                                                                  <  1:1830773/100‑7 (MQ=38)
                   ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 <  1:2670608/100‑1 (MQ=255)
                     acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCATCAGCTG‑TGGCAGGATGCGGGTg                                                                               <  2:2432687/100‑1 (MQ=255)
                     acGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTg                                                                               <  1:3792882/100‑1 (MQ=255)
                      cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGc                                                                              <  1:2533166/100‑1 (MQ=255)
                        aTAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCgg                                                                            >  2:528903/1‑100 (MQ=255)
                             aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACtt                                                                       <  1:2756817/100‑1 (MQ=255)
                             aCGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACtt                                                                       >  1:3181736/1‑100 (MQ=255)
                               gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                     <  2:3975675/100‑1 (MQ=255)
                               gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                     <  1:2113415/100‑1 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    <  1:257109/100‑1 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    >  2:2960354/1‑100 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    <  1:3519438/100‑1 (MQ=255)
                                caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    >  1:1281684/1‑100 (MQ=255)
                                  atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGGGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                  <  2:2720639/100‑1 (MQ=255)
                                      atGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                              >  2:3393639/1‑100 (MQ=255)
                                      atGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                              <  2:3504329/100‑1 (MQ=255)
                                      atGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATAcc                                                              <  1:3253427/100‑1 (MQ=255)
                                          cGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTg                                                          <  1:2365592/100‑1 (MQ=255)
                                            cTAATGCGAAAGTAAGGCTCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAg                                                        <  1:3347271/100‑1 (MQ=255)
                                              aaTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                      <  1:3371935/100‑1 (MQ=255)
                                              aaTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                      <  2:2915608/100‑1 (MQ=255)
                                              aaTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcg                                                      <  1:481943/100‑1 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    <  2:606611/100‑1 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:3928663/1‑100 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:3663488/1‑100 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:2139840/1‑100 (MQ=255)
                                                tGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:1591552/1‑100 (MQ=255)
                                                  cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                  >  1:591179/1‑100 (MQ=255)
                                                  cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                  <  1:1591556/100‑1 (MQ=255)
                                                  cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGAt                                                  <  2:653632/100‑1 (MQ=255)
                                                   gAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                 >  1:915060/1‑100 (MQ=255)
                                                   gAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATa                                                 <  2:1740905/100‑1 (MQ=255)
                                                    aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                <  1:3663507/100‑1 (MQ=255)
                                                       gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                             <  2:3178595/100‑1 (MQ=255)
                                                       gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                             <  2:2431250/100‑1 (MQ=255)
                                                           ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTTCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                         <  2:4001239/100‑1 (MQ=255)
                                                           ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                         <  2:3724657/100‑1 (MQ=255)
                                                           ggCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAcagcag                                         <  1:1801033/100‑1 (MQ=255)
                                                               cgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCCGCTGa                                     >  2:1940563/1‑100 (MQ=255)
                                                                 ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                   >  1:3662824/1‑100 (MQ=255)
                                                                 ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                   <  1:2242382/100‑1 (MQ=255)
                                                                 ccGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATg                                   >  2:3835973/1‑100 (MQ=255)
                                                                   ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCt                                 <  2:4089615/100‑1 (MQ=255)
                                                                      aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                              <  1:2404839/100‑1 (MQ=255)
                                                                       gtgtGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGATGAGCGCGATAGACAGCAGCTGATGCTTa                               >  2:4012966/1‑98 (MQ=255)
                                                                          tgCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGt                           <  2:3721294/99‑1 (MQ=255)
                                                                            ccGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACg                        <  1:1722930/100‑1 (MQ=255)
                                                                             cGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                       >  1:1500739/1‑100 (MQ=255)
                                                                             cGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACgg                       >  1:487057/1‑100 (MQ=255)
                                                                              gTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCTCGATAGACAGCAGCTTATGCTTAATGTACggg                      <  1:3500980/100‑1 (MQ=255)
                                                                               tttGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACg                        <  2:1597536/97‑1 (MQ=255)
                                                                                 tGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACGTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTa                   >  2:2318687/1‑100 (MQ=255)
                                                                                   ccTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACg                 >  1:2956120/1‑100 (MQ=255)
                                                                                    cccgcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACgg                >  1:3320210/8‑100 (MQ=38)
                                                                                     ttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa               >  1:4155237/1‑100 (MQ=255)
                                                                                       gcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa             >  1:775178/5‑100 (MQ=255)
                                                                                        tggcctgccAGCAGCTG‑TGGCAGGAGGCGGGTGCGGGACTTGGTCATACCGGTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            <  2:3902160/97‑1 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:2065203/4‑100 (MQ=255)
                                                                                        catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:1044551/4‑100 (MQ=255)
                                                                                        caccctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:3965045/3‑100 (MQ=255)
                                                                                         gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa             <  2:902960/98‑1 (MQ=255)
                                                                                         gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg           <  1:248450/100‑1 (MQ=255)
                                                                                          tcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc          >  1:87574/2‑100 (MQ=255)
                                                                                          ccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGc          >  1:278536/1‑100 (MQ=255)
                                                                                           cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc         >  2:2635118/1‑100 (MQ=255)
                                                                                           cctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGcc         >  1:1691203/1‑100 (MQ=255)
                                                                                             tgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCg       <  2:1350124/100‑1 (MQ=255)
                                                                                             tgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCg       >  1:2618270/1‑100 (MQ=255)
                                                                                              gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg      <  2:278536/100‑1 (MQ=255)
                                                                                              gccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCgg      <  2:1158306/100‑1 (MQ=255)
                                                                                                cagcagCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTg    >  2:1164802/1‑100 (MQ=255)
                                                                                                cagcagCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTg    <  1:154093/100‑1 (MQ=255)
                                                                                                cagcagCTG‑TGGCAGGAGGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTg    <  2:2308024/100‑1 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGTGTTACGGAAACGCCCGGTGAc  <  2:487954/100‑1 (MQ=255)
                                                                                                  gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAc  >  1:1093396/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 < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑