Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,809,621 C→A 23.8% 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,6210CA23.8% 527.2 / 149.2 223D60Y (GAT→TAT) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (82/87);  new base A (38/15);  total (120/103)
Fisher's exact test for biased strand distribution p-value = 4.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 233 total aligned reads displayed.

TGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGG  >  CP009273/2809472‑2809770
                                                                                                                                                     |                                                                                                                                                     
tGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATc                                                                                                                                                       >  2:941310/1‑150 (MQ=255)
tGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATc                                                                                                                                                       <  1:1142324/150‑1 (MQ=255)
   gATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATAGTGCGCGCCGTATAACGCTGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGgt                                                                                                                                                    <  2:1507377/150‑1 (MQ=255)
    aTGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGgta                                                                                                                                                   <  2:1742978/150‑1 (MQ=255)
             ggCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCAt                                                                                                                                          >  2:2362197/1‑150 (MQ=255)
                aTACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTtg                                                                                                                                       >  1:4317369/1‑150 (MQ=255)
                    ttAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCg                                                                                                                                   <  1:408371/150‑1 (MQ=255)
                      aaCGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGTTGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTc                                                                                                                                 <  1:3960833/150‑1 (MQ=255)
                      aaCGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTc                                                                                                                                 >  2:2558201/1‑150 (MQ=255)
                         ggCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGt                                                                                                                              >  2:994738/1‑150 (MQ=255)
                         ggCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGt                                                                                                                              >  1:4487880/1‑150 (MQ=255)
                         ggCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGt                                                                                                                              >  2:1368873/1‑150 (MQ=255)
                             acaTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGCGCgg                                                                                                                          >  2:3831407/1‑150 (MQ=255)
                              caTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCCCCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGa                                                                                                                         <  2:5186396/150‑1 (MQ=255)
                              caTCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGa                                                                                                                         <  2:2002497/150‑1 (MQ=255)
                                tCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAAc                                                                                                                       >  1:2683520/1‑150 (MQ=255)
                                 ccGCTCATCGCCCATAGTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAAcc                                                                                                                      <  1:4023205/150‑1 (MQ=255)
                                   gCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTa                                                                                                                    <  2:4580422/150‑1 (MQ=255)
                                     tCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAAt                                                                                                                  <  1:1925339/150‑1 (MQ=255)
                                     tCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAAt                                                                                                                  <  1:2706498/150‑1 (MQ=255)
                                        tCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCATTGCGGAACCTAATGCt                                                                                                               <  2:2282572/150‑1 (MQ=255)
                                         cGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCtt                                                                                                              >  1:958347/1‑150 (MQ=255)
                                          gCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTc                                                                                                             <  2:3498583/150‑1 (MQ=255)
                                           cccATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCa                                                                                                            >  2:443176/1‑150 (MQ=255)
                                            ccATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAg                                                                                                           >  1:752606/1‑150 (MQ=255)
                                            ccATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAg                                                                                                           <  2:525006/150‑1 (MQ=255)
                                             cATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAgg                                                                                                          <  1:634947/150‑1 (MQ=255)
                                             cATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAgg                                                                                                          >  2:3903803/1‑150 (MQ=255)
                                              atatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGa                                                                                                         >  2:1045351/1‑150 (MQ=255)
                                               tatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGAt                                                                                                        <  1:2229140/150‑1 (MQ=255)
                                               tatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGAt                                                                                                        <  2:3997774/150‑1 (MQ=255)
                                               tatTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGAt                                                                                                        >  2:2763043/1‑150 (MQ=255)
                                                  tGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGAc                                                                                                     <  2:3068652/150‑1 (MQ=255)
                                                   gcgcgcCGTATAACGATGCAGATCTCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGAcc                                                                                                    <  2:3754765/150‑1 (MQ=255)
                                                    cgcgcCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCt                                                                                                   <  1:1802837/150‑1 (MQ=255)
                                                       gcCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCtgtt                                                                                                <  2:2811506/150‑1 (MQ=255)
                                                         cGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCtgttgt                                                                                              >  1:4480222/1‑150 (MQ=255)
                                                          gTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTgttgtt                                                                                             <  1:1598467/150‑1 (MQ=255)
                                                          gTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTgttgtt                                                                                             <  2:2050286/150‑1 (MQ=255)
                                                            ataACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGc                                                                                           >  1:3419849/1‑150 (MQ=255)
                                                              aaCGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCa                                                                                         <  1:3342172/150‑1 (MQ=255)
                                                                       gATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGAGTGCCAt                                                                                >  2:3919508/1‑150 (MQ=255)
                                                                       gATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCAt                                                                                >  2:5058241/1‑150 (MQ=255)
                                                                       gATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGGAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTGGCCAGTGTGCCAt                                                                                >  2:3810617/1‑150 (MQ=255)
                                                                        aTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCACTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATc                                                                               <  2:2953494/150‑1 (MQ=255)
                                                                        aTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATc                                                                               <  1:1414894/150‑1 (MQ=255)
                                                                        aTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATc                                                                               >  2:2109138/1‑150 (MQ=255)
                                                                        aTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATc                                                                               >  2:2686787/1‑150 (MQ=255)
                                                                         tCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCa                                                                              <  1:2816299/150‑1 (MQ=255)
                                                                          cgcgCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTCTTGTTGCCAGTGTGCCATCAg                                                                             >  2:3608954/1‑150 (MQ=255)
                                                                          cgcgCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAg                                                                             >  2:3162472/1‑150 (MQ=255)
                                                                            cgcATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCa                                                                           >  1:824888/1‑150 (MQ=255)
                                                                              cATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCAATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCAtt                                                                         >  2:2870374/1‑150 (MQ=255)
                                                                                 aaaGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAac                                                                      <  1:5048643/150‑1 (MQ=255)
                                                                                 aaaGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAac                                                                      >  2:1589221/1‑150 (MQ=255)
                                                                                  aaGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAaca                                                                     <  2:825848/150‑1 (MQ=255)
                                                                                  aaGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAaca                                                                     >  2:1573533/1‑150 (MQ=255)
                                                                                     gTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGc                                                                  <  1:3234915/150‑1 (MQ=255)
                                                                                            tatGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATAAACACGCAAAGttt                                                           >  2:4775422/1‑150 (MQ=255)
                                                                                             atGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTc                                                          <  1:4436071/150‑1 (MQ=255)
                                                                                              tGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTcg                                                         <  2:623246/150‑1 (MQ=255)
                                                                                               gTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTcgc                                                        >  1:3840723/1‑150 (MQ=255)
                                                                                                  cAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCt                                                     >  2:1904390/1‑150 (MQ=255)
                                                                                                    atatCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTcc                                                   >  2:2661244/1‑150 (MQ=255)
                                                                                                      atCACCCTCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAg                                                 <  2:1124144/150‑1 (MQ=255)
                                                                                                      atCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAg                                                 >  2:1963690/1‑150 (MQ=255)
                                                                                                            aTCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACg                                           >  1:2572952/1‑150 (MQ=255)
                                                                                                             tCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGc                                          >  1:1166224/1‑150 (MQ=255)
                                                                                                               cACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTg                                        <  2:3225729/150‑1 (MQ=255)
                                                                                                                aCTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGt                                       <  1:1544139/150‑1 (MQ=255)
                                                                                                                     tgtAATGAATTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATccc                                  <  1:3903803/150‑1 (MQ=255)
                                                                                                                       tAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCt                                <  2:4026711/150‑1 (MQ=255)
                                                                                                                        aaTGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCtt                               <  2:3338475/150‑1 (MQ=255)
                                                                                                                         aTGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCttt                              <  1:4516911/150‑1 (MQ=255)
                                                                                                                            aaTTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAAc                           >  1:534827/1‑150 (MQ=255)
                                                                                                                            aaTTCCAGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAAc                           >  1:534899/1‑150 (MQ=255)
                                                                                                                                ccAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCt                       <  1:1494059/150‑1 (MQ=255)
                                                                                                                                  aGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGa                     <  1:605740/150‑1 (MQ=255)
                                                                                                                                  aGCAATGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGa                     >  1:51681/1‑150 (MQ=255)
                                                                                                                                    taaTGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAgg                   >  2:4555786/2‑150 (MQ=255)
                                                                                                                                    taaTGCTTCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAgg                   >  2:4555836/2‑150 (MQ=255)
                                                                                                                                     aaTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGa                  >  1:1304549/1‑150 (MQ=255)
                                                                                                                                     aaTGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGa                  >  2:5168058/1‑150 (MQ=255)
                                                                                                                                         cTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGtt              >  1:4956842/1‑150 (MQ=255)
                                                                                                                                         cTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGtt              >  1:4963787/1‑150 (MQ=255)
                                                                                                                                          ttCCGCAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGttt             >  2:32786/1‑150 (MQ=255)
                                                                                                                                           tCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGtttt            >  1:3301663/1‑150 (MQ=255)
                                                                                                                                           tCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGtttt            >  2:513936/1‑150 (MQ=255)
                                                                                                                                             cGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTc          >  1:1611301/1‑150 (MQ=255)
                                                                                                                                              gCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcc         <  1:3383068/150‑1 (MQ=255)
                                                                                                                                              gCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcc         <  1:3641817/150‑1 (MQ=255)
                                                                                                                                               cAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcca        <  1:2362197/150‑1 (MQ=255)
                                                                                                                                               cAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcca        >  2:4660241/1‑150 (MQ=255)
                                                                                                                                               cAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcca        >  2:3505208/1‑150 (MQ=255)
                                                                                                                                               cAAAATAGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTcca        <  1:631413/150‑1 (MQ=255)
                                                                                                                                                aaaaTCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccag       >  2:3882541/1‑150 (MQ=255)
                                                                                                                                                 aaaTCGGTAGAAATCCATTTGTGCGTCAGTGCGGAACCTAATGCATCAGGATGACCTGTTGTTGCCAGTGAGCCATCAGCATTAACACGCAAACTTTCGCGCTCCAGCCCACGCTGTAACCCCTTTAACGCCTGAGGATGTTTTTccagc      >  2:4308649/1‑150 (MQ=255)
                                                                                                                                                    tCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTccagccag   >  1:200393/1‑150 (MQ=255)
                                                                                                                                                     cGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCACAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAgccagg  >  1:4220877/1‑150 (MQ=255)
                                                                                                                                                     cGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAgccagg  <  2:173071/150‑1 (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: