Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,145 T→C 7.7% L35P (CTC→CCC)  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,1450TC 7.7% 228.0 / 7.4 78L35P (CTC→CCC) cydBcytochrome bd‑I ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base T (46/26);  new base C (2/4);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTC  >  NC_000913/772999‑773292
                                                                                                                                                  |                                                                                                                                                   
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                                                                       >  1:80602/1‑150 (MQ=255)
                                                                           gtggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCa                                                                       <  2:1499644/150‑1 (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                        <  1:306189/150‑1 (MQ=255)
                                                                                                                          gTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGa                        >  2:1081393/1‑150 (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)
                                                                                                                                            gCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGt      >  1:1168484/1‑150 (MQ=255)
                                                                                                                                            gCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGt      >  1:1162634/1‑150 (MQ=255)
                                                                                                                                            gCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGt      >  2:347390/1‑150 (MQ=255)
                                                                                                                                               tGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtct   >  2:1144561/1‑150 (MQ=255)
                                                                                                                                               tGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtct   >  2:156637/1‑150 (MQ=255)
                                                                                                                                                gCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTAGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtctc  <  1:1392060/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTC  >  NC_000913/772999‑773292

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