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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,840,3770CA 8.5% 744.4 / 29.0 210D281Y (GAT→TAT) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base C (112/80);  new base A (10/8);  total (122/88)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 214 total aligned reads displayed.

CGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGC  >  CP009273/1840231‑1840524
                                                                                                                                                  |                                                                                                                                                   
cgGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCg                                                                                                                                                    <  2:3498220/148‑1 (MQ=255)
  gATAAGTGATTAGCTTGTGCGTTTCGTACCGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCACAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGtt                                                                                                                                                  <  1:1823236/148‑1 (MQ=255)
   aTAAGTTATTAGCTTGTGCGTTTCGTAAAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAATCGAAAAAGTCAGCGTCGCGGATTCTGATTCCCGTTTATCGTta                                                                                                                                                 <  2:3219510/148‑1 (MQ=255)
   aTAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTta                                                                                                                                                 >  1:1292146/1‑148 (MQ=255)
    tAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTtat                                                                                                                                                <  2:1022470/148‑1 (MQ=255)
      aGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAg                                                                                                                                              <  2:432299/148‑1 (MQ=255)
       gTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGc                                                                                                                                             <  1:1896524/148‑1 (MQ=255)
       gTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGc                                                                                                                                             <  1:1896092/148‑1 (MQ=255)
         gATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTg                                                                                                                                           <  2:2274622/148‑1 (MQ=255)
          aTTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTgg                                                                                                                                          <  2:1722510/148‑1 (MQ=255)
           ttAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGt                                                                                                                                         >  1:2601020/1‑148 (MQ=255)
           ttAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCAAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGt                                                                                                                                         >  1:2550567/1‑148 (MQ=255)
             aGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGa                                                                                                                                       <  2:3104137/148‑1 (MQ=255)
              gCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGAc                                                                                                                                      >  2:1371739/1‑148 (MQ=255)
              gCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGAc                                                                                                                                      >  2:2366610/1‑148 (MQ=255)
                      gTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCg                                                                                                                              <  1:2253657/148‑1 (MQ=255)
                        ttCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGt                                                                                                                            >  1:2939474/1‑148 (MQ=255)
                           gTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGa                                                                                                                         <  2:1225007/148‑1 (MQ=255)
                             aCAGTTTCTTGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGAcc                                                                                                                       >  1:2346057/1‑148 (MQ=255)
                             aCAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGAcc                                                                                                                       >  1:70899/1‑148 (MQ=255)
                              cAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCa                                                                                                                      >  1:493341/1‑148 (MQ=255)
                               aGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCAc                                                                                                                     <  1:3472398/148‑1 (MQ=255)
                                 tttCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTa                                                                                                                   >  1:1100435/1‑148 (MQ=255)
                                 tttCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTa                                                                                                                   >  2:2675087/1‑148 (MQ=255)
                                 tttCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTa                                                                                                                   >  2:2795394/1‑148 (MQ=255)
                                  ttCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTaa                                                                                                                  >  2:1336122/1‑148 (MQ=255)
                                            tatCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAAc                                                                                                        >  1:1692230/1‑148 (MQ=255)
                                            tatCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAAc                                                                                                        >  1:1946723/1‑148 (MQ=255)
                                            tatCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAAc                                                                                                        <  1:43077/148‑1 (MQ=255)
                                             atCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAAcc                                                                                                       >  2:2091487/1‑148 (MQ=255)
                                              tCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACca                                                                                                      >  1:2477718/1‑148 (MQ=255)
                                               cAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACcac                                                                                                     >  1:3031243/1‑148 (MQ=255)
                                                  gCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATg                                                                                                  >  1:977389/1‑148 (MQ=255)
                                                  gCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATg                                                                                                  >  1:846343/1‑148 (MQ=255)
                                                   cACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGc                                                                                                 >  2:1268364/1‑148 (MQ=255)
                                                     cGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTc                                                                                               >  2:3161447/1‑148 (MQ=255)
                                                     cGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTc                                                                                               <  1:211749/148‑1 (MQ=255)
                                                     cGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTc                                                                                               <  1:482941/148‑1 (MQ=255)
                                                        tCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGc                                                                                            <  2:558265/148‑1 (MQ=255)
                                                            cgcACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGt                                                                                        <  1:262875/148‑1 (MQ=255)
                                                              cACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTgg                                                                                      <  2:1884958/148‑1 (MQ=255)
                                                              cACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTgg                                                                                      >  2:859818/1‑148 (MQ=255)
                                                               aCTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGa                                                                                     <  1:3185979/148‑1 (MQ=255)
                                                                 tCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACg                                                                                   <  2:1781823/148‑1 (MQ=255)
                                                                  cAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGa                                                                                  >  1:475584/1‑148 (MQ=255)
                                                                  cAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGa                                                                                  >  2:331249/1‑148 (MQ=255)
                                                                  cAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGa                                                                                  <  2:2857431/148‑1 (MQ=255)
                                                                  cAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGa                                                                                  >  1:492772/1‑148 (MQ=255)
                                                                      ccAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTa                                                                              <  1:1340372/148‑1 (MQ=255)
                                                                       cAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTaa                                                                             >  1:1939734/1‑148 (MQ=255)
                                                                       cAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTaa                                                                             >  2:2528083/1‑148 (MQ=255)
                                                                         aaaCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGGTTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACa                                                                           >  2:3417196/1‑148 (MQ=255)
                                                                          aaCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACaa                                                                          <  1:131435/148‑1 (MQ=255)
                                                                          aaCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACaa                                                                          <  1:1199388/148‑1 (MQ=255)
                                                                           aCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAg                                                                         >  1:926974/1‑148 (MQ=255)
                                                                           aCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAg                                                                         >  2:1828067/1‑148 (MQ=255)
                                                                           aCGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAg                                                                         >  2:2382989/1‑148 (MQ=255)
                                                                                 tGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCct                                                                   <  2:1985813/148‑1 (MQ=255)
                                                                                    cGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCcttct                                                                >  1:633315/1‑148 (MQ=255)
                                                                                     gCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCcttctt                                                               <  2:70899/148‑1 (MQ=255)
                                                                                          aaTCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATcc                                                          <  2:654515/148‑1 (MQ=255)
                                                                                          aaTCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCACCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATcc                                                          >  1:3174702/1‑148 (MQ=255)
                                                                                               gCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTa                                                     <  2:533005/148‑1 (MQ=255)
                                                                                               gCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTa                                                     >  1:3342768/1‑148 (MQ=255)
                                                                                               gCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTa                                                     <  1:2966185/148‑1 (MQ=255)
                                                                                                cAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTAc                                                    >  2:2920802/1‑148 (MQ=255)
                                                                                                cAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTAc                                                    >  2:862122/1‑148 (MQ=255)
                                                                                                cAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTAc                                                    >  2:2008944/1‑148 (MQ=255)
                                                                                                      tGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTca                                              <  2:1420203/148‑1 (MQ=255)
                                                                                                       gAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTcac                                             >  2:2974009/1‑148 (MQ=255)
                                                                                                       gAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTcac                                             >  1:1358981/1‑148 (MQ=255)
                                                                                                       gAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTcac                                             >  2:1201969/1‑148 (MQ=255)
                                                                                                              aaaaaGGCAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTc                                      >  2:215342/1‑148 (MQ=255)
                                                                                                               aaaaGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACACGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCg                                     >  1:982862/1‑148 (MQ=255)
                                                                                                                 aaGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCt                                   <  2:3031243/148‑1 (MQ=255)
                                                                                                                  aGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTc                                  >  1:2322878/1‑148 (MQ=255)
                                                                                                                   ggCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCg                                 >  2:2082434/1‑148 (MQ=255)
                                                                                                                    gCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCgg                                <  2:2475153/148‑1 (MQ=255)
                                                                                                                     cAGCGGCGCGGATTCTGATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGt                               >  2:2666786/1‑148 (MQ=255)
                                                                                                                            gcgGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGa                        <  1:2398664/148‑1 (MQ=255)
                                                                                                                             cgGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGAt                       <  1:1222130/148‑1 (MQ=255)
                                                                                                                              ggATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATa                      >  1:158082/1‑148 (MQ=255)
                                                                                                                                aTTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGc                    >  2:826855/1‑148 (MQ=255)
                                                                                                                                  tCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCgg                  >  1:1015535/1‑148 (MQ=255)
                                                                                                                                     gATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTaa               <  1:2565864/148‑1 (MQ=255)
                                                                                                                                     gATTCCCGTTTATAGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTaa               <  1:1259889/148‑1 (MQ=255)
                                                                                                                                      aTTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTaaa              >  2:2935698/1‑148 (MQ=255)
                                                                                                                                      aTTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTaaa              >  1:3051095/1‑148 (MQ=255)
                                                                                                                                       ttCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAAc             <  1:2641795/148‑1 (MQ=255)
                                                                                                                                       ttCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAAc             <  1:2640842/148‑1 (MQ=255)
                                                                                                                                        tCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAAcc            >  2:974443/1‑148 (MQ=255)
                                                                                                                                         cccGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACCTGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCg           <  2:1925327/148‑1 (MQ=255)
                                                                                                                                           cGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCt         <  2:1536890/148‑1 (MQ=255)
                                                                                                                                              ttATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCAt      <  2:201803/148‑1 (MQ=255)
                                                                                                                                              ttATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCAt      <  2:1100435/148‑1 (MQ=255)
                                                                                                                                              ttATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCAt      <  1:3125285/148‑1 (MQ=255)
                                                                                                                                              ttATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCAt      >  1:1755077/1‑148 (MQ=255)
                                                                                                                                               tATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATc     <  2:1878392/148‑1 (MQ=255)
                                                                                                                                                 tCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCgg   >  2:1979827/1‑148 (MQ=255)
                                                                                                                                                  cGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGc  <  2:1681469/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAACGTTTTGCCGCTTCAATCTGCAACGCTGAAAGCGAAAAAGGCAGCGGCGCGGATTCTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGC  >  CP009273/1840231‑1840524

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: ‑