Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,099 C→T 6.0% G237G (GGG→GGA)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,0990CT 6.0% 211.1 / 2.5 67G237G (GGG→GGA) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGC  >  NC_000913/1880951‑1881245
                                                                                                                                                    |                                                                                                                                                  
ttATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGc                                                                                                                                                    <  3:293771/149‑1 (MQ=255)
      gCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATc                                                                                                                                               >  4:35136/1‑148 (MQ=255)
      gCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATc                                                                                                                                               >  8:113843/1‑148 (MQ=255)
       cACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCga                                                                                                                                             <  5:144589/149‑1 (MQ=255)
          tttCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCCCCAGGAACGCAGAGCAATCCCCTGTCCATCgcgcc                                                                                                                                          >  2:226772/1‑145 (MQ=255)
          tttCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGTCCATCgaga                                                                                                                                           >  8:503803/1‑148 (MQ=255)
                   cGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGCACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAt                                                                                                                                 <  4:277210/149‑1 (MQ=255)
                   cGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAt                                                                                                                                 <  1:275445/149‑1 (MQ=255)
                           gCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGTCCATCGAGACACCACTGATGGACAAtt                                                                                                                         <  1:226772/149‑1 (MQ=255)
                           gCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGTCCATCGAGACACCACTGATGGACAAtt                                                                                                                         <  7:503803/149‑1 (MQ=255)
                            ccTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAAttt                                                                                                                        <  3:35136/149‑1 (MQ=255)
                             cTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTACCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtc                                                                                                                       >  6:13446/1‑149 (MQ=255)
                                gAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtctc                                                                                                                     <  3:447696/148‑1 (MQ=255)
                                gAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtctc                                                                                                                     >  8:272475/1‑148 (MQ=255)
                                gAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtctc                                                                                                                     <  8:389164/148‑1 (MQ=255)
                                gAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtct                                                                                                                      <  1:467929/147‑1 (MQ=255)
                                 aaaCATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCa                                                                                                                   >  1:208164/1‑149 (MQ=255)
                                     aTAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGAAGCACATGCACTAAAAGACCACAGCCAATGATTTCGCAAACATCCCACCAGCACCGCACAGCAATCCCCTGCCAATCGAGACACCACTGAAGGACAAATACTCCATGAt                                                                                                               >  8:361952/1‑149 (MQ=255)
                                      tAAACTGCCCAGACGTTCGCTGGCTTGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCATGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGAtt                                                                                                              <  5:121218/149‑1 (MQ=255)
                                       aaaCGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTa                                                                                                             <  4:181230/149‑1 (MQ=255)
                                       aaaCGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGAAATCCCCTGCCCATCGAGCCAACACTGATGGACAATTTCTCcatcat                                                                                                               >  6:386458/1‑147 (MQ=255)
                                         aCGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                                                                           >  5:283766/1‑149 (MQ=255)
                                         aCGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                                                                           <  7:113843/149‑1 (MQ=255)
                                         aCGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAga                                                                                                           <  2:45108/149‑1 (MQ=255)
                                              ccAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATgg                                                                                                                               >  3:86275/1‑124 (MQ=255)
                                              ccAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATgg                                                                                                                               <  4:86275/124‑1 (MQ=255)
                                                       cGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCg                                                                                              <  6:283766/148‑1 (MQ=255)
                                                          tggctggTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGCGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGa                                                                                          >  4:168/1‑149 (MQ=255)
                                                             ctggTAATATTCCGGTAGCACTTTCACTAAATGACCACTGCCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                                       <  5:386458/149‑1 (MQ=255)
                                                                gTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                                                    >  4:195964/1‑149 (MQ=255)
                                                                gTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTaa                                                                                    >  2:223979/1‑149 (MQ=255)
                                                                  gatatTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                                                  <  4:360265/148‑1 (MQ=255)
                                                                      ttCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCAt                                                                              <  1:500131/149‑1 (MQ=255)
                                                                      ttCCGGTAGCACTTGCACTAAATGACCACTGGCAATGATTTCGCAAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAAATTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCAt                                                                              >  1:287410/1‑149 (MQ=255)
                                                                         cGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggg                                                                           <  3:168/149‑1 (MQ=255)
                                                                                       cTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                                             >  4:120866/1‑149 (MQ=255)
                                                                                        tAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGCCACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGc                                                            >  7:459207/1‑149 (MQ=255)
                                                                                         aaaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAt                                                                                                                                 >  4:85495/1‑79 (MQ=255)
                                                                                         aaaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGAt                                                                                                                                 <  3:85495/79‑1 (MQ=255)
                                                                                            tGAACACTGGCAATGTTTTCGCTAACATCCCACCAGGAACTGAGAGCAATCCCCTTCCCATCGAGACACCACTGATGGACAATTTCTCCATGTTTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTaa                                                        <  2:287410/149‑1 (MQ=255)
                                                                                              aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAtt                                                      <  5:13446/149‑1 (MQ=255)
                                                                                                cACTGGCAATGTTTTCGCTAACATCCCACCAGGCACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                                                    >  3:143172/1‑149 (MQ=255)
                                                                                                cACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                                                    >  2:59362/1‑149 (MQ=255)
                                                                                                  cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                                                                             <  1:254491/122‑1 (MQ=255)
                                                                                                  cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATg                                                                             >  2:254491/1‑122 (MQ=255)
                                                                                                     gCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGa                                                >  8:368822/1‑148 (MQ=255)
                                                                                                       aaTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCAAGacac                                             >  2:400052/1‑149 (MQ=255)
                                                                                                           ttttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                                       >  2:267080/1‑103 (MQ=255)
                                                                                                           ttttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                                       <  1:267080/103‑1 (MQ=255)
                                                                                                             ttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaaa                                        <  7:368822/148‑1 (MQ=255)
                                                                                                                   aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCAATGATTAGACGAAAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGCCACCAAAAGgatgat                                 >  1:276496/1‑149 (MQ=255)
                                                                                                                    aCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATc                                <  1:223979/149‑1 (MQ=255)
                                                                                                                        cccACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                                                    >  6:336428/1‑125 (MQ=255)
                                                                                                                        cccACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGc                                                    <  5:336428/125‑1 (MQ=255)
                                                                                                                          caccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc                          >  7:27689/1‑149 (MQ=255)
                                                                                                                           accaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                                              >  4:345967/1‑112 (MQ=255)
                                                                                                                            ccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCtt                        <  4:143172/149‑1 (MQ=255)
                                                                                                                               ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTaa                     <  3:120866/149‑1 (MQ=255)
                                                                                                                                    gCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAAt                    >  1:452645/1‑145 (MQ=255)
                                                                                                                                    gCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAAt                    <  2:452645/145‑1 (MQ=255)
                                                                                                                                     cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCa               <  7:272475/149‑1 (MQ=255)
                                                                                                                                          cAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCAc          <  1:59362/149‑1 (MQ=255)
                                                                                                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCa               <  3:239328/141‑1 (MQ=255)
                                                                                                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACgg        >  7:407463/1‑148 (MQ=255)
                                                                                                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACgg        >  5:105163/1‑148 (MQ=255)
                                                                                                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCGTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCa               >  4:239328/1‑141 (MQ=255)
                                                                                                                                             tCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGATGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACgg        >  8:238978/1‑148 (MQ=255)
                                                                                                                                                ccTGCCCATCGAGACACCACTGTTGGACTATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAAt    <  7:361952/149‑1 (MQ=255)
                                                                                                                                                    cccATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATg   >  7:167257/1‑146 (MQ=255)
                                                                                                                                                    cccATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGGGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGc  <  1:253779/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCAGACGTTCGCTGGCTGGTAATATTCCGGTAGCACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTTTAATGACCAGGCACGGCAATGC  >  NC_000913/1880951‑1881245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑