Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,840,083 C→A 17.5% D379Y (GAT→TAT)  topB ← DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,840,0830CA17.5% 696.9 / 103.3 233D379Y (GAT→TAT) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base C (92/100);  new base A (24/17);  total (116/117)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCG  >  CP009273/1839934‑1840232
                                                                                                                                                     |                                                                                                                                                     
gTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATc                                                                                                                                                       <  1:742908/150‑1 (MQ=255)
   cAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAc                                                                                                                                                    >  2:3546815/1‑150 (MQ=255)
   cAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAc                                                                                                                                                    <  2:1700875/150‑1 (MQ=255)
   cAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAc                                                                                                                                                    <  2:3363886/150‑1 (MQ=255)
      ttCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCtt                                                                                                                                                 <  2:1503395/150‑1 (MQ=255)
      ttCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCtt                                                                                                                                                 >  1:3265060/1‑150 (MQ=255)
      tacGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGCTCGGAATGATGGCGTGGTGCGCATCGAccat                                                                                                                                                 >  1:3264289/3‑148 (MQ=255)
        cGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCtttt                                                                                                                                               >  1:584442/1‑150 (MQ=255)
          aTAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTg                                                                                                                                             <  2:107504/150‑1 (MQ=255)
          aTAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTg                                                                                                                                             <  2:2919370/150‑1 (MQ=255)
          aTAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTg                                                                                                                                             >  1:1783067/1‑150 (MQ=255)
            aaCACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTc                                                                                                                                           >  1:3639917/1‑150 (MQ=255)
            aaCACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTc                                                                                                                                           >  1:2354650/1‑150 (MQ=255)
             acacacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCa                                                                                                                                          <  2:2823522/150‑1 (MQ=255)
              cacacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCAt                                                                                                                                         <  2:1638217/150‑1 (MQ=255)
               acacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATc                                                                                                                                        <  1:3138591/150‑1 (MQ=255)
                cacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATcc                                                                                                                                       >  2:3985281/1‑150 (MQ=255)
                cacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATcc                                                                                                                                       >  1:2771628/1‑150 (MQ=255)
                 acTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATccc                                                                                                                                      <  2:339040/150‑1 (MQ=255)
                  cTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCa                                                                                                                                     <  1:2247308/150‑1 (MQ=255)
                  cTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCa                                                                                                                                     <  2:468532/150‑1 (MQ=255)
                     gCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACa                                                                                                                                  <  1:3413798/150‑1 (MQ=255)
                               gCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGt                                                                                                                        <  2:3172660/150‑1 (MQ=255)
                                cATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGta                                                                                                                       >  2:919137/1‑150 (MQ=255)
                                      ggCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTg                                                                                                                 <  1:4109894/150‑1 (MQ=255)
                                        cAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGa                                                                                                               <  1:2861560/150‑1 (MQ=255)
                                          gAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATc                                                                                                             <  2:1720544/150‑1 (MQ=255)
                                           aaTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCGGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCt                                                                                                            >  1:2087566/1‑150 (MQ=255)
                                           aaTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTCTATCTGGATCt                                                                                                            >  2:677158/1‑150 (MQ=255)
                                           aaTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCt                                                                                                            <  2:2239902/150‑1 (MQ=255)
                                            aTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTa                                                                                                           <  1:656400/150‑1 (MQ=255)
                                              tGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTAcc                                                                                                         >  1:4477328/1‑150 (MQ=255)
                                                cATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCAc                                                                                                       >  1:941330/1‑150 (MQ=255)
                                                cATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACAAc                                                                                                       >  1:3272033/1‑150 (MQ=255)
                                                   cAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTgg                                                                                                    >  1:117160/1‑150 (MQ=255)
                                                   cAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTgg                                                                                                    >  1:118956/1‑150 (MQ=255)
                                                   cAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGAGCATAGACCTTTTTGTCATCCCAACAGCGGTTTCGTATATCTGGATCTACCACTgg                                                                                                    >  1:3494291/1‑150 (MQ=255)
                                                     gATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCt                                                                                                  <  2:2322999/150‑1 (MQ=255)
                                                                gCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGa                                                                                       >  2:550453/1‑150 (MQ=255)
                                                                gCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGa                                                                                       >  2:3899350/1‑150 (MQ=255)
                                                                 cAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGAt                                                                                      <  1:1694213/150‑1 (MQ=255)
                                                                   aTCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                                                    >  2:2646177/1‑150 (MQ=255)
                                                                   aTCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                                                    <  1:4198121/150‑1 (MQ=255)
                                                                   aTCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                                                    <  2:261902/150‑1 (MQ=255)
                                                                   aTCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                                                    <  2:2527885/150‑1 (MQ=255)
                                                                    tCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCg                                                                                   <  2:245735/150‑1 (MQ=255)
                                                                    tCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCg                                                                                   >  1:1371301/1‑150 (MQ=255)
                                                                           ataGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATg                                                                            >  1:4311749/1‑150 (MQ=255)
                                                                            taGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACAGGCTGCGGCAACAGATCCGGTGCATGa                                                                           >  1:1052852/1‑150 (MQ=255)
                                                                              gACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAac                                                                         <  1:4428021/150‑1 (MQ=255)
                                                                                 cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACt                                                                      >  1:290944/1‑150 (MQ=255)
                                                                                  ttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTg                                                                     <  2:2771628/150‑1 (MQ=255)
                                                                                   tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGa                                                                    >  2:3849844/1‑150 (MQ=255)
                                                                                     gcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATg                                                                  >  1:1955487/1‑150 (MQ=255)
                                                                                          gTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAtt                                                             >  2:1636533/1‑150 (MQ=255)
                                                                                              tcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                         <  1:2002682/150‑1 (MQ=255)
                                                                                              tcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                         >  1:249460/1‑150 (MQ=255)
                                                                                               ccGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatca                                                        <  2:4153225/150‑1 (MQ=255)
                                                                                               ccGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatca                                                        >  1:1898554/1‑150 (MQ=255)
                                                                                               ccGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatca                                                        <  2:1877259/150‑1 (MQ=255)
                                                                                                   cAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACcgc                                                    >  1:739210/1‑150 (MQ=255)
                                                                                                     ggTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGt                                                  >  1:125231/1‑150 (MQ=255)
                                                                                                      gTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTg                                                 >  1:2514738/1‑150 (MQ=255)
                                                                                                       ttGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGTCATTCATCACCGCGTgg                                                <  2:3249413/150‑1 (MQ=255)
                                                                                                       ttGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTgg                                                <  1:3067955/150‑1 (MQ=255)
                                                                                                       ttGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTgg                                                <  2:4528026/150‑1 (MQ=255)
                                                                                                        tGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGc                                               >  2:2107062/1‑150 (MQ=255)
                                                                                                         gATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCg                                              >  1:3742965/1‑150 (MQ=255)
                                                                                                           tCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTc                                            <  1:3796230/150‑1 (MQ=255)
                                                                                                              cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCgg                                         <  1:3892962/150‑1 (MQ=255)
                                                                                                              cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCgg                                         >  2:1144028/1‑150 (MQ=255)
                                                                                                               aGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc                                        >  2:3755794/1‑150 (MQ=255)
                                                                                                               aGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc                                        >  1:1292308/1‑150 (MQ=255)
                                                                                                                gAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCa                                       >  2:1120590/1‑150 (MQ=255)
                                                                                                                  aCTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaa                                     <  1:1412415/150‑1 (MQ=255)
                                                                                                                  aCTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaa                                     >  1:1827487/1‑150 (MQ=255)
                                                                                                                   cTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                    >  1:4003163/1‑150 (MQ=255)
                                                                                                                    tCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAAt                                   <  2:4368636/150‑1 (MQ=255)
                                                                                                                     ccGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATg                                  <  1:566812/150‑1 (MQ=255)
                                                                                                                     ccGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATg                                  >  1:2571988/1‑150 (MQ=255)
                                                                                                                     ccGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTACATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATg                                  <  2:15936/150‑1 (MQ=255)
                                                                                                                     ccGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTACATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATg                                  <  2:122455/150‑1 (MQ=255)
                                                                                                                        tGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGttc                               <  2:1400911/150‑1 (MQ=255)
                                                                                                                          cGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGttctt                             <  1:280800/150‑1 (MQ=255)
                                                                                                                           ggTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGttcttc                            >  1:423784/1‑150 (MQ=255)
                                                                                                                              cGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTgg                         >  2:507046/1‑150 (MQ=255)
                                                                                                                              cGGAATGATGGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTgg                         >  2:1841055/1‑150 (MQ=255)
                                                                                                                               ggAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGc                        >  1:849117/1‑150 (MQ=255)
                                                                                                                                  atgatgGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCaaa                     >  1:3653264/1‑150 (MQ=255)
                                                                                                                                  atgatgGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCaaa                     >  2:3638966/1‑150 (MQ=255)
                                                                                                                                    gatgGCGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCACAATGTTCTTCTGGCAAATa                   >  1:3087846/1‑150 (MQ=255)
                                                                                                                                         cGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc              >  2:3364814/1‑150 (MQ=255)
                                                                                                                                         cGTGGTGCGCATAGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc              >  2:2428547/1‑150 (MQ=255)
                                                                                                                                           tggtgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCTTTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACaa            <  1:1066858/150‑1 (MQ=255)
                                                                                                                                           tggtgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACaa            >  2:1829999/1‑150 (MQ=255)
                                                                                                                                                 gcATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGaa      <  1:1541580/150‑1 (MQ=255)
                                                                                                                                                  cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAAc     >  1:2775921/1‑150 (MQ=255)
                                                                                                                                                  cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAAc     <  2:1999062/150‑1 (MQ=255)
                                                                                                                                                   aTCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAAcg    <  2:3232378/150‑1 (MQ=255)
                                                                                                                                                     cGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAAcgcg  >  1:48717/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCG  >  CP009273/1839934‑1840232

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: