Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,441 C→G 5.4% T234T (ACC→ACG 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,4410CG 5.4% 258.3 / 5.2 74T234T (ACC→ACGcbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base G (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCT  >  NC_000913/1038293‑1038589
                                                                                                                                                    |                                                                                                                                                    
tgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAcc                                                                                                                                                      >  2:253673/1‑149 (MQ=255)
 gCACACCGTAATGGCGGGCTACGTGACCGGGGCAATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGGGCGGTCGGGAGCGCCATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGAACTCTGGAGATTATCGGGAccc                                                                                                                                                     <  4:683843/149‑1 (MQ=255)
   acacCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTg                                                                                                                                                   >  1:780375/1‑149 (MQ=255)
          aaTGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCg                                                                                                                                            >  2:598682/1‑149 (MQ=255)
            tGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGga                                                                                                                                          <  2:816527/149‑1 (MQ=255)
            tGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGga                                                                                                                                          <  1:642041/149‑1 (MQ=255)
                 ggCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAg                                                                                                                                     >  3:551155/1‑149 (MQ=255)
                     aCGTGACCGGGGCTATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTCTCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCt                                                                                                                                 <  4:682395/149‑1 (MQ=255)
                        tGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCg                                                                                                                              <  1:642227/149‑1 (MQ=255)
                         gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGt                                                                                                                             >  4:403304/1‑149 (MQ=255)
                           ccGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTAt                                                                                                                           >  1:806928/1‑149 (MQ=255)
                            cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                                                                          <  3:555584/149‑1 (MQ=255)
                            cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                                                                          <  2:111848/149‑1 (MQ=255)
                            cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                                                                          >  1:634450/1‑149 (MQ=255)
                              gggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGaa                                                                                                                        <  1:253673/149‑1 (MQ=255)
                                   aTGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTcgc                                                                                                                   <  1:598682/149‑1 (MQ=255)
                                   aTGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTcgc                                                                                                                   >  2:722168/1‑149 (MQ=255)
                                    tGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTcgcg                                                                                                                  >  2:803372/1‑149 (MQ=255)
                                      tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCa                                                                                                                <  4:786194/149‑1 (MQ=255)
                                          ttatGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTa                                                                                                            >  3:161624/1‑149 (MQ=255)
                                              ggCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAAc                                                                                                        >  2:54169/1‑149 (MQ=255)
                                                 gATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGga                                                                                                                                          >  1:132875/1‑112 (MQ=255)
                                                 gATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGga                                                                                                                                          <  2:132875/112‑1 (MQ=255)
                                                 gATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCgg                                                                                                     >  1:239008/1‑149 (MQ=255)
                                                  aTCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGt                                                                                                    <  4:274929/149‑1 (MQ=255)
                                                   tCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCaa                                                                                                               <  3:412247/137‑1 (MQ=255)
                                                   tCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCaa                                                                                                               >  4:412247/1‑137 (MQ=255)
                                                      gcgcCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTaaaa                                                                                                <  3:84740/149‑1 (MQ=255)
                                                      gcgcCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTaaaa                                                                                                <  1:276069/149‑1 (MQ=255)
                                                           tGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTggc                                                                                           <  1:803372/149‑1 (MQ=255)
                                                              tATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACa                                                                                                                                      <  1:156487/103‑1 (MQ=255)
                                                              tATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACa                                                                                                                                      >  2:156487/1‑103 (MQ=255)
                                                                tttACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGa                                                                                      <  2:407093/149‑1 (MQ=255)
                                                                 ttACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGAt                                                                                     <  2:489350/149‑1 (MQ=255)
                                                                        gcgGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAggg                                                                              <  4:161624/149‑1 (MQ=255)
                                                                            aaacgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGt                                                                          >  1:330180/5‑149 (MQ=255)
                                                                             cGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTg                                                                         <  3:366921/149‑1 (MQ=255)
                                                                             cGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTg                                                                         <  4:51902/149‑1 (MQ=255)
                                                                              gggAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACGCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg                                                                        >  2:475806/1‑149 (MQ=255)
                                                                              gggAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACGCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg                                                                        <  3:34216/149‑1 (MQ=255)
                                                                              gggAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACGCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg                                                                        >  4:34216/1‑149 (MQ=255)
                                                                              gggAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACGCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg                                                                        <  1:475806/149‑1 (MQ=255)
                                                                                gAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCa                                                                      >  2:35457/1‑149 (MQ=255)
                                                                                   cgcAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAAc                                                                   <  2:806928/149‑1 (MQ=255)
                                                                                       aTGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGaa                                                               <  4:551155/149‑1 (MQ=255)
                                                                                        tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAac                                                              <  2:555649/149‑1 (MQ=255)
                                                                                        tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACGCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCTTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAac                                                              <  2:780375/149‑1 (MQ=255)
                                                                                         gTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc                                                             >  1:23086/1‑149 (MQ=255)
                                                                                          tCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacct                                                            >  2:350408/1‑149 (MQ=255)
                                                                                            gCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctgc                                                          >  4:558147/1‑149 (MQ=255)
                                                                                                   gCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAAcctgcacctgca                                                   >  4:655562/1‑149 (MQ=255)
                                                                                                     tCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                                                                                     >  2:58247/1‑49 (MQ=255)
                                                                                                     tCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                                                                                     <  1:58247/49‑1 (MQ=255)
                                                                                                             cATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCAtg                                         >  4:218865/1‑149 (MQ=255)
                                                                                                                 ggTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGt                                     <  3:558147/149‑1 (MQ=255)
                                                                                                                   ttCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTg                                   >  3:745715/1‑149 (MQ=255)
                                                                                                                   ttCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTg                                   >  2:246863/1‑149 (MQ=255)
                                                                                                                       gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTc                                                                                                                     >  4:687098/1‑63 (MQ=255)
                                                                                                                       gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTc                                                                                                                     <  3:687098/63‑1 (MQ=255)
                                                                                                                       gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGaa                                                               >  3:407028/1‑117 (MQ=255)
                                                                                                                       gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGaa                                                               <  4:407028/117‑1 (MQ=255)
                                                                                                                       gTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTg                               >  4:279015/1‑149 (MQ=255)
                                                                                                                          ttCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGcc                            >  4:107010/1‑149 (MQ=255)
                                                                                                                            cGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCgg                          >  4:231546/1‑149 (MQ=255)
                                                                                                                             ggTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGa                         >  2:35455/1‑149 (MQ=255)
                                                                                                                              gTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGaa                        >  3:390861/1‑149 (MQ=255)
                                                                                                                                    tgacGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGACGATGGAAGGGGAGTGGCAAACGGAacct                                                            <  2:629184/104‑1 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacct                                                            >  1:629184/1‑107 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGa                   >  2:674754/1‑148 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGa                   <  1:674754/148‑1 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGAt                  >  2:249485/1‑149 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGAt                  >  4:458770/1‑149 (MQ=255)
                                                                                                                                    tGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGAt                  <  3:403304/149‑1 (MQ=255)
                                                                                                                                     ggCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATc                 >  3:357826/1‑149 (MQ=255)
                                                                                                                                        gATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCt                                <  4:506360/131‑1 (MQ=255)
                                                                                                                                        gATTATCGGTACCCTGCAACTCGGAGACAGTTCGGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCt                                >  3:506360/1‑131 (MQ=255)
                                                                                                                                         aTTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAaga             <  2:239008/149‑1 (MQ=255)
                                                                                                                                          ttATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAagag            >  1:267348/1‑149 (MQ=255)
                                                                                                                                                    ccTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCt  >  3:171000/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCT  >  NC_000913/1038293‑1038589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: