Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,956,165 C→A 44.4% E475D (GAG→GAT gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,956,1650CA44.4% 248.8 / inf 286E475D (GAG→GATgppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (80/79);  new base A (67/59);  total (147/139)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 297 total aligned reads displayed.

GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGT  >  CP009273/3956017‑3956312
                                                                                                                                                    |                                                                                                                                                   
gCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAAtct                                                                                                                                                    <  1:1612760/150‑1 (MQ=255)
  gAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATctct                                                                                                                                                  <  2:4002131/150‑1 (MQ=255)
   aaTCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCtt                                                                                                                                                 <  1:1695387/150‑1 (MQ=255)
    aTCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCttt                                                                                                                                                <  2:1460235/150‑1 (MQ=255)
     tCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTa                                                                                                                                               <  1:824513/150‑1 (MQ=255)
      cTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTAc                                                                                                                                              >  2:4093071/1‑150 (MQ=255)
       tGAAAGCGATGATGGCGGCAAAACGAGGGACATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTAcc                                                                                                                                             >  1:1144025/1‑150 (MQ=255)
       tGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTAcc                                                                                                                                             <  2:1438762/150‑1 (MQ=255)
       tGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTAcc                                                                                                                                             <  2:2170766/150‑1 (MQ=255)
           aGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGc                                                                                                                                         >  1:1239040/1‑150 (MQ=255)
            gCGATGATGGCGGCAAACCGAGGGAAATAATCAAGCCCGGACAAGATGCGTCAGCATCGCATCCGGCACTTACTCAATAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCg                                                                                                                                        <  2:2377526/150‑1 (MQ=255)
              gatgatGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGa                                                                                                                                      <  2:102668/150‑1 (MQ=255)
                tgatgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATg                                                                                                                                    >  2:1114051/1‑150 (MQ=255)
                  atgGCGGCAAAACGATGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGtt                                                                                                                                  <  1:2785443/150‑1 (MQ=255)
                    ggcggcAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTgg                                                                                                                                <  2:1748416/150‑1 (MQ=255)
                    ggcggcAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTgg                                                                                                                                <  1:2802097/150‑1 (MQ=255)
                     gcggcAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTggg                                                                                                                               <  1:4595483/150‑1 (MQ=255)
                       ggcACAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGtt                                                                                                                             >  2:2447953/1‑150 (MQ=255)
                           aaaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGcc                                                                                                                         >  1:772857/1‑150 (MQ=255)
                            aaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCa                                                                                                                        >  2:4183418/1‑150 (MQ=255)
                             tcGAGGGATATAAGCATGTCCGGAAATGATGCGTAAGCTTCGCATCCGGCACTTACTCATTAATGCTCTTCCAGCGGCCTGTGGACATATCTCTGCCACTGGCTTTCCTGAGCAATAATTTCTTTTCCCCGCGGATGTTGGGTTAGCCaa                                                                                                                       <  2:261704/149‑1 (MQ=255)
                             aCGAGTGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCaa                                                                                                                       >  1:1868386/1‑150 (MQ=255)
                              cGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAAc                                                                                                                      >  2:4621342/1‑150 (MQ=255)
                               gAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAAcc                                                                                                                     <  2:720525/150‑1 (MQ=255)
                                  ggAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTg                                                                                                                  <  1:3867533/150‑1 (MQ=255)
                                  ggAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTg                                                                                                                  <  1:679585/150‑1 (MQ=255)
                                    aaaTAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCg                                                                                                                >  1:3375554/1‑150 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                            <  1:4269953/150‑1 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                            <  2:1348187/150‑1 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                            <  1:3994136/150‑1 (MQ=255)
                                          tCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCg                                                                                                          >  1:502553/1‑150 (MQ=255)
                                           cacggCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACGTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGt                                                                                                         <  2:4049659/147‑1 (MQ=255)
                                            aaGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTc                                                                                                        <  1:3924978/150‑1 (MQ=255)
                                               gCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAg                                                                                                     <  2:4552666/150‑1 (MQ=255)
                                                 cGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGt                                                                                                   >  2:2322072/1‑150 (MQ=255)
                                                  ggAAAAGATGCGTCAGCATCGCATGCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGTACATAGCTCTGCCTCTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTTGGTTAGCCAACCTTGCGGAAGCGTCAAGGTc                                                                                                  <  1:3561805/150‑1 (MQ=255)
                                                  ggAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTc                                                                                                  >  1:4304887/1‑150 (MQ=255)
                                                     aaaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAAc                                                                                               >  1:2428283/1‑150 (MQ=255)
                                                     aaaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAAc                                                                                               <  1:4087068/150‑1 (MQ=255)
                                                     aaaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAAc                                                                                               <  1:2654774/150‑1 (MQ=255)
                                                      aaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACa                                                                                              >  1:2059551/1‑150 (MQ=255)
                                                            cGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAt                                                                                        >  2:3367125/1‑150 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                       >  2:753128/1‑150 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                       <  1:884731/150‑1 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                       >  2:580177/1‑150 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                       >  1:1283896/1‑150 (MQ=255)
                                                               cAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGt                                                                                     <  1:3032163/150‑1 (MQ=255)
                                                                     cgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCt                                                                               <  1:4054932/150‑1 (MQ=255)
                                                                       catcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGt                                                                             >  1:1539324/1‑150 (MQ=255)
                                                                         tcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTaa                                                                           >  2:2210790/1‑150 (MQ=255)
                                                                                cTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtc                                                                    >  2:3647854/1‑150 (MQ=255)
                                                                                       aTTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCa                                                             <  2:2608690/150‑1 (MQ=255)
                                                                                        ttAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCAc                                                            >  1:1306300/1‑150 (MQ=255)
                                                                                         tAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACg                                                           <  1:450472/150‑1 (MQ=255)
                                                                                               aCTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAAAATCTTTACCCACCGGATGTTGGGTTAGCCAACCTTGCGGACGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTAATCTCTGGCACGAGAtag                                                     >  2:2596370/1‑148 (MQ=255)
                                                                                                 ttCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTc                                                   >  1:3621225/1‑150 (MQ=255)
                                                                                                 ttCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTc                                                   >  1:3698188/1‑150 (MQ=255)
                                                                                                  tCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCa                                                  <  2:4360951/150‑1 (MQ=255)
                                                                                                    cAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGAGGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTAATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACg                                                >  2:1992743/1‑150 (MQ=255)
                                                                                                     aGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACgg                                               <  2:4134951/150‑1 (MQ=255)
                                                                                                      gCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACCGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGc                                              <  2:1658836/150‑1 (MQ=255)
                                                                                                        ggCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGa                                            >  2:4067258/1‑150 (MQ=255)
                                                                                                        ggCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGa                                            <  1:248922/150‑1 (MQ=255)
                                                                                                        ggCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGa                                            >  2:2049628/1‑150 (MQ=255)
                                                                                                          ccAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACg                                          <  1:1915644/150‑1 (MQ=255)
                                                                                                            aGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGc                                        >  1:910276/1‑150 (MQ=255)
                                                                                                             gTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCt                                       >  1:2530226/1‑150 (MQ=255)
                                                                                                             gTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCt                                       <  2:3552398/150‑1 (MQ=255)
                                                                                                                gACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGc                                    <  2:4485891/150‑1 (MQ=255)
                                                                                                                gACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGc                                    >  2:2715533/1‑150 (MQ=255)
                                                                                                                gACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGc                                    <  2:3685493/150‑1 (MQ=255)
                                                                                                                   aTAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCaaa                                 >  2:3899989/1‑150 (MQ=255)
                                                                                                                      gCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAat                              >  2:851713/1‑150 (MQ=255)
                                                                                                                       ctctGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatg                             <  2:2958733/150‑1 (MQ=255)
                                                                                                                       ctctGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatg                             <  2:3089590/150‑1 (MQ=255)
                                                                                                                        tctGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatga                            >  2:2645054/1‑150 (MQ=255)
                                                                                                                         ctGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatgat                           >  2:2819545/1‑150 (MQ=255)
                                                                                                                         ctGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatgat                           <  2:3148706/150‑1 (MQ=255)
                                                                                                                          tGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatgatg                          <  2:676574/150‑1 (MQ=255)
                                                                                                                           gCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATgg                         <  1:4535802/150‑1 (MQ=255)
                                                                                                                           gCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATgg                         <  1:179878/150‑1 (MQ=255)
                                                                                                                           gCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATgg                         >  1:2292159/1‑150 (MQ=255)
                                                                                                                            ccACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGc                        <  1:1381404/150‑1 (MQ=255)
                                                                                                                              aCTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCa                      <  1:1362576/150‑1 (MQ=255)
                                                                                                                                  gCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCg                  >  1:4398390/1‑150 (MQ=255)
                                                                                                                                   cTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgt                 >  1:2125107/1‑150 (MQ=255)
                                                                                                                                      tCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtagt              >  1:1428507/1‑150 (MQ=255)
                                                                                                                                         tGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTaaa           >  2:4137655/1‑150 (MQ=255)
                                                                                                                                          gAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAGGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAAc          >  1:1046394/1‑150 (MQ=255)
                                                                                                                                          gAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAAc          <  1:4183418/150‑1 (MQ=255)
                                                                                                                                          gAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAAc          <  1:2074625/150‑1 (MQ=255)
                                                                                                                                           aGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACg         >  2:3413964/1‑150 (MQ=255)
                                                                                                                                            gCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACTGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGTTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACgg        <  2:4320075/150‑1 (MQ=255)
                                                                                                                                            gCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACgg        <  2:154689/150‑1 (MQ=255)
                                                                                                                                            gCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACgg        <  2:636281/150‑1 (MQ=255)
                                                                                                                                                taatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCaga    >  2:912264/1‑150 (MQ=255)
                                                                                                                                                taatATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAACTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCaga    >  2:2075608/1‑150 (MQ=255)
                                                                                                                                                 aatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCagag   <  1:3197468/150‑1 (MQ=255)
                                                                                                                                                 aatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCagag   >  2:3293918/1‑150 (MQ=255)
                                                                                                                                                  atCTCTTTACCCAGCGTATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTAATGGTTAGCATGTAATGTCATCTCTGGCACGAGATCGTCACGGCTACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAgt  <  1:3514679/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGT  >  CP009273/3956017‑3956312

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: