Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,138 G→A 8.3% G33S (GGC→AGC)  cydB → cytochrome bd‑I ubiquinol oxidase subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,1380GA 8.3% 208.7 / 7.1 72G33S (GGC→AGC) cydBcytochrome bd‑I ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base G (39/27);  new base A (2/4);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

TGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCG  >  NC_000913/772990‑773287
                                                                                                                                                    |                                                                                                                                                     
tGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTggg                                                                                                                                                      >  1:93592/1‑150 (MQ=255)
tGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTggg                                                                                                                                                      <  2:72454/150‑1 (MQ=255)
  aGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCa                                                                                                                                                    <  2:783462/150‑1 (MQ=255)
  aGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGGTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCa                                                                                                                                                    <  1:1304246/150‑1 (MQ=255)
     aGTCTTCCACGACTACTCAGCCGGCACGCTAAGAGAGGAGTCGTCAAATGGTCGATTATTAAGTTTTGCGTTTTATCTGGTGGGTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGc                                                                                                                                                 <  2:1300821/150‑1 (MQ=255)
         ttCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAc                                                                                                                                             >  2:248332/1‑150 (MQ=255)
          tCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCTTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCTcc                                                                                                                                            <  2:747078/150‑1 (MQ=255)
           ccACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCAccc                                                                                                                                           >  1:952824/1‑150 (MQ=255)
            cACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCg                                                                                                                                          <  2:733775/150‑1 (MQ=255)
              cGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGGCGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                        <  2:1594423/150‑1 (MQ=255)
              cGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAccccca                                                                                                                                        <  1:872740/150‑4 (MQ=255)
               gACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccccat                                                                                                                                       <  1:1601180/150‑5 (MQ=255)
                actactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccccatg                                                                                                                                      <  1:1217643/150‑6 (MQ=255)
                  tactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccccatgtc                                                                                                                                    <  1:248332/150‑8 (MQ=255)
                    ctCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCg                                                                                                                                  >  1:1590255/1‑150 (MQ=255)
                      cAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCCTGCTCACCCGTTTCCTCGGt                                                                                                                                >  1:294308/1‑150 (MQ=255)
                      cAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                <  2:947899/150‑1 (MQ=255)
                       aGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                               <  2:14651/150‑1 (MQ=255)
                               cGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGac                                                                                                                       >  2:394768/1‑150 (MQ=255)
                               cGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGac                                                                                                                       >  2:562455/1‑150 (MQ=255)
                                  tAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACAcc                                                                                                                    >  1:163892/1‑150 (MQ=255)
                                    aggcAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTTACGGTTTCGACATGGGGGTAGGCATGCTCACCCGTTTCCTAGGTCGTAACGACACAGa                                                                                                                  >  2:694188/4‑150 (MQ=255)
                                     gACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAg                                                                                                                 >  2:242657/1‑150 (MQ=255)
                                      aCAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCGTTTCCTCGGTCGTAACGCCAACGAGCg                                                                                                               >  2:70782/1‑150 (MQ=255)
                                         ggAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTc                                                                                                             >  1:781016/1‑150 (MQ=255)
                                          gAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCg                                                                                                            >  1:53146/1‑150 (MQ=255)
                                               gtcAAAAGATCGATTATGAAGTATTGCGTTTAATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTa                                                                                                       <  2:1600181/150‑1 (MQ=255)
                                                    aTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAtt                                                                                                  <  2:516581/150‑1 (MQ=255)
                                                          gATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTcc                                                                                            >  1:1598902/1‑150 (MQ=255)
                                                             tATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCAtt                                                                                         <  2:591665/150‑1 (MQ=255)
                                                              aTGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTg                                                                                        <  1:635909/150‑1 (MQ=255)
                                                               tgtaGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGc                                                                                       <  1:694188/147‑1 (MQ=255)
                                                                gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGca                                                                                      >  2:291787/1‑150 (MQ=255)
                                                                gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTACAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGca                                                                                      <  2:334710/150‑1 (MQ=255)
                                                                       tGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACt                                                                               >  2:620482/1‑150 (MQ=255)
                                                                         cGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTgg                                                                             <  1:985625/150‑1 (MQ=255)
                                                                            tttATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGAc                                                                          <  1:562455/150‑1 (MQ=255)
                                                                               aTCAGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGt                                                                       >  1:650287/1‑150 (MQ=255)
                                                                                tCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTa                                                                      >  2:252343/1‑150 (MQ=255)
                                                                                 cTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTaa                                                                     >  2:1055964/1‑150 (MQ=255)
                                                                                  tggtggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAAc                                                                    >  1:744305/1‑150 (MQ=255)
                                                                                   ggtggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAAcc                                                                   <  2:952824/150‑1 (MQ=255)
                                                                                    gtggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCa                                                                  <  2:1499644/150‑1 (MQ=255)
                                                                                    gtggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCa                                                                  >  1:80602/1‑150 (MQ=255)
                                                                                        ctgctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                              >  2:306189/1‑150 (MQ=255)
                                                                                         tgctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGttt                                                             >  1:1135178/1‑150 (MQ=255)
                                                                                          gctgGTTTGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTg                                                            <  2:1408766/150‑1 (MQ=255)
                                                                                             ggTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCt                                                         >  1:536371/1‑150 (MQ=255)
                                                                                              gTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTg                                                        <  2:781016/150‑1 (MQ=255)
                                                                                                tGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGAt                                                      >  2:1437026/1‑150 (MQ=255)
                                                                                                 ggCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATc                                                     <  2:53146/150‑1 (MQ=255)
                                                                                                   cGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCAc                                                   <  1:727868/150‑1 (MQ=255)
                                                                                                       ctgctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcgcg                                               >  2:980743/1‑150 (MQ=255)
                                                                                                            gATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                          >  1:219901/1‑150 (MQ=255)
                                                                                                            gATTGGTTTTGCAGTCACTGACGGTTTCGAAATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACAGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                          >  2:656194/1‑150 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGggcggc                                         >  2:1262651/1‑150 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGggcggc                                         <  1:980743/150‑1 (MQ=255)
                                                                                                               tGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGgcgc                                       <  1:1525830/150‑1 (MQ=255)
                                                                                                                 gTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCAc                                     >  1:1032805/1‑150 (MQ=255)
                                                                                                                    ttGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGGTTCCTCGGTCGTAACGACACCGAGCGGCGAATAATGATTAACTCCATTGCACCACACTGGGACGGAAACCAGGTTTGGCTGATCACCGCGGGCGGCGCActct                                  <  1:1615756/150‑1 (MQ=255)
                                                                                                                       cAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtg                               >  2:317279/1‑150 (MQ=255)
                                                                                                                         gTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgct                             >  1:273436/1‑150 (MQ=255)
                                                                                                                            aCTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCCCACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGcc                          >  1:1057040/1‑150 (MQ=255)
                                                                                                                                 cGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGcc                     <  1:79204/150‑1 (MQ=255)
                                                                                                                                  ggTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCg                    >  2:15790/1‑150 (MQ=255)
                                                                                                                                   gTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGa                   >  2:1081393/1‑150 (MQ=255)
                                                                                                                                   gTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGa                   <  1:306189/150‑1 (MQ=255)
                                                                                                                                   gTTTCGACAAGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGa                   >  2:821177/1‑150 (MQ=255)
                                                                                                                                        gACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCAACGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTc              >  1:110505/1‑150 (MQ=255)
                                                                                                                                        aggaaacgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCCCCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTc              >  2:723845/8‑150 (MQ=255)
                                                                                                                                          cATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTAGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCAATCTTTGCTGCCTGGCCGATGGTCTa            >  2:1392060/1‑150 (MQ=255)
                                                                                                                                            tGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATg          >  1:1237102/1‑150 (MQ=255)
                                                                                                                                             gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCATGTTTTGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGc         >  2:1538696/1‑150 (MQ=255)
                                                                                                                                              cgggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGcc        >  1:749187/2‑150 (MQ=255)
                                                                                                                                                   gggCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGc   <  1:774342/150‑1 (MQ=255)
                                                                                                                                                   gggCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGc   >  2:1411060/1‑150 (MQ=255)
                                                                                                                                                    ggCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCg  >  2:578103/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCG  >  NC_000913/772990‑773287

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