Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,557 C→T 6.0% A273V (GCC→GTC)  cbdA → cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,5570CT 6.0% 210.7 / 2.5 66A273V (GCC→GTC) cbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAA  >  NC_000913/1038411‑1038696
                                                                                                                                                  |                                                                                                                                           
cGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCt                                                                                                                                           >  6:153741/1‑149 (MQ=255)
cGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCt                                                                                                                                           <  5:154883/149‑1 (MQ=255)
 gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGACCCTGCACCTGCACCGTTCCATGTGGGTGCCTg                                                                                                                                          >  2:325285/1‑149 (MQ=255)
       ggTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCgg                                                                                                                                     >  4:214671/1‑148 (MQ=255)
        gTACTCTTGCGAGTAAATGAACACTTCAACACGAAGACAGAACTGCGTTAGAAGTCCCCCAACTACAACCGGTAAAACTGGCGGAGAAGGAAGCGGAGAGGCAATCGTAAACATCACCTGCACCGAGACAAGTGATTGCGAGGCCGGaa                                                                                                                                   >  4:319983/1‑149 (MQ=255)
              tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAgg                                                                                                                               <  4:395991/147‑1 (MQ=255)
               ggCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTCGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGAt                                                                                                                             >  3:340790/1‑148 (MQ=255)
                gCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATc                                                                                                                            <  8:397987/148‑1 (MQ=255)
                   aTTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCaa                                                                                                                          >  5:336893/1‑147 (MQ=255)
                       tCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGTCTGGCCGGAACAGGATCAAGAGCg                                                                                                                     >  7:232199/1‑148 (MQ=255)
                       tCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCg                                                                                                                     >  6:249769/1‑148 (MQ=255)
                       tCGGTACCCTGCAACTCGGAGACAGGTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGGGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGt                                                                                                                    <  7:405037/149‑1 (MQ=255)
                       tCGGTACCATGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGTCTGGCCGGAACAGGATCAAGAGCg                                                                                                                     >  7:226333/1‑148 (MQ=255)
                        cGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGt                                                                                                                    >  5:61921/1‑148 (MQ=255)
                           tACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACg                                                                                                                <  5:249769/149‑1 (MQ=255)
                              ccTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCt                                                                                                             >  6:180721/1‑149 (MQ=255)
                                  ctaCTCGGAGTCAGTTTTGGGTATGAAGTCGCGCAAGTTCTTCCGGTAAAACTGGCGGCGATGGATGGGGAGTGGCAAACGGAACCTGCAGGTGCACCGTTCCATGTGGTTGTCTGGCCGGAACAGGATCAAGTGCGTAACGCCttt                                                                                                           <  8:232199/145‑1 (MQ=255)
                                  cAACTCGGAGACAGTTCTTCGTATGAAGTCGCGAAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGTCTGGCCGGAACAGGATCAAGAGCGTAACTCCttt                                                                                                           <  8:226333/147‑1 (MQ=255)
                                      tCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCt                                                                                                      >  3:383774/1‑148 (MQ=255)
                                         gagaCAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCaa                                                                                                   >  1:346004/1‑148 (MQ=255)
                                              aGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTcc                                                                                             <  6:91405/149‑1 (MQ=255)
                                              aGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              >  4:248605/1‑148 (MQ=255)
                                              aGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              <  3:248605/148‑1 (MQ=255)
                                                     cGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCgctg                                                                                      >  8:64144/1‑149 (MQ=255)
                                                     cGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCgctg                                                                                      >  1:83714/1‑149 (MQ=255)
                                                         tGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAg                                                                                  <  1:325285/149‑1 (MQ=255)
                                                         tGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAg                                                                                  <  6:336893/149‑1 (MQ=255)
                                                           aaGTCGCGCAAGTACAACCGGTAAACCTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAgg                                                                                 <  6:128921/148‑1 (MQ=255)
                                                                gcgcAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCt                                                                           >  3:195556/1‑149 (MQ=255)
                                                                 cgcAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTc                                                                          <  5:180721/149‑1 (MQ=255)
                                                                  gcAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCg                                                                         >  2:354071/1‑149 (MQ=255)
                                                                      gTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcac                                                                     >  5:155875/1‑149 (MQ=255)
                                                                       taaaaCCGGTAAAACTGGCGGCGATGGAAGGGGTGTGGCAAACGGCACCTGCACCAGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGAATCCTCGCcact                                                                    >  8:22866/4‑149 (MQ=255)
                                                                         cAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcactca                                                                  <  2:346004/149‑1 (MQ=255)
                                                                          aaCCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcactcac                                                                 >  8:129610/1‑149 (MQ=255)
                                                                          aaCCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcactcac                                                                 >  8:62989/1‑149 (MQ=255)
                                                                          aaCCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcactca                                                                  >  8:399922/1‑148 (MQ=255)
                                                                           aCCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcact                                                                    >  8:346557/1‑145 (MQ=255)
                                                                           aCCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCcact                                                                    <  7:346557/145‑1 (MQ=255)
                                                                                   aaCTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGa                                                         >  1:238337/1‑148 (MQ=255)
                                                                                   aaCTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGAt                                                        <  1:395953/149‑1 (MQ=255)
                                                                                        gcggcgATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATaaa                                                     <  7:62989/147‑1 (MQ=255)
                                                                                             gATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCGTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGc                                              <  6:155875/149‑1 (MQ=255)
                                                                                                  acggggATTGGCACACGGACCCTTCACCTGCACAGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGATGCTAGGGATCCTCGACACTCACTCATTAGATAAACCCGTGCCgtgt                                         <  2:408609/147‑1 (MQ=255)
                                                                                                  aaGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGt                                         <  5:153741/149‑1 (MQ=255)
                                                                                                            gCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAAtt                               >  1:187217/1‑149 (MQ=255)
                                                                                                             cAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAAttt                              <  7:129610/149‑1 (MQ=255)
                                                                                                             cAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCAATCACTCATTAGATAACCCCGTGCCGGGTCTGAAGAAttt                              <  4:340790/149‑1 (MQ=255)
                                                                                                              aaaCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              <  1:44376/84‑1 (MQ=255)
                                                                                                              aaaCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              >  2:44376/1‑84 (MQ=255)
                                                                                                              aaaCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              >  4:224151/1‑84 (MQ=255)
                                                                                                              aaaCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTc                                                                                              <  3:224151/84‑1 (MQ=255)
                                                                                                               aaCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGTCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCGGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAAACCGTGCCGGGTCTGAAGAATTTGa                            <  3:319983/149‑1 (MQ=255)
                                                                                                                aCGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGAt                           <  3:214671/149‑1 (MQ=255)
                                                                                                                 cGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATg                          <  1:354071/149‑1 (MQ=255)
                                                                                                                    aaCCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCt                       <  4:383774/149‑1 (MQ=255)
                                                                                                                       ctgcacctgcacCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGa                     >  3:68964/1‑148 (MQ=255)
                                                                                                                          cacctgcaccGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAAc                  <  4:195556/148‑1 (MQ=255)
                                                                                                                              tgcaccGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTAcc             >  6:386031/1‑149 (MQ=255)
                                                                                                                                caccGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCa           <  2:187217/149‑1 (MQ=255)
                                                                                                                                caccGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGAATCCTCGCCACTCACTCATTAGATAAACCCGGGCCGGGTCTGAAGAATTTGATGGCTGACACCTAcc             <  7:22866/147‑1 (MQ=255)
                                                                                                                                  ccGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACg         >  5:149463/1‑149 (MQ=255)
                                                                                                                                   cGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACACACGc        >  7:204596/1‑149 (MQ=255)
                                                                                                                                         aTGTGGTTGCCTGGCCGGAACAGGATCAAGATCGTAACGCCTTTGCCCTCAAAATTCAGGCGCTGCTAGGGATCCACGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGATTTTGATGGCTGAAACCTCCCCACGCTTGCaa  >  4:266257/1‑149 (MQ=255)
                                                                                                                                             ggTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCTTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGaaa                   <  1:148723/128‑1 (MQ=255)
                                                                                                                                             ggTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCATTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGaaa                   >  2:148723/1‑128 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAA  >  NC_000913/1038411‑1038696

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