Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,234 C→T 100% E27K (GAA→AAA)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,2340CT98.8% 1154.2 / 1.3 344E27K (GAA→AAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (162/178);  total (165/179)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 346 total aligned reads displayed.

TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACT  >  CP009273/647087‑647380
                                                                                                                                                   |                                                                                                                                                  
   gcgcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCa                                                                                                                                               <  2:2289074/150‑1 (MQ=255)
      ccAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGa                                                                                                                                            <  1:3653557/150‑1 (MQ=255)
      ccAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGa                                                                                                                                            >  1:1974626/1‑150 (MQ=255)
       cAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGaa                                                                                                                                           >  2:2703802/1‑150 (MQ=255)
        aGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAt                                                                                                                                          >  1:5076123/1‑150 (MQ=255)
        aGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAt                                                                                                                                          >  2:1672048/1‑150 (MQ=255)
        aGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAt                                                                                                                                          <  2:4393810/150‑1 (MQ=255)
         gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                         <  2:3934292/150‑1 (MQ=255)
         gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                         >  1:2357870/1‑150 (MQ=255)
          tCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTg                                                                                                                                        <  2:3037524/150‑1 (MQ=255)
           ccTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTgg                                                                                                                                       <  1:1008560/150‑1 (MQ=255)
            cTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTggg                                                                                                                                      <  1:4533356/150‑1 (MQ=255)
              tCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCg                                                                                                                                    <  2:343738/150‑1 (MQ=255)
                  cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAttt                                                                                                                                <  2:2628696/150‑1 (MQ=255)
                   gcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTc                                                                                                                               <  1:1212451/150‑1 (MQ=255)
                    cggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTcc                                                                                                                              >  2:205523/1‑150 (MQ=255)
                    cggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTcc                                                                                                                              <  1:2405083/150‑1 (MQ=255)
                    cggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTcc                                                                                                                              <  1:5234425/150‑1 (MQ=255)
                      gacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGc                                                                                                                            >  2:4567522/1‑150 (MQ=255)
                      gacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGc                                                                                                                            >  2:3340007/1‑150 (MQ=255)
                        cggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCa                                                                                                                          >  2:1799823/1‑150 (MQ=255)
                         ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                         <  1:1467018/150‑1 (MQ=255)
                         ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                         <  2:1506241/150‑1 (MQ=255)
                          gaTTCACTGATAGCAACGCATTTAAAAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAta                                                                                                                        >  2:3533108/1‑149 (MQ=255)
                           aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttt                                                                                                                       >  2:473268/1‑150 (MQ=255)
                           aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttt                                                                                                                       >  2:5000202/1‑150 (MQ=255)
                             tCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttttt                                                                                                                     >  2:4357678/1‑150 (MQ=255)
                              cACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                    <  1:1401347/150‑1 (MQ=255)
                                cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTa                                                                                                                  >  1:4319115/1‑150 (MQ=255)
                                    tAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                              <  1:4493486/150‑1 (MQ=255)
                                     aGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTc                                                                                                             >  1:4903614/1‑150 (MQ=255)
                                     aGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTc                                                                                                             >  2:1352812/1‑150 (MQ=255)
                                         aCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGaa                                                                                                         <  1:1915278/150‑1 (MQ=255)
                                          cGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAAc                                                                                                        <  1:2022814/150‑1 (MQ=255)
                                           gCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                       <  2:2002280/150‑1 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                     >  2:4062516/1‑150 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                     >  2:67246/1‑150 (MQ=255)
                                               ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                   <  2:2550739/150‑1 (MQ=255)
                                                 aataatATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAc                                                                                                 >  2:779764/1‑150 (MQ=255)
                                                   taataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcgc                                                                                               >  2:125168/1‑150 (MQ=255)
                                                       atTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGt                                                                                           >  2:144331/1‑150 (MQ=255)
                                                         ttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGa                                                                                         >  1:1586572/1‑150 (MQ=255)
                                                             ggCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAat                                                                                     <  2:5140824/150‑1 (MQ=255)
                                                                aaTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATc                                                                                  <  1:1017487/150‑1 (MQ=255)
                                                                   tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                               <  1:335006/150‑1 (MQ=255)
                                                                   tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                               <  2:1561771/150‑1 (MQ=255)
                                                                     cACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGc                                                                             >  2:1827044/1‑150 (MQ=255)
                                                                      aCCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                            >  2:1576300/1‑150 (MQ=255)
                                                                           aCGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGTAAATATCATTGCCGAACa                                                                       <  1:4737597/150‑1 (MQ=255)
                                                                            cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                      <  2:1398870/150‑1 (MQ=255)
                                                                            cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                      <  1:3117842/150‑1 (MQ=255)
                                                                              cgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                                    >  2:4935922/1‑150 (MQ=255)
                                                                               gATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                                   >  1:2886669/1‑150 (MQ=255)
                                                                                 tAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACAtaataa                                                                 <  1:921389/150‑1 (MQ=255)
                                                                                  aaGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                                >  1:4008677/1‑150 (MQ=255)
                                                                                   aGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaat                                                               <  1:1807697/150‑1 (MQ=255)
                                                                                     cTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAg                                                             >  1:2662091/1‑150 (MQ=255)
                                                                                       tttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTa                                                           <  1:2703802/150‑1 (MQ=255)
                                                                                        ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                          <  1:4659415/150‑1 (MQ=255)
                                                                                        ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                          <  2:2435170/150‑1 (MQ=255)
                                                                                        ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                          <  1:4069361/150‑1 (MQ=255)
                                                                                         tCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATc                                                         <  2:1728183/150‑1 (MQ=255)
                                                                                           aTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                       >  2:3666709/1‑150 (MQ=255)
                                                                                           aTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                       >  2:5261719/1‑150 (MQ=255)
                                                                                               cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                   >  1:1666307/1‑150 (MQ=255)
                                                                                                 ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGt                                                 >  1:1520781/1‑150 (MQ=255)
                                                                                                  gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                <  2:56522/150‑1 (MQ=255)
                                                                                                      cGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtg                                            >  1:4588490/1‑150 (MQ=255)
                                                                                                      cGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtg                                            <  2:4665267/150‑1 (MQ=255)
                                                                                                        gTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGTAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgtt                                          >  2:2965402/1‑150 (MQ=255)
                                                                                                         taataaATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                         <  2:3961830/150‑1 (MQ=255)
                                                                                                          aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                        >  2:3769213/1‑150 (MQ=255)
                                                                                                          aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                        >  1:3910300/1‑150 (MQ=255)
                                                                                                           ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATgttgtt                                       <  1:736519/150‑1 (MQ=255)
                                                                                                             aaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTAt                                     >  2:4210874/1‑150 (MQ=255)
                                                                                                              tataCTTAGACTGTGGTTTCAACGATCAAATCATTTTTATGCAGGAATTGGGCGATTTCAGCCATTTTTTTATTTTCTTAACTTTTTACGCTGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                    <  2:3121948/149‑1 (MQ=255)
                                                                                                              aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                    >  1:2529156/1‑150 (MQ=255)
                                                                                                              aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                    <  2:4015964/150‑1 (MQ=255)
                                                                                                              aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                    >  1:3737155/1‑150 (MQ=255)
                                                                                                                 aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                 <  2:967778/150‑1 (MQ=255)
                                                                                                                  cTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAtt                                <  1:3513643/150‑1 (MQ=255)
                                                                                                                      gACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCa                            >  1:4437216/1‑150 (MQ=255)
                                                                                                                        cTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAAc                          <  1:57726/150‑1 (MQ=255)
                                                                                                                         tgtgGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                         >  1:2854756/1‑150 (MQ=255)
                                                                                                                         tgatgtgtAAACGCTCAAATCATATATAAGCAGGAATTGGGCGATTACCGCCATTTTTTTATTTTATGAACGTTTTACGAGGGTGAAAATATCATTGCCTAACATAATAAAAAGTATCAAGAAGGTGAATGTTCTTATCGATTTTCAAct                         >  2:4183670/5‑149 (MQ=255)
                                                                                                                                tCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttg                  <  1:2915323/150‑1 (MQ=255)
                                                                                                                                 cAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAATTGCAACGAAtgttgt                 >  1:5161776/1‑150 (MQ=255)
                                                                                                                                  aaCGCTCAAATCATTTTTATGCAGGAATTGGGCGATCTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATgttgtt                >  2:3733362/1‑150 (MQ=255)
                                                                                                                                  aaCGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATgttgtt                >  2:2520768/1‑150 (MQ=255)
                                                                                                                                    cGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa              <  1:1425534/150‑1 (MQ=255)
                                                                                                                                       tCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg           >  1:2574378/1‑150 (MQ=255)
                                                                                                                                        cAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGt          <  2:2357870/150‑1 (MQ=255)
                                                                                                                                         ataTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt         <  1:3195639/148‑1 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt         >  1:1839259/1‑150 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt         >  2:1579598/1‑150 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt         <  1:3265890/150‑1 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt         >  2:3186594/1‑150 (MQ=255)
                                                                                                                                          aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTg        <  2:217639/150‑1 (MQ=255)
                                                                                                                                          aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTg        <  2:3621406/150‑1 (MQ=255)
                                                                                                                                               ttttTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAAc   <  1:3137359/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACT  >  CP009273/647087‑647380

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: ‑