Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,809,621 C→A 15.6% D60Y (GAT→TAT)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,809,6210CA15.6% 802.4 / 102.1 258D60Y (GAT→TAT) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (123/94);  new base A (19/22);  total (142/116)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGG  >  CP009273/2809472‑2809770
                                                                                                                                                     |                                                                                                                                                     
tGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATc                                                                                                                                                       <  1:1133496/150‑1 (MQ=255)
tGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATc                                                                                                                                                       <  2:1573952/150‑1 (MQ=255)
      gTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGgtagt                                                                                                                                                 <  1:1192430/150‑1 (MQ=255)
        aGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTaa                                                                                                                                               <  1:1538312/150‑1 (MQ=255)
         gCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAAt                                                                                                                                              >  1:233625/1‑150 (MQ=255)
          cATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATc                                                                                                                                             <  2:1761708/150‑1 (MQ=255)
          cATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATc                                                                                                                                             <  2:3464965/150‑1 (MQ=255)
           aTGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATcc                                                                                                                                            >  1:4647579/1‑150 (MQ=255)
            tGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCa                                                                                                                                           >  2:1337803/1‑150 (MQ=255)
            tGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCa                                                                                                                                           >  2:3731184/1‑150 (MQ=255)
             ggCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCAt                                                                                                                                          >  2:1779327/1‑150 (MQ=255)
               cATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCAttt                                                                                                                                        >  2:3168773/1‑150 (MQ=255)
               cATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCAttt                                                                                                                                        >  2:3326838/1‑150 (MQ=255)
                  aCTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTtgtg                                                                                                                                     <  1:394232/150‑1 (MQ=255)
                   cTTCACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGc                                                                                                                                    <  2:2878880/150‑1 (MQ=255)
                   cTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGc                                                                                                                                    <  1:2568558/150‑1 (MQ=255)
                   cTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGc                                                                                                                                    <  1:343755/150‑1 (MQ=255)
                    ttAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCg                                                                                                                                   >  2:3671657/1‑150 (MQ=255)
                    ttAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCg                                                                                                                                   <  1:1057768/150‑1 (MQ=255)
                      aaCGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTc                                                                                                                                 >  1:4670548/1‑150 (MQ=255)
                      aaCGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTc                                                                                                                                 >  2:742877/1‑150 (MQ=255)
                         ggCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGt                                                                                                                              >  1:2459922/1‑150 (MQ=255)
                          gCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTg                                                                                                                             <  2:1494347/150‑1 (MQ=255)
                          gCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTg                                                                                                                             <  2:4500345/150‑1 (MQ=255)
                            cacaTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCg                                                                                                                           >  2:3853693/1‑150 (MQ=255)
                            cacaTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCg                                                                                                                           >  2:2386296/1‑150 (MQ=255)
                            cacaTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCg                                                                                                                           <  1:1894033/150‑1 (MQ=255)
                              caTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGa                                                                                                                         >  1:250605/1‑150 (MQ=255)
                               aTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGaa                                                                                                                        >  2:1493960/1‑150 (MQ=255)
                               aTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGaa                                                                                                                        >  2:4716421/1‑150 (MQ=255)
                               aTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGaa                                                                                                                        <  2:4280546/150‑1 (MQ=255)
                               aTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGaa                                                                                                                        <  2:685826/150‑1 (MQ=255)
                                tCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAAc                                                                                                                       >  2:3882416/1‑150 (MQ=255)
                                 ccGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAAcc                                                                                                                      <  1:1376682/150‑1 (MQ=255)
                                  cGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCt                                                                                                                     <  2:3336564/150‑1 (MQ=255)
                                    cTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTaa                                                                                                                   >  2:328323/1‑150 (MQ=255)
                                     taaTCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGATTCAGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTTCGGAACCTAAt                                                                                                                  >  1:2627338/3‑150 (MQ=255)
                                     tCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAAt                                                                                                                  >  1:348105/1‑150 (MQ=255)
                                     tCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAAt                                                                                                                  >  1:1730569/1‑150 (MQ=255)
                                       aTCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGc                                                                                                                >  2:2288044/1‑150 (MQ=255)
                                       aTCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATACCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGc                                                                                                                <  2:2384820/150‑1 (MQ=255)
                                        tCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCt                                                                                                               >  1:3308754/1‑150 (MQ=255)
                                          gCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTc                                                                                                             >  1:4475058/1‑150 (MQ=255)
                                          gCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTc                                                                                                             >  1:4473552/1‑150 (MQ=255)
                                            ccATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAg                                                                                                           >  2:558271/1‑150 (MQ=255)
                                              atatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGa                                                                                                         >  2:1979016/1‑150 (MQ=255)
                                               tatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAAGATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGAt                                                                                                        <  2:4173473/150‑1 (MQ=255)
                                                  tGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGAc                                                                                                     >  2:271518/1‑150 (MQ=255)
                                                    cgcgcCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCt                                                                                                   <  2:161959/150‑1 (MQ=255)
                                                    cgcgcCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCt                                                                                                   <  1:1077504/150‑1 (MQ=255)
                                                       gcCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCtgtt                                                                                                <  1:4587853/150‑1 (MQ=255)
                                                         cGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCtgttgt                                                                                              >  1:2585873/1‑150 (MQ=255)
                                                           tataACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGttgttg                                                                                            <  1:1685250/150‑1 (MQ=255)
                                                           tataACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGttgttg                                                                                            >  1:4373726/1‑150 (MQ=255)
                                                             taACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGcc                                                                                          >  1:604726/1‑150 (MQ=255)
                                                              aaCGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCa                                                                                         <  2:4543760/150‑1 (MQ=255)
                                                                cGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGGAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTCCCAgt                                                                                       >  2:2944982/1‑150 (MQ=255)
                                                                 gATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAgtg                                                                                      >  2:700364/1‑150 (MQ=255)
                                                                  aTGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAgtgt                                                                                     <  2:3542090/150‑1 (MQ=255)
                                                                     cAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGcc                                                                                  >  2:2576739/1‑150 (MQ=255)
                                                                       gATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCAt                                                                                >  1:813174/1‑150 (MQ=255)
                                                                       gATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCAt                                                                                >  1:1259018/1‑150 (MQ=255)
                                                                        aTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATc                                                                               <  2:1119467/150‑1 (MQ=255)
                                                                         tCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTACTGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCa                                                                              <  2:821892/150‑1 (MQ=255)
                                                                         tCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCa                                                                              <  2:1310544/150‑1 (MQ=255)
                                                                             gcATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCAt                                                                          >  1:172410/1‑150 (MQ=255)
                                                                              cATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCAtt                                                                         >  2:3471511/1‑150 (MQ=255)
                                                                               aTAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTa                                                                        >  2:2348363/1‑150 (MQ=255)
                                                                               aTAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTa                                                                        >  1:3708410/1‑150 (MQ=255)
                                                                                  aaGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAaca                                                                     >  1:2080153/1‑150 (MQ=255)
                                                                                     gTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGc                                                                  <  1:761445/150‑1 (MQ=255)
                                                                                        aGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCaaa                                                               >  2:3568255/1‑150 (MQ=255)
                                                                                         gCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAg                                                              <  1:216520/150‑1 (MQ=255)
                                                                                          cATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGt                                                             <  1:3671657/150‑1 (MQ=255)
                                                                                           atatGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGtt                                                            <  1:3381276/150‑1 (MQ=255)
                                                                                            tatGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGttt                                                           <  1:4096330/150‑1 (MQ=255)
                                                                                                   aaTATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTc                                                    >  2:1668026/1‑150 (MQ=255)
                                                                                                        cACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGcc                                               >  1:3861833/1‑150 (MQ=255)
                                                                                                        cACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGcc                                               >  2:559902/1‑150 (MQ=255)
                                                                                                            aTCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACg                                           <  2:3420530/150‑1 (MQ=255)
                                                                                                            aTCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACg                                           <  2:4171113/150‑1 (MQ=255)
                                                                                                                        aaTGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCtt                               <  1:859453/150‑1 (MQ=255)
                                                                                                                         aTGAATTCCAGCAATGCTTCCGCAAAAACGGTAGTAATCCATTTGTGCGTCAGTGCTGAACCAAAAGCTTCAGGATGACCTGTTGTTGCCAGTGTTCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCttt                              >  2:463936/1‑150 (MQ=255)
                                                                                                                          tGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTa                             <  1:503378/150‑1 (MQ=255)
                                                                                                                             aTTCCAGCAATGCTTCCGCAACATCGGTAGTAATCCTTTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAAct                          <  2:2627338/150‑2 (MQ=255)
                                                                                                                             aTTCCAGCAAGGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCGTTAACg                          >  2:437704/1‑150 (MQ=255)
                                                                                                                              ttCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGc                         >  1:1962991/1‑150 (MQ=255)
                                                                                                                                   gCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAg                    <  1:1117737/150‑1 (MQ=255)
                                                                                                                                    cAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAgg                   >  2:3998469/1‑150 (MQ=255)
                                                                                                                                    cAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAgg                   >  2:2269479/1‑150 (MQ=255)
                                                                                                                                     aaTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGa                  >  1:3138557/1‑150 (MQ=255)
                                                                                                                                      aTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGAt                 <  2:4151971/150‑1 (MQ=255)
                                                                                                                                      aTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGAt                 >  2:3021606/1‑150 (MQ=255)
                                                                                                                                      aTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGAt                 <  1:1493960/150‑1 (MQ=255)
                                                                                                                                       ggCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATg                >  1:4016502/2‑150 (MQ=255)
                                                                                                                                           tCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGtttt            >  1:2232192/1‑150 (MQ=255)
                                                                                                                                                aaaaTCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccag       <  1:2851863/150‑1 (MQ=255)
                                                                                                                                                aaaaTCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccag       <  1:3326838/150‑1 (MQ=255)
                                                                                                                                                aaaaTCGGTAGTAATCCATTTGTGCCTCAGTGCGGAACCTAATGCTTCAGTATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCACCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccag       >  1:3602064/1‑150 (MQ=255)
                                                                                                                                                  aaTCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccagcc     >  1:1795225/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGG  >  CP009273/2809472‑2809770

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: