Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,982 G→T 90.2% T226T (ACG→ACT topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,9820GT90.2% 376.4 / 34.9 150T226T (ACG→ACTtopADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base T (56/79);  total (63/87)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

TGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAA  >  CP009273/1325837‑1326126
                                                                                                                                                 |                                                                                                                                                
tGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAc                                                                                                                                                >  1:1113594/1‑148 (MQ=255)
    ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACg                                                                                                                                            >  2:1037131/1‑148 (MQ=255)
    ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACg                                                                                                                                            <  1:1851591/148‑1 (MQ=255)
     tCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGc                                                                                                                                           <  2:2784742/148‑1 (MQ=255)
     tCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGc                                                                                                                                           <  2:2784686/148‑1 (MQ=255)
        ccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCAt                                                                                                                                        >  2:2912298/1‑148 (MQ=255)
         cgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                       <  1:1686021/148‑1 (MQ=255)
           ctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTg                                                                                                                                     >  1:3018270/1‑148 (MQ=255)
           ctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTAGCCAGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg                                                                                                                                     >  1:1399077/1‑148 (MQ=255)
             gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                   >  1:1794447/1‑148 (MQ=255)
             gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                   >  1:1865810/1‑148 (MQ=255)
             gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGt                                                                                                                                   >  1:1794254/1‑148 (MQ=255)
              ctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTg                                                                                                                                  >  2:822864/1‑148 (MQ=255)
                  ggAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGc                                                                                                                              >  2:1180888/1‑148 (MQ=255)
                  ggAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGc                                                                                                                              >  2:1180887/1‑148 (MQ=255)
                  ggAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGc                                                                                                                              >  2:802145/1‑148 (MQ=255)
                   gAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCg                                                                                                                             <  1:3010776/148‑1 (MQ=255)
                       aaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                         <  2:2723314/148‑1 (MQ=255)
                       aaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCGGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                         >  2:565318/1‑148 (MQ=255)
                          aTCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGt                                                                                                                      <  1:1107098/148‑1 (MQ=255)
                          aTCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTATGGGAAGTAGATGCCAGCACTACCACGCCAACTGGTGAAGCGTTGGCGt                                                                                                                      >  1:1791661/1‑148 (MQ=255)
                              cTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACa                                                                                                                  <  2:1295360/148‑1 (MQ=255)
                               tCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAg                                                                                                                 >  1:2201828/1‑148 (MQ=255)
                               tCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAg                                                                                                                 <  2:743875/148‑1 (MQ=255)
                                cGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAgg                                                                                                                <  1:1503117/148‑1 (MQ=255)
                                 gTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGt                                                                                                               <  2:3009405/148‑1 (MQ=255)
                                        gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcat                                                                                                        >  2:1029301/1‑148 (MQ=255)
                                            gCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGa                                                                                                    <  2:775302/148‑1 (MQ=255)
                                             ccGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGaa                                                                                                   <  1:3061314/148‑1 (MQ=255)
                                                gTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGa                                                                                                <  1:666504/148‑1 (MQ=255)
                                                gTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGa                                                                                                >  1:2704262/1‑148 (MQ=255)
                                                 tCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTAGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAc                                                                                               >  1:2078457/1‑148 (MQ=255)
                                                  cGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACa                                                                                              >  2:2856407/1‑148 (MQ=255)
                                                    tgtgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACaaa                                                                                            <  2:3301487/148‑1 (MQ=255)
                                                    tgtgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACaaa                                                                                            <  2:945388/148‑1 (MQ=255)
                                                       gCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAAcgg                                                                                         <  2:2081189/148‑3 (MQ=255)
                                                       gCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAAcgg                                                                                         <  2:2081179/148‑3 (MQ=255)
                                                         aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGtt                                                                                       <  2:3018270/148‑1 (MQ=255)
                                                         aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGtt                                                                                       <  2:2858203/148‑1 (MQ=255)
                                                         aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGtt                                                                                       <  1:1013855/148‑1 (MQ=255)
                                                          gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtc                                                                                      >  2:438540/1‑148 (MQ=255)
                                                          gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtc                                                                                      >  2:438761/1‑148 (MQ=255)
                                                             ggTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgt                                                                                   <  1:1885148/148‑1 (MQ=255)
                                                                     ttCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                           >  2:1950971/1‑148 (MQ=255)
                                                                     ttCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                           <  1:438761/148‑1 (MQ=255)
                                                                     ttCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                           <  1:2053235/148‑1 (MQ=255)
                                                                        gCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaacaa                                                                        <  2:3072156/148‑1 (MQ=255)
                                                                         ccTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACGGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACaaa                                                                       <  2:1773131/148‑1 (MQ=255)
                                                                            ggtggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGaa                                                                    >  1:830693/1‑148 (MQ=255)
                                                                             gtggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAAc                                                                   >  2:2604287/1‑148 (MQ=255)
                                                                                gtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaaa                                                                <  2:3274885/148‑1 (MQ=255)
                                                                                gtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaaa                                                                <  2:2479057/148‑1 (MQ=255)
                                                                                      cgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAgg                                                          <  1:57925/148‑1 (MQ=255)
                                                                                      cgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAgg                                                          <  2:1000419/148‑1 (MQ=255)
                                                                                        tgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCt                                                        >  2:3087372/1‑148 (MQ=255)
                                                                                         gagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTg                                                       <  1:2842134/148‑1 (MQ=255)
                                                                                            cgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCgg                                                    >  1:1395639/1‑148 (MQ=255)
                                                                                              tgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTa                                                  <  1:2274276/148‑1 (MQ=255)
                                                                                               gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTaa                                                 >  1:412995/1‑148 (MQ=255)
                                                                                                aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                >  2:347355/1‑148 (MQ=255)
                                                                                                aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                >  2:347348/1‑148 (MQ=255)
                                                                                                aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                >  1:2712266/1‑148 (MQ=255)
                                                                                                aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                <  2:520062/148‑1 (MQ=255)
                                                                                                 aaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                               <  1:762749/148‑1 (MQ=255)
                                                                                                          gTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaa                                      <  2:97328/148‑1 (MQ=255)
                                                                                                           ttCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTCGTGAAGCTTTGGCGTTACAGGTGACTCATAAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAAAAAACTAAGGCTGCGCTAAGTCTGCTGGaaa                                     >  2:1282001/1‑148 (MQ=255)
                                                                                                            tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaa                                    >  2:3103981/1‑148 (MQ=255)
                                                                                                            tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaa                                    <  2:1694691/148‑1 (MQ=255)
                                                                                                            tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaa                                    >  2:1540450/1‑148 (MQ=255)
                                                                                                             cGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaaa                                   <  1:2977354/148‑1 (MQ=255)
                                                                                                              gTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAg                                  <  2:1794254/148‑1 (MQ=255)
                                                                                                              gTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAg                                  <  2:2943907/148‑1 (MQ=255)
                                                                                                              gTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAg                                  <  2:1865810/148‑1 (MQ=255)
                                                                                                                gCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcg                                <  2:1531738/148‑1 (MQ=255)
                                                                                                                  cGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgcg                              <  1:1029301/148‑1 (MQ=255)
                                                                                                                   ggAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgcgc                             <  1:2231171/148‑1 (MQ=255)
                                                                                                                    gAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCt                            <  1:412847/148‑1 (MQ=255)
                                                                                                                     aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                           <  2:2704475/148‑1 (MQ=255)
                                                                                                                     aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                           >  1:796954/1‑148 (MQ=255)
                                                                                                                     aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                           >  1:2880779/1‑148 (MQ=255)
                                                                                                                      agagTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTAc                          <  2:2540433/148‑1 (MQ=255)
                                                                                                                        agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAg                        <  1:226672/148‑1 (MQ=255)
                                                                                                                        agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAg                        <  1:227244/148‑1 (MQ=255)
                                                                                                                          ttCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCGACAGCg                      <  1:1409902/148‑1 (MQ=255)
                                                                                                                           tCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGt                     >  1:1997197/1‑148 (MQ=255)
                                                                                                                            cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTg                    <  1:3206595/148‑1 (MQ=255)
                                                                                                                               ggAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTg                 >  2:1279714/1‑148 (MQ=255)
                                                                                                                                 aaGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGa               >  1:2568941/1‑148 (MQ=255)
                                                                                                                                  aGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa              <  1:802145/148‑1 (MQ=255)
                                                                                                                                   gTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAAc             >  1:2649488/1‑148 (MQ=255)
                                                                                                                                     cGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGt           <  2:2427648/148‑1 (MQ=255)
                                                                                                                                       aTGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGa         <  1:1180888/148‑1 (MQ=255)
                                                                                                                                       aTGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGa         <  1:1180887/148‑1 (MQ=255)
                                                                                                                                         gCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAg       >  2:3139688/1‑148 (MQ=255)
                                                                                                                                         gCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAg       <  1:2942613/148‑1 (MQ=255)
                                                                                                                                         gCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAg       <  1:1295581/148‑1 (MQ=255)
                                                                                                                                         gCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAg       >  2:1719539/1‑148 (MQ=255)
                                                                                                                                          ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGa      <  1:1570656/148‑1 (MQ=255)
                                                                                                                                              cACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa  >  1:2328530/1‑148 (MQ=255)
                                                                                                                                              cACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa  >  2:1373034/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAA  >  CP009273/1325837‑1326126

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: