Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,519,134 C→A 6.4% L309L (CTG→CTT sgcC ← putative PTS system EIIC permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,519,1340CA 6.4% 631.4 / 14.7 170L309L (CTG→CTTsgcCputative PTS system EIIC permease component
Reads supporting (aligned to +/- strand):  ref base C (77/82);  new base A (4/7);  total (81/89)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 175 total aligned reads displayed.

GATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGG  >  CP009273/4518987‑4519281
                                                                                                                                                   |                                                                                                                                                   
gATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATc                                                                                                                                                     <  2:980168/148‑1 (MQ=255)
gATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATc                                                                                                                                                     <  2:1394285/148‑1 (MQ=255)
 aTGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCa                                                                                                                                                    <  2:1626885/148‑1 (MQ=255)
 aTGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCa                                                                                                                                                    >  1:3209679/1‑148 (MQ=255)
  tGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAg                                                                                                                                                   >  2:2147085/1‑148 (MQ=255)
   gACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGt                                                                                                                                                  >  2:1667806/1‑148 (MQ=255)
      cccGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGc                                                                                                                                               >  1:30585/1‑148 (MQ=255)
        cGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAg                                                                                                                                             <  1:433553/148‑1 (MQ=255)
          cTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGAc                                                                                                                                           >  1:3090474/1‑148 (MQ=255)
          cTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATAAGTAGCAGAc                                                                                                                                           <  1:2860490/148‑1 (MQ=255)
            tATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCg                                                                                                                                         >  1:739609/1‑148 (MQ=255)
             aTCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGa                                                                                                                                        <  2:2857053/148‑1 (MQ=255)
             aTCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGa                                                                                                                                        <  2:2138371/148‑1 (MQ=255)
                aaCGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAg                                                                                                                                     <  2:2910324/148‑1 (MQ=255)
                  cGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCa                                                                                                                                   <  1:808522/148‑1 (MQ=255)
                  cGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATAAGTAGCAGACCGACAGCa                                                                                                                                   >  1:2825981/1‑148 (MQ=255)
                       tGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAAtggt                                                                                                                              >  2:222172/1‑148 (MQ=255)
                       tGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAAtggt                                                                                                                              >  2:222167/1‑148 (MQ=255)
                        gACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATCAACATAATGGGTATCAGTAGCAGACCGACAGCAAtggtg                                                                                                                             >  2:2102165/1‑148 (MQ=255)
                         aCCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAAtggtgg                                                                                                                            <  1:3126816/148‑1 (MQ=255)
                            aGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTg                                                                                                                         >  1:1100482/1‑148 (MQ=255)
                              gTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGa                                                                                                                       <  2:739703/148‑1 (MQ=255)
                              gTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGa                                                                                                                       >  1:1032143/1‑148 (MQ=255)
                              gTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGACAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGa                                                                                                                       >  1:2143806/1‑148 (MQ=255)
                                       gTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGc                                                                                                              >  1:434884/1‑148 (MQ=255)
                                          gATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATAAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTc                                                                                                           >  2:310462/1‑148 (MQ=255)
                                              aCCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCAcag                                                                                                       >  2:2393179/1‑148 (MQ=255)
                                               ccGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCAcagc                                                                                                      >  2:691303/1‑148 (MQ=255)
                                                cGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCAcagca                                                                                                     <  2:1503808/148‑1 (MQ=255)
                                                cGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCAcagca                                                                                                     <  2:1503453/148‑1 (MQ=255)
                                                     ccATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATc                                                                                                <  2:844113/148‑1 (MQ=255)
                                                        tACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAg                                                                                             <  2:2012598/148‑1 (MQ=255)
                                                         aCAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGa                                                                                            >  1:1081231/1‑148 (MQ=255)
                                                         aCAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGa                                                                                            <  1:2179175/148‑1 (MQ=255)
                                                         aCAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGa                                                                                            <  1:189814/148‑1 (MQ=255)
                                                             aaaaaGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATTCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATTCCCCAGCGTCACAGCAGTATCcagacaga                                                                                        >  1:898167/3‑148 (MQ=255)
                                                             aTAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGa                                                                                        >  2:3266086/1‑148 (MQ=255)
                                                             aTAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGa                                                                                        >  2:3262884/1‑148 (MQ=255)
                                                                 aGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATNATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGAtata                                                                                    >  1:1680547/1‑148 (MQ=255)
                                                                 aGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGAtata                                                                                    <  2:185444/148‑1 (MQ=255)
                                                                 aGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGAtata                                                                                    <  1:2223014/148‑1 (MQ=255)
                                                                   aaCGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATaca                                                                                  >  2:343499/1‑148 (MQ=255)
                                                                    aCGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATacac                                                                                 <  1:2957576/148‑1 (MQ=255)
                                                                        tGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCa                                                                             <  1:1196/148‑1 (MQ=255)
                                                                        tGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCa                                                                             <  2:2188186/148‑1 (MQ=255)
                                                                        tGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATAAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCa                                                                             <  2:3255668/148‑1 (MQ=255)
                                                                            aCCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTc                                                                         <  1:222094/148‑1 (MQ=255)
                                                                            aCCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTc                                                                         <  1:222172/148‑1 (MQ=255)
                                                                                gCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCttt                                                                     <  2:1320830/148‑1 (MQ=255)
                                                                                 cAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTg                                                                    >  1:148274/1‑148 (MQ=255)
                                                                                 cAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTg                                                                    <  2:1301204/148‑1 (MQ=255)
                                                                                 cAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTg                                                                    <  2:1778260/148‑1 (MQ=255)
                                                                                    gTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGaaa                                                                 >  1:1739794/1‑148 (MQ=255)
                                                                                    gTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGaaa                                                                 >  1:2998257/1‑148 (MQ=255)
                                                                                        gCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATAtt                                                             <  1:861975/148‑1 (MQ=255)
                                                                                         cGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATAttt                                                            <  2:687475/148‑1 (MQ=255)
                                                                                          gAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTc                                                           >  1:1168497/1‑148 (MQ=255)
                                                                                           aGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGTGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGACCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCt                                                          >  1:2731512/1‑148 (MQ=255)
                                                                                            ggggTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTg                                                         >  1:2052960/1‑148 (MQ=255)
                                                                                             gggTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTgg                                                        <  2:2961565/148‑1 (MQ=255)
                                                                                               gTACCACTTTGGTGCCCGGTAAAATGCGGGCCAGAATTAACATAATGGGTATCCGTCGCAGACCGACAGCACTGCTGGTTGGATGCCCCAGCTTCAAAGCAGTTTCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGgaa                                                      <  2:836606/148‑1 (MQ=255)
                                                                                                 aaCACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGgaaga                                                    <  2:1587039/148‑1 (MQ=255)
                                                                                                 aaCACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGgaaga                                                    <  2:2537871/148‑1 (MQ=255)
                                                                                                 aaCACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGgaaga                                                    <  1:748679/148‑1 (MQ=255)
                                                                                                  acacTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGgaagaa                                                   >  2:641165/1‑148 (MQ=255)
                                                                                                    acTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGTAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACt                                                 <  1:1392045/148‑1 (MQ=255)
                                                                                                     cTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACtt                                                >  1:2636286/1‑148 (MQ=255)
                                                                                                     cTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACtt                                                >  2:905247/1‑148 (MQ=255)
                                                                                                      tttGTTGCCCGGTAAAATGCTGGCCAGAATTTACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACttt                                               <  1:2092633/148‑1 (MQ=255)
                                                                                                      tttGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACttt                                               <  1:2092649/148‑1 (MQ=255)
                                                                                                      tttGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACttt                                               <  2:2566057/148‑1 (MQ=255)
                                                                                                       ttgttgCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTc                                              <  2:2034790/148‑1 (MQ=255)
                                                                                                       ttgttgCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTc                                              <  2:1699944/148‑1 (MQ=255)
                                                                                                        tgttgCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCt                                             <  1:2324471/148‑1 (MQ=255)
                                                                                                         gttgCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCtt                                            <  2:217796/148‑1 (MQ=255)
                                                                                                         gttgCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCtt                                            <  1:2393179/148‑1 (MQ=255)
                                                                                                             cccGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCAc                                        <  2:1119688/148‑1 (MQ=255)
                                                                                                                  tAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCg                                   >  1:808910/1‑148 (MQ=255)
                                                                                                                    aaaTGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGga                                 >  1:2232119/1‑148 (MQ=255)
                                                                                                                    aaaTGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGga                                 >  1:2232009/1‑148 (MQ=255)
                                                                                                                       tGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGAt                              <  2:1032143/148‑1 (MQ=255)
                                                                                                                        gCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATg                             >  2:105228/1‑148 (MQ=255)
                                                                                                                             ccAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGgcag                        <  1:759141/148‑1 (MQ=255)
                                                                                                                                gAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGcagcag                     >  2:2096848/1‑148 (MQ=255)
                                                                                                                                 aaTTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAgg                    >  1:2918572/1‑148 (MQ=255)
                                                                                                                                   ttAACATACTGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGcc                  <  2:30585/148‑1 (MQ=255)
                                                                                                                                   ttAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGcc                  <  1:2931667/148‑1 (MQ=255)
                                                                                                                                   ttAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGcc                  >  1:2778389/1‑148 (MQ=255)
                                                                                                                                       cATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTca              >  2:1473979/1‑148 (MQ=255)
                                                                                                                                             gggTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGTCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATc        <  2:2307235/148‑1 (MQ=255)
                                                                                                                                              ggTATCAGTAGCAGACCGACTGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCa       <  2:2731512/148‑1 (MQ=255)
                                                                                                                                              ggTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCa       <  1:3266086/148‑1 (MQ=255)
                                                                                                                                              ggTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCa       <  1:3262900/148‑1 (MQ=255)
                                                                                                                                              ggTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCa       >  1:2481182/1‑148 (MQ=255)
                                                                                                                                              ggTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCa       >  1:2481208/1‑148 (MQ=255)
                                                                                                                                               gTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAg      >  1:2580993/1‑148 (MQ=255)
                                                                                                                                               gTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAg      <  1:2101850/148‑1 (MQ=255)
                                                                                                                                                tATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAgg     <  1:1183053/148‑1 (MQ=255)
                                                                                                                                                 aTCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGc    >  2:1180450/1‑148 (MQ=255)
                                                                                                                                                   cAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCgg  >  1:338793/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTAAAATGCTGGCCAGAATTAACATAATGGGTATCAGTAGCAGACCGACAGCAATGGTGGTTGGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGG  >  CP009273/4518987‑4519281

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: