Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,776,799 C→A 9.0% R354S (CGT→AGT)  betA‑II → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,776,7990CA 9.0% 311.1 / 13.1 100R354S (CGT→AGT) betA‑IIcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (41/49);  new base A (5/4);  total (47/53)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

ACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAA  >  NC_002947/5776653‑5776944
                                                                                                                                                  |                                                                                                                                                  
aCCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTATGGTGGAACCAGCCGGCGATCTGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTc                                                                                                                                                 <  1:877371/150‑1 (MQ=255)
 ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTcc                                                                                                                                                <  2:2237415/150‑1 (MQ=255)
  cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCg                                                                                                                                               >  2:2221712/1‑150 (MQ=255)
   tGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCgg                                                                                                                                              >  2:4848/1‑150 (MQ=255)
    gCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGa                                                                                                                                             >  1:2233127/1‑150 (MQ=255)
     cAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACGAGTCCGGat                                                                                                                                            <  1:2290651/150‑2 (MQ=255)
          cGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       <  2:1351016/150‑1 (MQ=255)
           gCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAg                                                                                                                                      <  1:2425740/150‑1 (MQ=255)
               gCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGc                                                                                                                                  <  2:2012983/150‑1 (MQ=255)
                    cAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAAc                                                                                                                             >  1:239912/1‑150 (MQ=255)
                         ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCa                                                                                                                        >  1:584760/1‑150 (MQ=255)
                         ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGAGCCTAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCa                                                                                                                        >  2:1908091/1‑150 (MQ=255)
                          gTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  1:972987/1‑150 (MQ=255)
                                     cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            >  1:1516754/1‑150 (MQ=255)
                                     cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            <  1:4848/150‑1 (MQ=255)
                                     cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            >  1:1883543/1‑150 (MQ=255)
                                      ccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTg                                                                                                           <  2:239912/150‑1 (MQ=255)
                                       cGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAAAGATCGGAAGGCCGGGGGGGTGATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGc                                                                                                          >  1:2252328/1‑150 (MQ=255)
                                        gTCGCTGCGGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGcc                                                                                                         >  1:2038811/1‑150 (MQ=255)
                                            ctgctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGtt                                                                                                     <  2:2164213/150‑1 (MQ=255)
                                              gctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGc                                                                                                   <  2:1688273/150‑1 (MQ=255)
                                                tgtgGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 <  1:590429/150‑1 (MQ=255)
                                                  tggtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAtt                                                                                               >  1:1335900/1‑150 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTa                                                                                                                     <  1:818100/126‑1 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTa                                                                                                                     >  2:818100/1‑126 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTc                                                                                                              >  1:753011/1‑133 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTc                                                                                                              <  2:753011/133‑1 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                                                             >  2:2388996/1‑150 (MQ=255)
                                                    gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTCCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                                                             <  1:553312/150‑1 (MQ=255)
                                                     tggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                                                            <  2:741819/150‑1 (MQ=255)
                                                     tggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                                                            >  1:494171/1‑150 (MQ=255)
                                                          ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACaa                                                                                       <  2:1889129/150‑1 (MQ=255)
                                                          ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACaa                                                                                       >  1:1889129/1‑150 (MQ=255)
                                                          ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCC‑GGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAAc                                                                                      <  1:689978/150‑1 (MQ=255)
                                                           cAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTANAAc                                                                                      <  2:1516754/150‑1 (MQ=255)
                                                            aGCCGGCGATCGGTGCCGAGAGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGTCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACg                                                                                     >  2:742652/1‑150 (MQ=255)
                                                               cGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCt                                                                                  >  1:1524669/1‑150 (MQ=255)
                                                                ggCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTc                                                                                 <  1:239890/150‑1 (MQ=255)
                                                                   gATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCaa                                                                              >  2:2330696/1‑150 (MQ=255)
                                                                      cGGTGCCGAGCGGCAGTCCCATCGCACGGGGATCATCGCCAGCAAGCAGTTCTAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATgg                                                                           <  2:856756/150‑1 (MQ=255)
                                                                       ggTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGc                                                                          <  2:1524669/150‑1 (MQ=255)
                                                                        gTGCGGAGTGGTTGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCg                                                                         <  1:1216313/150‑1 (MQ=255)
                                                                         tGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGt                                                                        <  2:1614786/150‑1 (MQ=255)
                                                                            cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       <  1:1491281/100‑1 (MQ=255)
                                                                            cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  2:1491281/1‑100 (MQ=255)
                                                                            cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     <  2:2038811/150‑1 (MQ=255)
                                                                             gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     <  2:816454/149‑1 (MQ=255)
                                                                             gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     >  1:816454/1‑149 (MQ=255)
                                                                             gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    <  1:2449777/150‑1 (MQ=255)
                                                                             gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    >  2:18129/1‑150 (MQ=255)
                                                                              aGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGTCGGGGGGGTTCATGCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAgg                                                                   >  2:929957/1‑150 (MQ=255)
                                                                              aGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATCAACTACAACGGCTCCAATGGCGTGAAgg                                                                   <  1:423861/150‑1 (MQ=255)
                                                                                 ggATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAAc                                                                <  2:1784590/150‑1 (MQ=255)
                                                                                  gATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACa                                                               >  1:2006933/1‑150 (MQ=255)
                                                                                    tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             <  2:1883543/150‑1 (MQ=255)
                                                                                    tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             >  1:1996756/1‑150 (MQ=255)
                                                                                    tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             <  1:1285113/150‑1 (MQ=255)
                                                                                      ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGc                                                           <  1:67418/150‑1 (MQ=255)
                                                                                         aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       <  1:2155317/71‑1 (MQ=255)
                                                                                         aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       >  2:2155317/1‑71 (MQ=255)
                                                                                         aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc                                                        >  1:1006839/1‑150 (MQ=255)
                                                                                           tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCa                                                      >  1:1032988/1‑150 (MQ=255)
                                                                                           tGGCGCCGATACCGGTGCCAGCAACCAGATCGAGGCCGGGGGGGGCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACAACCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCa                                                      >  2:588741/1‑150 (MQ=255)
                                                                                            ggCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCCCTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAg                                                     <  2:2252328/150‑1 (MQ=255)
                                                                                            ggCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAg                                                     <  1:2221712/150‑1 (MQ=255)
                                                                                             gCACCGGTATCGGCGCCAGCAACCAGTTCGAGCCCGGGGGGTTCATCCGGACTAGTTCGTAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAgg                                                    <  1:1206570/150‑1 (MQ=255)
                                                                                              cACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   >  2:1852082/1‑150 (MQ=255)
                                                                                              cACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   <  1:1852082/150‑1 (MQ=255)
                                                                                               aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGAATAATCCGGACTGGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcg                                                  >  1:1005103/1‑150 (MQ=255)
                                                                                                     aTCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATg                                            <  2:2006933/150‑1 (MQ=255)
                                                                                                      tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 >  2:956590/1‑96 (MQ=255)
                                                                                                      tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 <  1:956590/96‑1 (MQ=255)
                                                                                                         gtgcCAGCAACCGGTTCGGGGCCGGGGGGTTCATCCGGGCTCGTCCGGAGTTCAAGTGACCGGACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACGACGGCTCGAATGGCGTGAAGGAACATGGCTTCCAGTCGCACATG‑GGCt                                       <  2:1032988/148‑1 (MQ=255)
                                                                                                         gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     <  2:999219/121‑1 (MQ=255)
                                                                                                         gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     >  1:999219/1‑121 (MQ=255)
                                                                                                         gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCt                                       <  2:972987/150‑1 (MQ=255)
                                                                                                             cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCAt                                   <  2:1006839/150‑1 (MQ=255)
                                                                                                             cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCAt                                   >  2:468666/1‑150 (MQ=255)
                                                                                                                 aaCCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATgcgc                               >  2:508220/1‑150 (MQ=255)
                                                                                                                   ccAGTTCGAGGCAAGGGGGGTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTc                             >  2:1798323/1‑150 (MQ=255)
                                                                                                                     aGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           >  1:2436820/1‑150 (MQ=255)
                                                                                                                          gAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    <  1:302322/105‑1 (MQ=255)
                                                                                                                          gAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    >  2:302322/1‑105 (MQ=255)
                                                                                                                                 gggggTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGcc               >  2:2324827/1‑150 (MQ=255)
                                                                                                                                      ttCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATc          <  2:2436820/150‑1 (MQ=255)
                                                                                                                                        actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCa                                    >  1:709643/3‑122 (MQ=255)
                                                                                                                                        actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCa                                    <  2:709643/120‑1 (MQ=255)
                                                                                                                                        actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCa        >  2:1767917/3‑150 (MQ=255)
                                                                                                                                         ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       >  1:1277414/2‑150 (MQ=255)
                                                                                                                                         ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGC‑GATGAGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       >  1:1308061/2‑150 (MQ=255)
                                                                                                                                         aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  2:689484/150‑1 (MQ=255)
                                                                                                                                         aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  1:18129/150‑1 (MQ=255)
                                                                                                                                         aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  2:1861348/150‑1 (MQ=255)
                                                                                                                                          tCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           >  2:1708348/1‑129 (MQ=255)
                                                                                                                                          tCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           <  1:1708348/129‑1 (MQ=255)
                                                                                                                                             ggACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCa   >  2:1827976/1‑150 (MQ=255)
                                                                                                                                              gACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCaa  <  1:2324827/150‑1 (MQ=255)
                                                                                                                                              gACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCaa  <  1:1798323/150‑1 (MQ=255)
                                                                                                                                               aCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTAc                                                                                                                    <  2:2442880/36‑1 (MQ=255)
                                                                                                                                               aCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTAc                                                                                                                    >  1:2442880/1‑36 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAA  >  NC_002947/5776653‑5776944

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