Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,646 G→T 9.8% R385S (CGC→AGC)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,6460GT 9.8% 1093.5 / 57.3 303R385S (CGC→AGC) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (124/148);  new base T (12/19);  total (136/167)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

ATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808497‑2808795
                                                                                                                                                     |                                                                                                                                                     
atcaCACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCg                                                                                                                                                       >  1:2450125/1‑150 (MQ=255)
  cacacaCTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTAc                                                                                                                                                     >  1:757454/1‑150 (MQ=255)
   acacacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGAcc                                                                                                                                                    <  2:1809469/150‑1 (MQ=255)
    cacacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCt                                                                                                                                                   <  1:942470/150‑1 (MQ=255)
     acacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCtt                                                                                                                                                  >  1:2239124/1‑150 (MQ=255)
      cacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTc                                                                                                                                                 <  2:2417881/150‑1 (MQ=255)
      cacTTTCTGATACGCTTCGCCGCCGTGAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTc                                                                                                                                                 >  1:1596419/1‑150 (MQ=255)
       acTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCg                                                                                                                                                >  1:4452536/1‑150 (MQ=255)
        cgttCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGa                                                                                                                                               <  2:2316616/148‑1 (MQ=255)
        cTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGa                                                                                                                                               >  2:321504/1‑150 (MQ=255)
           tCTGAGACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGa                                                                                                                                            <  2:3679817/150‑1 (MQ=255)
            cTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGAt                                                                                                                                           >  1:3162400/1‑150 (MQ=255)
             tGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATc                                                                                                                                          >  2:3849122/1‑150 (MQ=255)
               aTACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCAc                                                                                                                                        <  2:219396/150‑1 (MQ=255)
                    cTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGt                                                                                                                                   <  1:4629672/150‑1 (MQ=255)
                    cTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGt                                                                                                                                   >  2:3900029/1‑150 (MQ=255)
                      tCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTc                                                                                                                                 <  2:577781/150‑1 (MQ=255)
                       cgccgcCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTcc                                                                                                                                <  2:2534959/150‑1 (MQ=255)
                       cgccgcCGTTAATACAATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTcc                                                                                                                                >  1:3467409/1‑150 (MQ=255)
                        gccgccGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCa                                                                                                                               <  2:4216273/150‑1 (MQ=255)
                        gccgccGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCa                                                                                                                               <  1:1822998/150‑1 (MQ=255)
                        gccgccGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCa                                                                                                                               >  1:1365201/1‑150 (MQ=255)
                         ccgccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAg                                                                                                                              <  2:1669111/150‑1 (MQ=255)
                           gccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGtt                                                                                                                            <  1:3003702/150‑1 (MQ=255)
                            ccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTa                                                                                                                           <  2:2374114/150‑1 (MQ=255)
                                  aTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCgt                                                                                                                     <  2:2429725/150‑1 (MQ=255)
                                      tATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTAc                                                                                                                 >  1:595054/1‑150 (MQ=255)
                                        tCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAg                                                                                                               <  2:3294350/150‑1 (MQ=255)
                                           aGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCa                                                                                                            >  2:491756/1‑150 (MQ=255)
                                               tttGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGtt                                                                                                        <  1:667906/150‑1 (MQ=255)
                                                   cgcgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCt                                                                                                    <  1:3220866/150‑1 (MQ=255)
                                                     cgACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTAc                                                                                                  >  1:1038819/1‑150 (MQ=255)
                                                      aaCGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACt                                                                                                 <  1:4289525/149‑1 (MQ=255)
                                                          cgTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACTGGAACTGTGCGGTTTCGCAGCCGATACACATCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCAGGTTCCAGTTAACGCGTGTAAAGGCAAGTTCGCTACTGCTc                                                                                             <  1:2945740/150‑1 (MQ=255)
                                                          cgTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTc                                                                                             <  1:190197/150‑1 (MQ=255)
                                                           gTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCa                                                                                            <  2:2607265/150‑1 (MQ=255)
                                                             ttCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAtt                                                                                          <  1:3667318/150‑1 (MQ=255)
                                                             ttCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAtt                                                                                          >  2:3134789/1‑150 (MQ=255)
                                                             ttCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAtt                                                                                          <  1:3920690/150‑1 (MQ=255)
                                                              tCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAttt                                                                                         >  2:1113902/1‑150 (MQ=255)
                                                                aGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTcc                                                                                       >  1:3090247/1‑150 (MQ=255)
                                                                aGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTcc                                                                                       <  2:1738646/150‑1 (MQ=255)
                                                                 gATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGCTCGCTACTGCTCATTTCCg                                                                                      <  1:2677372/150‑1 (MQ=255)
                                                                         aaCAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATca                                                                              <  2:2670243/150‑1 (MQ=255)
                                                                           cAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagc                                                                            <  2:2391664/150‑1 (MQ=255)
                                                                           cAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagc                                                                            >  1:277541/1‑150 (MQ=255)
                                                                            aGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagcc                                                                           >  2:3082275/1‑150 (MQ=255)
                                                                              aTCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagccag                                                                         >  1:1494828/1‑150 (MQ=255)
                                                                              aTCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagccag                                                                         <  1:796851/150‑1 (MQ=255)
                                                                                cTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAgcg                                                                       >  1:3713224/1‑150 (MQ=255)
                                                                                     cccACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCacac                                                                  <  1:3846359/150‑1 (MQ=255)
                                                                                     cccACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCacac                                                                  >  2:3618183/1‑150 (MQ=255)
                                                                                        aCCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAg                                                               <  2:2922865/150‑1 (MQ=255)
                                                                                         ccTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGa                                                              >  1:177635/1‑150 (MQ=255)
                                                                                           tGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGAcc                                                            <  1:877677/150‑1 (MQ=255)
                                                                                           tGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGAcc                                                            <  2:3780685/150‑1 (MQ=255)
                                                                                            gCGGTAACGGGAACTGTGCGGTTTCGCAGCCGAGACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCa                                                           <  2:154634/150‑1 (MQ=255)
                                                                                             cGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCAt                                                          >  2:1404438/1‑150 (MQ=255)
                                                                                             cGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCAt                                                          <  2:2768118/150‑1 (MQ=255)
                                                                                             cGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCAt                                                          <  1:2102735/150‑1 (MQ=255)
                                                                                              ggTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATa                                                         >  1:2488477/1‑150 (MQ=255)
                                                                                               gTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaa                                                        <  2:2977233/150‑1 (MQ=255)
                                                                                                tAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaaa                                                       <  1:3030525/150‑1 (MQ=255)
                                                                                                tAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaaa                                                       >  2:4419590/1‑150 (MQ=255)
                                                                                                aaaCGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaaa                                                       >  2:2347336/2‑150 (MQ=255)
                                                                                                  aCGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACa                                                     <  1:2874787/150‑1 (MQ=255)
                                                                                                  aCGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACa                                                     <  1:3805082/150‑1 (MQ=255)
                                                                                                     ggAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGt                                                  <  1:2019387/150‑1 (MQ=255)
                                                                                                      gAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTc                                                 <  1:3248985/150‑1 (MQ=255)
                                                                                                      gAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTc                                                 <  1:520680/150‑1 (MQ=255)
                                                                                                             gCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAAt                                          <  2:4303694/150‑1 (MQ=255)
                                                                                                              cGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAgtcatc                                         >  2:823610/1‑145 (MQ=255)
                                                                                                              cGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATc                                         >  1:4345569/1‑150 (MQ=255)
                                                                                                               ggTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCg                                        <  2:4673749/150‑1 (MQ=255)
                                                                                                                    cGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACct                                   <  1:2282147/150‑1 (MQ=255)
                                                                                                                    cGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACct                                   <  2:1492903/150‑1 (MQ=255)
                                                                                                                       aGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGACTCGCACctgct                                >  2:2811218/1‑150 (MQ=255)
                                                                                                                       aGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctgct                                >  1:4530975/1‑150 (MQ=255)
                                                                                                                        gCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctgctg                               >  2:1943776/1‑150 (MQ=255)
                                                                                                                        gCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGGTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCCGCGCACACCAGACCAGAAACAGGTCGAGGAATCGCACctgctg                               >  1:1559025/1‑150 (MQ=255)
                                                                                                                          cGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGtt                             <  1:3311384/150‑1 (MQ=255)
                                                                                                                               cccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCt                        <  2:2371034/150‑1 (MQ=255)
                                                                                                                                ccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTa                       >  2:675564/1‑150 (MQ=255)
                                                                                                                                ccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTa                       <  2:1658963/150‑1 (MQ=255)
                                                                                                                                 cAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTac                      >  1:3421016/1‑150 (MQ=255)
                                                                                                                                   gCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTacac                    <  1:4487042/150‑1 (MQ=255)
                                                                                                                                     gTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCa                  >  1:1731779/1‑150 (MQ=255)
                                                                                                                                       cAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAAt                >  1:1104246/1‑150 (MQ=255)
                                                                                                                                       cAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAAt                <  1:2371894/150‑1 (MQ=255)
                                                                                                                                        aGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATc               <  1:1691203/150‑1 (MQ=255)
                                                                                                                                        aGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATc               <  1:4638843/150‑1 (MQ=255)
                                                                                                                                           cccGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGc            >  2:3478764/1‑150 (MQ=255)
                                                                                                                                           cccGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGc            >  1:4146594/1‑150 (MQ=255)
                                                                                                                                            ccGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCg           >  1:3267817/1‑150 (MQ=255)
                                                                                                                                             cGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCga          <  2:2766197/150‑1 (MQ=255)
                                                                                                                                             cGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCga          <  2:1631902/150‑1 (MQ=255)
                                                                                                                                              ggTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgag         >  2:1227312/1‑150 (MQ=255)
                                                                                                                                               gTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTACCGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgaga        <  2:3170932/150‑1 (MQ=255)
                                                                                                                                                tttGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGaa       <  2:2830643/150‑1 (MQ=255)
                                                                                                                                                    cGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACggg   >  2:628568/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808497‑2808795

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: