Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,144,163 C→A 11.9% A205S (GCT→TCT)  argE ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,144,1630CA11.9% 1153.6 / 84.2 332A205S (GCT→TCT) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (151/141);  new base A (19/21);  total (170/162)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

CAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTT  >  CP009273/4144015‑4144309
                                                                                                                                                    |                                                                                                                                                  
cAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCg                                                                                                                                                   >  2:4363932/1‑150 (MQ=255)
 aaaTACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGt                                                                                                                                                  >  2:2148898/1‑150 (MQ=255)
  aaTACGGTTAGAATCGTCGCCACAGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCAATGTAAACGCTTCGTAGTGATAACGTTCTTTAAGGTTATCGCGCAATTGCAAAATATGCCCGATTGCGTCGTGCATTAGTTCGATAGAGtt                                                                                                                                                 <  2:152929/150‑1 (MQ=255)
  aaTACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGtt                                                                                                                                                 <  1:158337/150‑1 (MQ=255)
    tACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTaa                                                                                                                                               <  1:1737060/150‑1 (MQ=255)
       ggTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTc                                                                                                                                            >  2:3588393/1‑150 (MQ=255)
        gTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTcc                                                                                                                                           <  2:1816879/150‑1 (MQ=255)
           aGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCcgcg                                                                                                                                        >  1:277092/1‑150 (MQ=255)
                cGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTg                                                                                                                                   >  1:1202991/1‑150 (MQ=255)
                 gTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTgg                                                                                                                                  <  2:2623179/150‑1 (MQ=255)
                  tCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGa                                                                                                                                 <  1:1492330/150‑1 (MQ=255)
                   cGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGAt                                                                                                                                >  1:3099870/1‑150 (MQ=255)
                    gCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATc                                                                                                                               >  2:1301112/1‑150 (MQ=255)
                    gCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATc                                                                                                                               >  1:2257550/1‑150 (MQ=255)
                      cACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGc                                                                                                                             >  1:2659302/1‑150 (MQ=255)
                       aCCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCt                                                                                                                            <  1:600902/150‑1 (MQ=255)
                            gAATATGCCCGAGGGTGAGCGTAGGGGATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGt                                                                                                                       <  2:2517692/150‑1 (MQ=255)
                              atatGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGc                                                                                                                     <  2:3728998/150‑1 (MQ=255)
                               tatGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGcc                                                                                                                    <  1:4316220/150‑1 (MQ=255)
                                atGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGccc                                                                                                                   <  2:2098556/150‑1 (MQ=255)
                                 tGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGcccc                                                                                                                  <  2:688587/150‑1 (MQ=255)
                                 tGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGcccc                                                                                                                  <  2:1651703/150‑1 (MQ=255)
                                    ccGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGAc                                                                                                               >  2:1219219/1‑150 (MQ=255)
                                      gAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTg                                                                                                             >  1:3362566/1‑150 (MQ=255)
                                         gTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGcc                                                                                                          >  1:4029740/1‑150 (MQ=255)
                                         gTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGcc                                                                                                          <  2:2261028/150‑1 (MQ=255)
                                          ttGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGccc                                                                                                         <  1:625340/150‑1 (MQ=255)
                                          ttGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGccc                                                                                                         >  2:881708/1‑150 (MQ=255)
                                           tGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCt                                                                                                        <  1:1761819/150‑1 (MQ=255)
                                            gAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTg                                                                                                       <  2:2002258/150‑1 (MQ=255)
                                             aGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGa                                                                                                      <  1:4632950/150‑1 (MQ=255)
                                                gTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATa                                                                                                   <  2:2563952/150‑1 (MQ=255)
                                                 tAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATAc                                                                                                  <  2:3808191/150‑1 (MQ=255)
                                                 tAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATAc                                                                                                  >  1:107924/1‑150 (MQ=255)
                                                 tAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATAc                                                                                                  >  2:2330140/1‑150 (MQ=255)
                                                      tATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATg                                                                                             <  2:3262682/150‑1 (MQ=255)
                                                            aCGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTa                                                                                       <  1:3591212/150‑1 (MQ=255)
                                                            aCGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTa                                                                                       >  2:93835/1‑150 (MQ=255)
                                                            aCGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTa                                                                                       >  2:3012024/1‑150 (MQ=255)
                                                               gTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTagag                                                                                    <  2:1791496/150‑1 (MQ=255)
                                                                 aaaCGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGat                                                                                  <  2:73500/150‑1 (MQ=255)
                                                                   aCGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGatat                                                                                <  2:3973674/150‑1 (MQ=255)
                                                                    cGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATg                                                                               <  1:384755/150‑1 (MQ=255)
                                                                    cGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATg                                                                               <  1:2846260/150‑1 (MQ=255)
                                                                     gCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGa                                                                              >  2:544458/1‑150 (MQ=255)
                                                                     gCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGa                                                                              >  2:1726091/1‑150 (MQ=255)
                                                                         cGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCtt                                                                          <  2:3853404/150‑1 (MQ=255)
                                                                           tAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTa                                                                        <  1:3284523/150‑1 (MQ=255)
                                                                            aGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTAt                                                                       >  1:4401693/1‑150 (MQ=255)
                                                                            aGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTAt                                                                       >  1:190276/1‑150 (MQ=255)
                                                                            aGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCGTTAt                                                                       >  1:4084902/1‑150 (MQ=255)
                                                                                 tAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTgcg                                                                  <  2:3544821/150‑1 (MQ=255)
                                                                                     gTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTa                                                              >  2:2125804/1‑150 (MQ=255)
                                                                                      ttCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTAc                                                             <  1:994639/150‑1 (MQ=255)
                                                                                       tCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTAcc                                                            <  2:495482/150‑1 (MQ=255)
                                                                                         tttCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCgg                                                          >  1:3456228/1‑150 (MQ=255)
                                                                                            cAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTg                                                       <  2:3143412/150‑1 (MQ=255)
                                                                                             aGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGt                                                      >  2:2439082/1‑150 (MQ=255)
                                                                                              ggTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTa                                                     <  1:421436/150‑1 (MQ=255)
                                                                                                  aTCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGa                                                 >  1:3689176/1‑150 (MQ=255)
                                                                                                  aTCGCGCAATTGCAAAATAGGCCCGATGGCGTCGTGCATTAGTGCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGACCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGa                                                 <  2:2802351/150‑1 (MQ=255)
                                                                                                      cgcAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTc                                             <  1:2976025/150‑1 (MQ=255)
                                                                                                      cgcAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTc                                             >  1:3199213/1‑150 (MQ=255)
                                                                                                        cAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCgg                                           >  1:1970467/1‑150 (MQ=255)
                                                                                                         aaTTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGt                                          >  2:4158253/1‑150 (MQ=255)
                                                                                                          atggCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGtt                                         <  1:764885/147‑1 (MQ=255)
                                                                                                           ttGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTc                                        <  2:2919946/150‑1 (MQ=255)
                                                                                                           ttGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTc                                        <  2:1124317/150‑1 (MQ=255)
                                                                                                            tGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCg                                       <  2:2953015/150‑1 (MQ=255)
                                                                                                               aaaaTATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCa                                    <  1:134055/150‑1 (MQ=255)
                                                                                                                aaaTATGCCCGCTGGCGTCGTGCATTAGTTCGATAGCGTTATCTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCaa                                   <  2:4559115/150‑1 (MQ=255)
                                                                                                                aaaTATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCaa                                   >  2:149557/1‑150 (MQ=255)
                                                                                                                  ctatGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAatg                                 <  1:4468312/149‑1 (MQ=255)
                                                                                                                     tGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAatgatg                              >  1:2209731/1‑150 (MQ=255)
                                                                                                                       cccGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGgc                            >  1:1007015/1‑150 (MQ=255)
                                                                                                                          gATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCa                         >  1:2896065/1‑150 (MQ=255)
                                                                                                                             ggCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATc                      >  2:3906103/1‑150 (MQ=255)
                                                                                                                             ggCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATc                      >  1:2518448/1‑150 (MQ=255)
                                                                                                                             ggCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATc                      >  2:3444208/1‑150 (MQ=255)
                                                                                                                              gCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATcc                     <  2:54984/150‑1 (MQ=255)
                                                                                                                              gCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATcc                     <  2:759150/150‑1 (MQ=255)
                                                                                                                                gtcgtGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCgg                   >  2:1351219/1‑150 (MQ=255)
                                                                                                                                  cgtGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGgc                 >  2:4195512/1‑150 (MQ=255)
                                                                                                                                      cATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAg             <  2:1409242/150‑1 (MQ=255)
                                                                                                                                       aTTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAgg            >  2:1285877/1‑150 (MQ=255)
                                                                                                                                         tAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCACTCCGGGCGCAGGGc          >  2:4445484/1‑150 (MQ=255)
                                                                                                                                         tAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGc          <  1:3548117/150‑1 (MQ=255)
                                                                                                                                          aGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCg         >  1:2350587/1‑150 (MQ=255)
                                                                                                                                          aGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCg         >  1:3800632/1‑150 (MQ=255)
                                                                                                                                          aGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCg         <  1:2493052/150‑1 (MQ=255)
                                                                                                                                          aGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCg         <  2:2673824/150‑1 (MQ=255)
                                                                                                                                          aGTTCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCg         <  1:3426229/150‑1 (MQ=255)
                                                                                                                                            ttCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGt       >  2:1876442/1‑150 (MQ=255)
                                                                                                                                            ttCGATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGt       <  2:277092/150‑1 (MQ=255)
                                                                                                                                             tCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGGAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTa      <  1:617547/150‑1 (MQ=255)
                                                                                                                                               gATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt    >  1:1667049/1‑150 (MQ=255)
                                                                                                                                               gATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt    >  2:866212/1‑150 (MQ=255)
                                                                                                                                               gATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt    >  1:3965576/1‑150 (MQ=255)
                                                                                                                                               gATAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTAGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt    >  1:2077623/1‑150 (MQ=255)
                                                                                                                                                 tAGAGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGttt  <  1:1650205/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAATACGGTTAGAAGCGTCGCCACCGTGAATATGCCCGAGGTTGAGCGTAGGGTATGGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTT  >  CP009273/4144015‑4144309

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: