Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,073,741 A→T 13.8% L283Q (CTG→CAG)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,073,7410AT13.8% 611.9 / 75.8 217L283Q (CTG→CAG) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (80/107);  new base T (17/13);  total (97/120)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 222 total aligned reads displayed.

ATGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073594‑4073887
                                                                                                                                                   |                                                                                                                                                  
aTGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCa                                                                                                                                                    >  2:2270803/1‑148 (MQ=255)
   gtgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagag                                                                                                                                                 <  2:232783/148‑1 (MQ=255)
   gtgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagag                                                                                                                                                 <  2:2594687/148‑1 (MQ=255)
    tgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagaga                                                                                                                                                <  1:1256661/148‑1 (MQ=255)
         tcACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGa                                                                                                                                           >  2:505506/1‑148 (MQ=255)
         tcACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGa                                                                                                                                           <  2:1642453/148‑1 (MQ=255)
          cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAg                                                                                                                                          >  1:1230994/1‑148 (MQ=255)
           acacCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAgg                                                                                                                                         >  2:782170/1‑148 (MQ=255)
            cacCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGt                                                                                                                                        >  2:308176/1‑148 (MQ=255)
            cacCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGt                                                                                                                                        >  2:2235696/1‑148 (MQ=255)
              cccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCa                                                                                                                                      <  2:1829816/148‑1 (MQ=255)
               ccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCAt                                                                                                                                     <  2:1660501/148‑1 (MQ=255)
               ccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCAt                                                                                                                                     >  1:1046644/1‑148 (MQ=255)
               ccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCAt                                                                                                                                     <  1:783834/148‑1 (MQ=255)
                cTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATc                                                                                                                                    <  1:1988108/148‑1 (MQ=255)
                  aTGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTg                                                                                                                                  <  1:1120220/148‑1 (MQ=255)
                     aGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTa                                                                                                                               <  1:2129265/148‑1 (MQ=255)
                     aGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTa                                                                                                                               <  1:2129531/148‑1 (MQ=255)
                          tactCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGACACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCt                                                                                                                          <  1:1409495/145‑1 (MQ=255)
                              ccGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGa                                                                                                                      <  2:1990297/148‑1 (MQ=255)
                                 ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   >  1:1828777/1‑148 (MQ=255)
                                 ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   >  2:625492/1‑148 (MQ=255)
                                     aCCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTgg                                                                                                               >  1:1039873/1‑148 (MQ=255)
                                      ccTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGa                                                                                                              <  1:2965906/148‑1 (MQ=255)
                                      ccTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGa                                                                                                              <  1:2801576/148‑1 (MQ=255)
                                        tGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATc                                                                                                            >  1:1345606/1‑148 (MQ=255)
                                         gcgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCt                                                                                                           <  1:2107585/148‑1 (MQ=255)
                                            ccACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtttt                                                                                                        <  1:2927301/148‑1 (MQ=255)
                                            ccACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCAGAGCCTTCTTGCTCCAGTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtttt                                                                                                        <  1:3021449/148‑1 (MQ=255)
                                             cACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCttttt                                                                                                       <  2:3592070/148‑1 (MQ=255)
                                             cACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCttttt                                                                                                       <  1:1368785/148‑1 (MQ=255)
                                              aCCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTc                                                                                                      <  1:1647119/148‑1 (MQ=255)
                                              aCCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTc                                                                                                      <  2:1439602/148‑1 (MQ=255)
                                               cccTGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  2:69114/148‑1 (MQ=255)
                                               cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  2:1585099/148‑1 (MQ=255)
                                                ccGTAAAACGTAAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCt                                                                                                    <  2:153159/148‑1 (MQ=255)
                                                ccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCt                                                                                                    >  2:775682/1‑148 (MQ=255)
                                                ccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCt                                                                                                    >  1:577419/1‑148 (MQ=255)
                                                 cGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTg                                                                                                   >  1:726332/1‑148 (MQ=255)
                                                    aaaaCGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATa                                                                                                >  2:1870132/1‑148 (MQ=255)
                                                     aaaCGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATag                                                                                               <  1:2782116/148‑1 (MQ=255)
                                                     aaaCGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATag                                                                                               <  2:2821664/148‑1 (MQ=255)
                                                         gTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGaa                                                                                           <  1:1321174/148‑1 (MQ=255)
                                                           cAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAtt                                                                                         <  2:91622/148‑1 (MQ=255)
                                                           cAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAGTTTTGGATCTTTTTCCTGATAGAGAAtt                                                                                         >  1:2493577/1‑148 (MQ=255)
                                                              agaTTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTc                                                                                      <  2:3754559/148‑1 (MQ=255)
                                                                aTTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCTAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCGGGATAGAGAATTTTCAg                                                                                    >  1:1332712/1‑148 (MQ=255)
                                                                 ttACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAgg                                                                                   <  2:1334169/148‑1 (MQ=255)
                                                                  tACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGt                                                                                  >  1:2363430/1‑148 (MQ=255)
                                                                  tACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGt                                                                                  >  2:2855877/1‑148 (MQ=255)
                                                                   aCTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTa                                                                                 >  1:1519189/1‑148 (MQ=255)
                                                                   aCTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTa                                                                                 <  2:3304951/148‑1 (MQ=255)
                                                                   aCTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTa                                                                                 <  2:1078784/148‑1 (MQ=255)
                                                                    ctctCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTaa                                                                                >  2:3407116/1‑148 (MQ=255)
                                                                    ctctCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTaa                                                                                >  2:3407044/1‑148 (MQ=255)
                                                                     tctcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAg                                                                               >  1:11677/1‑148 (MQ=255)
                                                                       tcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAgtg                                                                             <  1:3105095/148‑1 (MQ=255)
                                                                        ccTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAgtgt                                                                            <  1:3704655/148‑1 (MQ=255)
                                                                             aCCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAg                                                                       >  2:129272/1‑148 (MQ=255)
                                                                             aCCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAg                                                                       >  1:381887/1‑148 (MQ=255)
                                                                             aCCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAg                                                                       >  2:2058491/1‑148 (MQ=255)
                                                                                cGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCa                                                                    <  1:1187252/148‑1 (MQ=255)
                                                                                cGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCa                                                                    <  1:1007741/148‑1 (MQ=255)
                                                                                 ggAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCAt                                                                   <  1:1191282/148‑1 (MQ=255)
                                                                                        ggAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCt                                                            >  2:3739330/1‑148 (MQ=255)
                                                                                        ggAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCt                                                            <  1:3192211/148‑1 (MQ=255)
                                                                                        ggAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCt                                                            <  1:46499/148‑1 (MQ=255)
                                                                                           aCGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTc                                                         >  1:2281423/1‑148 (MQ=255)
                                                                                           aCGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTc                                                         >  2:2097550/1‑148 (MQ=255)
                                                                                            cGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCg                                                        >  1:2453854/1‑148 (MQ=255)
                                                                                              ggTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGc                                                      <  1:1760104/148‑1 (MQ=255)
                                                                                                  gATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTa                                                  <  1:2854819/148‑1 (MQ=255)
                                                                                                    tGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAAt                                                <  2:1363217/148‑1 (MQ=255)
                                                                                                          tGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGc                                          >  1:540440/1‑148 (MQ=255)
                                                                                                             gggCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTg                                       <  2:1486712/148‑1 (MQ=255)
                                                                                                             gggCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTg                                       <  2:2978391/148‑1 (MQ=255)
                                                                                                              ggCTTAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGt                                      <  2:2804178/148‑1 (MQ=255)
                                                                                                              ggCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGt                                      <  2:2843897/148‑1 (MQ=255)
                                                                                                               gCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGtt                                     <  2:481098/148‑1 (MQ=255)
                                                                                                                  tAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGt                                  <  1:2518010/148‑1 (MQ=255)
                                                                                                                    gCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCa                                >  2:3815869/1‑148 (MQ=255)
                                                                                                                    gCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCa                                >  1:172961/1‑148 (MQ=255)
                                                                                                                     ccTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCAt                               >  1:2820552/1‑148 (MQ=255)
                                                                                                                        tcttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCa                            >  1:2314569/1‑148 (MQ=255)
                                                                                                                         cttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAg                           >  2:2954651/1‑148 (MQ=255)
                                                                                                                          ttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTACGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGc                          >  1:894664/1‑148 (MQ=255)
                                                                                                                               ccATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCa                     >  2:2357162/1‑148 (MQ=255)
                                                                                                                                  ttCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  <  1:3623320/148‑1 (MQ=255)
                                                                                                                                  ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  <  2:2021502/148‑1 (MQ=255)
                                                                                                                                  ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  <  2:1259400/148‑1 (MQ=255)
                                                                                                                                  ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCATTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  <  2:2740979/148‑1 (MQ=255)
                                                                                                                                   tCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                 <  2:3647004/148‑1 (MQ=255)
                                                                                                                                     gTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATg               <  1:2237521/148‑1 (MQ=255)
                                                                                                                                     gTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATg               <  2:3038663/148‑1 (MQ=255)
                                                                                                                                       gCGTCGAGTATCAGAGAGGAGAGGTCATCTGCTACGGCTTCGTTTTTTGGATCTTTTTCCTTCTAGAGAATTTTAAGGTAAGTTTCGCAGTAATCGCAGCTTTAGGCTTTACTCGCGGCCTGTTCGTCATACAGCGACCAGTAATGCa             <  2:3254313/148‑1 (MQ=255)
                                                                                                                                         gTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGt           <  2:3494850/148‑1 (MQ=255)
                                                                                                                                           ccAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGttt         >  2:3764382/1‑148 (MQ=255)
                                                                                                                                           ccAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGttt         >  2:2453913/1‑148 (MQ=255)
                                                                                                                                              gTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGcc      >  1:2324672/1‑148 (MQ=255)
                                                                                                                                                  cAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGctc  >  1:2593778/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073594‑4073887

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: