Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,931,021 C→A 100% A123E (GCG→GAG)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,931,0210CA99.7% 985.7 / ‑4.4 293A123E (GCG→GAG) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (147/145);  total (147/146)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 296 total aligned reads displayed.

GCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGC  >  CP009273/3930874‑3931166
                                                                                                                                                   |                                                                                                                                                 
gCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgagc                                                                                                                                                 >  1:2588355/1‑150 (MQ=255)
  cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtg                                                                                                                                               <  1:2424518/150‑1 (MQ=255)
   agcagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                              <  1:4881881/150‑1 (MQ=255)
   agcagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                              >  2:1697486/1‑150 (MQ=255)
     cagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtggtg                                                                                                                                            >  2:595365/1‑150 (MQ=255)
      agcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtggtgg                                                                                                                                           >  2:771173/1‑150 (MQ=255)
         tCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAg                                                                                                                                        >  1:2511746/1‑150 (MQ=255)
          cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgc                                                                                                                                       <  1:1716433/150‑1 (MQ=255)
             cACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCa                                                                                                                                    <  1:4563468/150‑1 (MQ=255)
             cACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCa                                                                                                                                    >  2:4547716/1‑150 (MQ=255)
              aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGATCTGGTGGAAGCGCaa                                                                                                                                   >  1:3133072/1‑150 (MQ=255)
               cTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAAc                                                                                                                                  <  1:937230/150‑1 (MQ=255)
                 gATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGt                                                                                                                                >  1:3320959/1‑150 (MQ=255)
                   tAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                              <  1:379785/150‑1 (MQ=255)
                     aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGc                                                                                                                            <  1:3001052/150‑1 (MQ=255)
                       aTTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                          <  2:3399017/150‑1 (MQ=255)
                          gATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTAtt                                                                                                                       <  2:4984670/150‑1 (MQ=255)
                               tACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGTCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                                                                  <  1:2321399/150‑1 (MQ=255)
                                aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAc                                                                                                                 <  1:4687781/150‑1 (MQ=255)
                                aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAc                                                                                                                 >  2:3324955/1‑150 (MQ=255)
                                aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAc                                                                                                                 <  1:5024894/150‑1 (MQ=255)
                                 cTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAcg                                                                                                                >  2:3856620/1‑150 (MQ=255)
                                  tCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAcgc                                                                                                               >  1:1256631/1‑150 (MQ=255)
                                     ggTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAAAGCTGCCCTTTCCCCGGAGATGGTGGAAGCGCAACGTGAGCGTATTGACAACGCGTc                                                                                                            >  1:3821529/1‑150 (MQ=255)
                                          gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCat                                                                                                       <  2:3443257/150‑1 (MQ=255)
                                          gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCat                                                                                                       >  1:1798714/1‑150 (MQ=255)
                                            gTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatta                                                                                                     >  2:355670/1‑150 (MQ=255)
                                             tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCattat                                                                                                    >  1:4149596/1‑150 (MQ=255)
                                                 cAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                                                                >  1:2186604/1‑150 (MQ=255)
                                                  aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATg                                                                                               >  2:4307627/1‑150 (MQ=255)
                                                    cggCGACTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                                                             <  2:4838084/149‑1 (MQ=255)
                                                    aGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                                                             >  1:1692309/1‑150 (MQ=255)
                                                       cGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCt                                                                                          <  2:4471702/150‑1 (MQ=255)
                                                        gAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                                                         >  2:5061888/1‑150 (MQ=255)
                                                        gAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                                                         <  1:3083207/150‑1 (MQ=255)
                                                        gAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCAGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATTCAGCTg                                                                                         >  1:2913930/1‑150 (MQ=255)
                                                         aaTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                                                        <  1:595365/150‑1 (MQ=255)
                                                           tCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                      <  2:2288168/150‑1 (MQ=255)
                                                           tCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                      >  1:3901618/1‑150 (MQ=255)
                                                           tCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                      >  1:1483095/1‑150 (MQ=255)
                                                           tCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCAGCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGac                                                                                      >  2:123391/1‑149 (MQ=255)
                                                             aaCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    >  1:355671/1‑150 (MQ=255)
                                                             aaCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    >  1:721025/1‑150 (MQ=255)
                                                               cAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcac                                                                                  >  2:1625835/1‑150 (MQ=255)
                                                               cAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcac                                                                                  >  2:5051805/1‑150 (MQ=255)
                                                                 ggTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacca                                                                                <  2:2988645/150‑1 (MQ=255)
                                                                  gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcaccac                                                                               <  1:2335726/150‑1 (MQ=255)
                                                                   tgtgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACt                                                                              >  2:3159757/1‑150 (MQ=255)
                                                                   tgtgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACt                                                                              >  2:1853676/1‑150 (MQ=255)
                                                                    gtgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTc                                                                             >  2:1803986/1‑150 (MQ=255)
                                                                      ggCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGa                                                                           >  2:1904061/1‑150 (MQ=255)
                                                                       gcgcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTTATGCAGCTGGAATCACCACTCGaa                                                                          <  1:4259683/150‑1 (MQ=255)
                                                                         tcTGATTTTTGTTAATTGCGTAGGTGAGAATTTCATAGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGTTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                                                        <  2:2641268/149‑1 (MQ=255)
                                                                         gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                                                        >  1:3990620/1‑150 (MQ=255)
                                                                          cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                       <  2:4295436/150‑1 (MQ=255)
                                                                             aTTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                                    <  1:1113812/150‑1 (MQ=255)
                                                                             aTTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                                    <  1:2824971/150‑1 (MQ=255)
                                                                             aTTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                                    <  2:4584018/150‑1 (MQ=255)
                                                                              tttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGa                                                                   >  1:790788/1‑150 (MQ=255)
                                                                               ttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAt                                                                  <  2:3078852/150‑1 (MQ=255)
                                                                                tttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                                                                 <  1:3697994/150‑1 (MQ=255)
                                                                                tttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                                                                 >  2:4701717/1‑150 (MQ=255)
                                                                                  tGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGc                                                               <  2:2381458/150‑1 (MQ=255)
                                                                                    ttAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAg                                                             <  2:2511746/150‑1 (MQ=255)
                                                                                          gCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGa                                                       <  2:232797/150‑1 (MQ=255)
                                                                                                  gagaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGcc                                               >  1:2307855/1‑150 (MQ=255)
                                                                                                  gagaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGcc                                               >  1:1983776/1‑150 (MQ=255)
                                                                                                   agaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc                                              >  2:88862/1‑150 (MQ=255)
                                                                                                    gaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCa                                             <  1:2693083/150‑1 (MQ=255)
                                                                                                        gTCATCGGTATTCATGCCGGCGCTAATGCTGCCCATTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaa                                         >  1:1670485/1‑150 (MQ=255)
                                                                                                          cATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaaaa                                       >  1:1845169/1‑150 (MQ=255)
                                                                                                            tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATa                                     <  1:647563/150‑1 (MQ=255)
                                                                                                              ggTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAg                                   >  1:1918012/1‑150 (MQ=255)
                                                                                                              ggTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAg                                   >  2:3749411/1‑150 (MQ=255)
                                                                                                               gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGa                                  >  1:2105557/1‑150 (MQ=255)
                                                                                                                tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGAc                                 >  2:1134525/1‑150 (MQ=255)
                                                                                                                tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGAGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGAc                                 <  2:334940/150‑1 (MQ=255)
                                                                                                                   tCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTAt                              <  2:1349483/150‑1 (MQ=255)
                                                                                                                       gCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGtt                          <  1:1628786/150‑1 (MQ=255)
                                                                                                                        ccGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTg                         <  1:350755/150‑1 (MQ=255)
                                                                                                                          ggCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTgcg                       >  1:1190683/1‑150 (MQ=255)
                                                                                                                             gcTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt                    <  1:705125/150‑1 (MQ=255)
                                                                                                                              cTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                   <  2:2034646/150‑1 (MQ=255)
                                                                                                                              cTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                   <  2:4144643/150‑1 (MQ=255)
                                                                                                                              cTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                   >  2:1097800/1‑150 (MQ=255)
                                                                                                                              cTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                   >  2:1568035/1‑150 (MQ=255)
                                                                                                                               tAATTCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAAGGTGAGCGTATTGCCAAAGCGTCAGCATTTTTAATGCAGATGGAATCACCTCTCGAAAGTGTGATTGCAGCGGAGTAAATCGCCCTTCAAAATAAGACTATCGTTGCGCTTaa                  <  2:1936032/150‑1 (MQ=255)
                                                                                                                               tAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTaa                  >  1:806510/1‑150 (MQ=255)
                                                                                                                                aaTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAc                 >  1:1440633/1‑150 (MQ=255)
                                                                                                                                 aTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGGTGCGCTTAAcc                <  2:280387/150‑1 (MQ=255)
                                                                                                                                    ctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccgg             >  1:3184644/1‑150 (MQ=255)
                                                                                                                                      gcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggct           <  2:3307104/150‑1 (MQ=255)
                                                                                                                                      gcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggct           <  2:721792/150‑1 (MQ=255)
                                                                                                                                       cccTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggctc          <  2:3679372/150‑1 (MQ=255)
                                                                                                                                             ccccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGACTCACCACCCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCcggctccggctc    >  1:4614688/1‑150 (MQ=255)
                                                                                                                                             ccccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCcggctccggctc    >  1:216364/1‑150 (MQ=255)
                                                                                                                                              cccGGATCTGGTGGAATCGCAACGTGAGCGTATTGCAAACGAGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTcg   <  2:2190343/150‑1 (MQ=255)
                                                                                                                                              cccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTcg   <  1:448018/150‑1 (MQ=255)
                                                                                                                                              cccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTcg   >  1:1376634/1‑150 (MQ=255)
                                                                                                                                               ccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTcgc  >  2:3666378/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGC  >  CP009273/3930874‑3931166

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: