Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,661,354 C→A 5.8% A41A (GCG→GCT ynfK ← putative dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,661,3540CA 5.8% 751.6 / 13.4 197A41A (GCG→GCTynfKputative dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base C (103/82);  new base A (10/2);  total (113/84)
Fisher's exact test for biased strand distribution p-value = 7.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

ATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCAC  >  CP009273/1661208‑1661501
                                                                                                                                                  |                                                                                                                                                   
  tgGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCt                                                                                                                                                  >  1:1840531/1‑148 (MQ=255)
      cAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGCGTGTCTCTTTGCTCCCCTTCGCTACgg                                                                                                                                              >  2:1852843/1‑148 (MQ=255)
      cAACTGTGCGCCACGCTACTTTCTTCTTCACTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACgg                                                                                                                                              >  2:3065245/1‑148 (MQ=255)
          tgtgCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGttt                                                                                                                                          <  2:2447656/148‑1 (MQ=255)
                ccACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTNTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATcc                                                                                                                                    >  1:1679334/1‑148 (MQ=255)
                      tACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAAc                                                                                                                              >  1:564505/1‑148 (MQ=255)
                      tACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAAc                                                                                                                              >  2:2790307/1‑148 (MQ=255)
                      tACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAAc                                                                                                                              >  2:1993819/1‑148 (MQ=255)
                        cTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCg                                                                                                                            <  1:2451851/148‑1 (MQ=255)
                         tttCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGt                                                                                                                           <  2:2547907/148‑1 (MQ=255)
                         tttCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGt                                                                                                                           <  2:2548426/148‑1 (MQ=255)
                         tttCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGt                                                                                                                           <  2:796057/148‑1 (MQ=255)
                          ttcttcttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGtt                                                                                                                          <  1:962912/148‑1 (MQ=255)
                           tcttcttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGttt                                                                                                                         <  1:2519187/148‑1 (MQ=255)
                            cttcttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGtttt                                                                                                                        >  2:3076064/1‑148 (MQ=255)
                            cttcttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGtttt                                                                                                                        <  1:3162284/148‑1 (MQ=255)
                               cttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTcc                                                                                                                     <  1:973427/148‑1 (MQ=255)
                                ttcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTccc                                                                                                                    <  1:80118/148‑1 (MQ=255)
                                 tcGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTccat                                                                                                                   <  2:863985/148‑3 (MQ=255)
                                    cTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTggg                                                                                                                <  1:2911030/148‑1 (MQ=255)
                                      tAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTggg                                                                                                                >  1:2973256/1‑146 (MQ=255)
                                      tAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTggg                                                                                                                <  2:2973256/146‑1 (MQ=255)
                                       aaCGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAgg                                                                                                             <  2:1994476/148‑1 (MQ=255)
                                          gcgATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTa                                                                                                          <  1:2330749/148‑1 (MQ=255)
                                          gcgATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTa                                                                                                          <  2:2395851/148‑1 (MQ=255)
                                           cgATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTaa                                                                                                         <  1:3039111/148‑1 (MQ=255)
                                             aTAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACg                                                                                                       >  2:2025211/1‑148 (MQ=255)
                                              tAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGc                                                                                                      >  1:2495719/1‑148 (MQ=255)
                                              tAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGc                                                                                                      >  2:1260273/1‑148 (MQ=255)
                                               aGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCt                                                                                                     >  1:2670761/1‑148 (MQ=255)
                                                 gATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTg                                                                                                   <  1:522379/148‑1 (MQ=255)
                                                      aCTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCa                                                                                              >  1:3278097/1‑148 (MQ=255)
                                                      aCTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCa                                                                                              >  1:531471/1‑148 (MQ=255)
                                                         gCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATg                                                                                           <  1:1142984/148‑1 (MQ=255)
                                                         gCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATg                                                                                           >  1:614842/1‑148 (MQ=255)
                                                             cATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGcgcg                                                                                       <  1:1260962/148‑1 (MQ=255)
                                                              aTAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCgg                                                                                      >  1:1365683/1‑148 (MQ=255)
                                                              aTAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCgg                                                                                      <  2:2703188/148‑1 (MQ=255)
                                                              aTAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCgg                                                                                      >  1:2571102/1‑148 (MQ=255)
                                                                aaGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGa                                                                                    >  2:50963/1‑148 (MQ=255)
                                                                aaGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGa                                                                                    >  1:1612581/1‑148 (MQ=255)
                                                                aaGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGa                                                                                    >  2:441292/1‑148 (MQ=255)
                                                                 aGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGaa                                                                                   >  2:2149939/1‑148 (MQ=255)
                                                                  ggCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGaaa                                                                                  >  1:756161/1‑148 (MQ=255)
                                                                  ggCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGaaa                                                                                  >  2:1341384/1‑148 (MQ=255)
                                                                  ggCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTAATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGaaa                                                                                  >  2:1342108/1‑148 (MQ=255)
                                                                   gCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAac                                                                                 <  1:2847454/148‑1 (MQ=255)
                                                                    cAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAacc                                                                                >  2:1343297/1‑148 (MQ=255)
                                                                         cGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGt                                                                           >  1:146404/1‑148 (MQ=255)
                                                                            tGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGtttt                                                                        >  2:2087654/1‑148 (MQ=255)
                                                                               ttGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTccc                                                                     >  1:2358586/1‑148 (MQ=255)
                                                                                tGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCt                                                                    >  1:1012705/1‑148 (MQ=255)
                                                                                  aaaCACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTAc                                                                  >  1:2785251/1‑148 (MQ=255)
                                                                                  aaaCACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTAc                                                                  >  2:1133212/1‑148 (MQ=255)
                                                                                   aaCACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACa                                                                 >  2:2601864/1‑148 (MQ=255)
                                                                                      acTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGaa                                                              <  2:155783/148‑1 (MQ=255)
                                                                                       ctctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAg                                                             <  1:2287852/148‑1 (MQ=255)
                                                                                       ctctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAg                                                             <  1:3195567/148‑1 (MQ=255)
                                                                                        tctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAgt                                                            >  1:483773/1‑148 (MQ=255)
                                                                                         ctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAgtg                                                           <  1:1496435/148‑1 (MQ=255)
                                                                                         ctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAgtg                                                           >  2:2195972/1‑148 (MQ=255)
                                                                                         ctGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAgtg                                                           >  2:2193650/1‑148 (MQ=255)
                                                                                           gCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTATCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTc                                                         <  2:285106/148‑1 (MQ=255)
                                                                                           gCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTc                                                         <  2:531445/148‑1 (MQ=255)
                                                                                           gCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTc                                                         >  2:2480390/1‑148 (MQ=255)
                                                                                            caccacCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCt                                                        <  2:1734740/145‑1 (MQ=255)
                                                                                               cacCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTa                                                     >  2:2342836/1‑148 (MQ=255)
                                                                                               cacCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTa                                                     >  2:1013888/1‑148 (MQ=255)
                                                                                                 ccAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTAcc                                                   >  2:3111564/1‑148 (MQ=255)
                                                                                                 ccAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTAcc                                                   >  2:1470675/1‑148 (MQ=255)
                                                                                                  cAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCg                                                  >  2:1219903/1‑148 (MQ=255)
                                                                                                     ggCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGta                                               <  1:2010052/148‑1 (MQ=255)
                                                                                                      gCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGtaa                                              <  1:369199/148‑1 (MQ=255)
                                                                                                        aTCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGtaata                                            >  1:3214331/1‑148 (MQ=255)
                                                                                                          cTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATaaa                                          <  1:1185358/148‑1 (MQ=255)
                                                                                                          cTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATaaa                                          >  1:3224970/1‑148 (MQ=255)
                                                                                                          cTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATaaa                                          <  2:841583/148‑1 (MQ=255)
                                                                                                          cTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATaaa                                          >  1:3223600/1‑148 (MQ=255)
                                                                                                            ttattaCGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGa                                        <  2:2913515/148‑1 (MQ=255)
                                                                                                            ttattaCGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGa                                        >  1:374126/1‑148 (MQ=255)
                                                                                                             tattaCGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGaa                                       >  1:2554291/1‑148 (MQ=255)
                                                                                                             tattaCGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGaa                                       >  1:2689806/1‑148 (MQ=255)
                                                                                                               ttaCGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAAc                                     >  2:369060/1‑148 (MQ=255)
                                                                                                                  cGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCt                                  >  2:1967079/1‑148 (MQ=255)
                                                                                                                    cAGCCCTTCGGGTGTCTCTTTGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTc                                <  1:441292/148‑1 (MQ=255)
                                                                                                                      gCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAg                              >  1:3302192/1‑148 (MQ=255)
                                                                                                                        ccTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCa                            >  1:1856843/1‑148 (MQ=255)
                                                                                                                        ccTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCa                            <  2:1575155/148‑1 (MQ=255)
                                                                                                                        ccTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCa                            >  1:168878/1‑148 (MQ=255)
                                                                                                                         cTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCAt                           >  1:3191760/1‑148 (MQ=255)
                                                                                                                         cTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCAt                           <  1:410909/148‑1 (MQ=255)
                                                                                                                               gtgtCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTAGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTaa                     >  2:630195/1‑148 (MQ=255)
                                                                                                                                 gtCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACt                   >  2:1373586/1‑148 (MQ=255)
                                                                                                                                 gtCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACt                   <  1:1993819/148‑1 (MQ=255)
                                                                                                                                    tctTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCa                <  1:2070292/148‑1 (MQ=255)
                                                                                                                                      tttGCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTCCAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCAcc              >  1:2179417/1‑148 (MQ=255)
                                                                                                                                         gCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGtt           <  2:742101/148‑1 (MQ=255)
                                                                                                                                         gCTCCCCTTAGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGtt           >  1:812259/1‑148 (MQ=255)
                                                                                                                                          cTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTa          <  2:1152055/148‑1 (MQ=255)
                                                                                                                                               cTTCGCTACGGGTTTATATCCCGCAACCGATTTTCCATGGGAGGCTAACGCTTGTAGCAAAGCGCGGGAAACCACCGTTTTAACTACAGAAGAGTCTGTACCTGTAATAAATAAACGCTTCAGCATCACTAACTCCACCGATATGCtt     >  1:2977473/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGATTAACTGCTTCATAAGGCAGTTCGATGGTTGAAACACTCTGCAACACCAGGGCATCTTTATTACGCAGCCCTTCGGGTGTCTCTTTGCTCCCCTTCGCTACGGGTTTATATCCCGCAACCGTTTTTCCCTGGGAGGCTAACGCTTGTAGCAATGCGCGGGAAACCACCGTTTTCCCTACAGAAGTGTCTGTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCAC  >  CP009273/1661208‑1661501

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: