Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,492,696 A→C 6.4% T52P (ACC→CCC)  PP_2185 → formate dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,492,6960AC 6.4% 663.5 / 26.1 202T52P (ACC→CCC) PP_2185formate dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (101/88);  new base C (6/7);  total (107/95)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 203 total aligned reads displayed.

CAATTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGAAACCCCG  >  NC_002947/2492548‑2492845
                                                                                                                                                    |                                                                                                                                                     
cAATTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCAc                                                                                                                                                      <  1:5447575/150‑1 (MQ=255)
 aaTTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCAcc                                                                                                                                                     >  2:3738352/1‑150 (MQ=255)
 aaTTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCAcc                                                                                                                                                     >  1:707168/1‑150 (MQ=255)
  aTTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCa                                                                                                                                                    <  2:5337045/150‑1 (MQ=255)
  aTTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCa                                                                                                                                                    >  2:4431013/1‑150 (MQ=255)
   tttCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAg                                                                                                                                                   >  1:3709309/1‑150 (MQ=255)
    ttcttcGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGc                                                                                                                                                  <  2:1905549/150‑1 (MQ=255)
      cttcGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCAt                                                                                                                                                >  1:3272684/1‑150 (MQ=255)
       ttcGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGTGCCCAGCAtg                                                                                                                                               <  1:5414937/150‑2 (MQ=255)
        tcGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGTGCCCAGCAtgg                                                                                                                                              <  2:3218077/150‑3 (MQ=255)
          gACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATcccc                                                                                                                                            >  1:2315717/1‑150 (MQ=255)
            ccccGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCaa                                                                                                                                          >  1:1533778/1‑150 (MQ=255)
                   tCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTgc                                                                                                                                   <  1:629222/150‑1 (MQ=255)
                       aTAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGGCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCa                                                                                                                               <  1:4995760/150‑1 (MQ=255)
                             gcagcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACa                                                                                                                         <  2:152771/150‑1 (MQ=255)
                             gcagcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACa                                                                                                                         <  2:3679344/150‑1 (MQ=255)
                               agcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGc                                                                                                                       <  1:5100272/150‑1 (MQ=255)
                               agcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGc                                                                                                                       <  2:3272684/150‑1 (MQ=255)
                               agcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGc                                                                                                                       >  1:3490422/1‑150 (MQ=255)
                                gcGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGcc                                                                                                                      <  2:1589232/150‑1 (MQ=255)
                                     cTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGaa                                                                                                                 <  2:1079587/150‑1 (MQ=255)
                                     cTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGaa                                                                                                                 <  1:2148325/150‑1 (MQ=255)
                                       cGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGc                                                                                                               >  2:1377759/1‑150 (MQ=255)
                                        ggCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGcc                                                                                                              <  1:3343068/150‑1 (MQ=255)
                                          cACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCtt                                                                                                            <  2:4023557/150‑1 (MQ=255)
                                             cccGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCgg                                                                                                         >  2:4456146/1‑150 (MQ=255)
                                             cccGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCgg                                                                                                         >  1:4519211/1‑150 (MQ=255)
                                             cccGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCgg                                                                                                         >  2:4927983/1‑150 (MQ=255)
                                             cccGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCgg                                                                                                         <  2:3889352/150‑1 (MQ=255)
                                              ccGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGc                                                                                                        <  2:4618676/150‑1 (MQ=255)
                                                ggCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTc                                                                                                      >  2:2479528/1‑150 (MQ=255)
                                                  cccGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCt                                                                                                    >  2:3490598/1‑150 (MQ=255)
                                                  cccGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCt                                                                                                    <  1:4061931/150‑1 (MQ=255)
                                                   ccGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGTGCCCAGCAtggcgg                                                                                                                                           >  2:814818/1‑105 (MQ=255)
                                                       gAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGc                                                                                               <  1:4246231/150‑1 (MQ=255)
                                                                  gCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtggt                                                                                    <  1:2867230/150‑1 (MQ=255)
                                                                   cAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtggtg                                                                                   <  1:3856288/150‑1 (MQ=255)
                                                                      gTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGaa                                                                                <  1:3193608/150‑1 (MQ=255)
                                                                       tCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGaaa                                                                               <  2:2197264/150‑1 (MQ=255)
                                                                         aGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATc                                                                             <  1:2288136/150‑1 (MQ=255)
                                                                            cTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGaa                                                                          <  1:2598954/150‑1 (MQ=255)
                                                                            cTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGaa                                                                          <  1:1053967/150‑1 (MQ=255)
                                                                             tGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCCCACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAg                                                                         <  2:3673215/150‑1 (MQ=255)
                                                                             tGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAg                                                                         <  2:4207635/150‑1 (MQ=255)
                                                                              gAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAgg                                                                        <  1:1158390/150‑1 (MQ=255)
                                                                               aGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGc                                                                       <  2:4785180/150‑1 (MQ=255)
                                                                                 cATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCAGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCAt                                                                     >  2:5105469/1‑150 (MQ=255)
                                                                                   tCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGcc                                                                                                 >  1:1837270/1‑120 (MQ=255)
                                                                                   tCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGcc                                                                                                 <  2:1837270/120‑1 (MQ=255)
                                                                                          ccGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTa                                                            >  1:5152805/1‑150 (MQ=255)
                                                                                          aggcgcCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGc                                                                                                                       >  1:3636738/3‑91 (MQ=255)
                                                                                           cgcgCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTAc                                                           >  2:949361/1‑150 (MQ=255)
                                                                                            gcgcCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTAcc                                                          <  1:2818635/150‑1 (MQ=255)
                                                                                             cgcCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTAccc                                                         >  2:2928827/1‑150 (MQ=255)
                                                                                              gcCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCg                                                        <  1:1725486/150‑1 (MQ=255)
                                                                                                cATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGc                                                      >  1:466273/1‑150 (MQ=255)
                                                                                                 aTCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCa                                                     >  1:5286279/1‑150 (MQ=255)
                                                                                                 aTCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCa                                                     >  1:3380472/1‑150 (MQ=255)
                                                                                                    aGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGaa                                                                                                                 <  2:1221422/87‑1 (MQ=255)
                                                                                                    aGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATcc                                                  >  2:2792440/1‑150 (MQ=255)
                                                                                                    aGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACAAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATcc                                                  >  2:1466432/1‑150 (MQ=255)
                                                                                                     gCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCt                                                 >  1:1503150/1‑150 (MQ=255)
                                                                                                       gTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGc                                               <  2:1551848/150‑1 (MQ=255)
                                                                                                       gTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGc                                               >  1:4615661/1‑150 (MQ=255)
                                                                                                         gCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCACACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGcac                                             <  2:2295361/150‑1 (MQ=255)
                                                                                                         gCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtg                                                                                      >  1:3021126/1‑109 (MQ=255)
                                                                                                         gCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtg                                                                                      <  2:3021122/109‑1 (MQ=255)
                                                                                                         gCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCgg                                                               >  2:284193/1‑132 (MQ=255)
                                                                                                         gCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGcac                                             <  2:3143008/150‑1 (MQ=255)
                                                                                                          ccTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGcacc                                            <  1:898878/150‑1 (MQ=255)
                                                                                                                ggCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCa                                      <  1:4090262/150‑1 (MQ=255)
                                                                                                                ggCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCa                                      <  1:2849513/150‑1 (MQ=255)
                                                                                                                ggCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCa                                      <  1:4472580/150‑1 (MQ=255)
                                                                                                                  cACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGt                                    <  2:4367180/150‑1 (MQ=255)
                                                                                                                   aCCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtg                                                                                      <  2:669949/99‑1 (MQ=255)
                                                                                                                    ccTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCa                                  >  2:3959855/1‑150 (MQ=255)
                                                                                                                     cTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCGTCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCg                                                                                                >  1:1459867/1‑87 (MQ=255)
                                                                                                                     cTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCGTCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCg                                                                                                <  2:1459867/87‑1 (MQ=255)
                                                                                                                         gTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGaa                             <  2:4285119/150‑1 (MQ=255)
                                                                                                                          tCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAg                            <  1:3997537/150‑1 (MQ=255)
                                                                                                                                gcgcCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCAtgg                      >  1:2008224/1‑150 (MQ=255)
                                                                                                                                 cgccgccgcCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGGCAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCAtggt                     <  1:3525363/150‑1 (MQ=255)
                                                                                                                                 cgccgccgcCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTCCCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCAtggt                     >  2:2360076/1‑150 (MQ=255)
                                                                                                                                    cgccgccATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATggtggt                  <  2:4476612/150‑1 (MQ=255)
                                                                                                                                      ccgccATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAtggtg                                                                                   >  1:973058/1‑83 (MQ=255)
                                                                                                                                      ccgccATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTgc                <  2:3609614/150‑1 (MQ=255)
                                                                                                                                       cgccATGCTGGGCACCAGCATCCCCAAACTGTCCGCCACCGACAGCCTCGAAGCCGTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGGGGTgcg               >  2:2327313/1‑150 (MQ=255)
                                                                                                                                          gaTGCTGGTGCCCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGaa                                                                                >  1:4893985/2‑82 (MQ=255)
                                                                                                                                          gaTGCTGGTGCCCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGaa                                                                                <  2:4894016/81‑1 (MQ=255)
                                                                                                                                          gaTGCTGGTGCCCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGaa                                                                                <  2:4893985/81‑1 (MQ=255)
                                                                                                                                          gaTGCTGGTGCCCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCAGGGTGGaa                                                                                >  1:4894016/2‑82 (MQ=255)
                                                                                                                                          cATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCAc            >  2:3860080/1‑150 (MQ=255)
                                                                                                                                          cATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTGCGGCTCCTGTCGCATGTGCATGGTGGACATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCAc            >  2:171238/1‑150 (MQ=255)
                                                                                                                                             gCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGa         >  1:4350766/1‑150 (MQ=255)
                                                                                                                                              cTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGaa        >  1:401697/1‑150 (MQ=255)
                                                                                                                                              cTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGaa        >  1:3175871/1‑150 (MQ=255)
                                                                                                                                               tGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGaaa       >  2:257930/1‑150 (MQ=255)
                                                                                                                                               tGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGaaa       >  2:4493909/1‑150 (MQ=255)
                                                                                                                                               tGGGCACCAGCATCCCCAAACTGTGCACCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGaaa       >  1:980568/1‑150 (MQ=255)
                                                                                                                                                  gCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGAAAccc    >  1:127822/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATTTCTTCGACCCCGCTTCTGATAAAAGCAGCGACCTCGGCACCCCGGCCCGCGAAAGCGACGTGCAGGTCAGCCTGAGCATCGACGGCCGCGCCATCAGCGTGCCTGCTGGCACCTCGGTCATGCGCGCCGCCGCCATGCTGGGCACCAGCATCCCCAAACTGTGCGCCACCGACAGCCTCGAAGCCTTCGGCTCCTGCCGCATGTGCATGGTGGAAATCGAAGGCATGCGCGGCTACCCGGCATCCTGCACCACGCCAGTCAGCGAAGGCATGGTGGTGCGCACCGAAACCCCG  >  NC_002947/2492548‑2492845

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