Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 564,319 C→A 5.8% intergenic (+216/‑39) emrE → / → ybcK DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273564,3190CA 5.8% 973.2 / 20.3 265intergenic (+216/‑39)emrE/ybcKDLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base C (142/107);  new base A (14/2);  total (156/109)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 281 total aligned reads displayed.

TTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTAAATAGTG  >  CP009273/564172‑564466
                                                                                                                                                   |                                                                                                                                                   
       tGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCACACCACTAAgg                                                                                                                                              >  1:92451/1‑148 (MQ=255)
          cTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGttt                                                                                                                                           <  1:2757703/148‑1 (MQ=255)
              aTTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATaa                                                                                                                                       >  2:669791/1‑148 (MQ=255)
              aTTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAAAAAAAATTATATCCCACCACTAAGGTTTATaa                                                                                                                                       >  2:2919790/1‑148 (MQ=255)
               ttttttATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATaaa                                                                                                                                      >  2:1755953/1‑148 (MQ=255)
               ttttttATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATaaa                                                                                                                                      <  2:837881/148‑1 (MQ=255)
               ttttttATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACAACTAAGGTTTATaaa                                                                                                                                      >  2:3471606/1‑148 (MQ=255)
                tttttATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATaaaa                                                                                                                                     >  1:787282/1‑148 (MQ=255)
                   ttATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCa                                                                                                                                  >  2:1185449/1‑148 (MQ=255)
                     aTGTATTCTTTTGTGCTTTCAAGATTAGGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATa                                                                                                                                <  1:2112927/148‑1 (MQ=255)
                      tGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATAc                                                                                                                               >  1:227358/1‑148 (MQ=255)
                          ttCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTa                                                                                                                           >  1:1598035/1‑148 (MQ=255)
                          ttCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTa                                                                                                                           >  1:857539/1‑148 (MQ=255)
                          ttCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTa                                                                                                                           <  2:1376738/148‑1 (MQ=255)
                              tttGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAg                                                                                                                       >  1:1519374/1‑148 (MQ=255)
                               ttGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAgg                                                                                                                      >  2:1797218/1‑148 (MQ=255)
                                 gtgCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACAACTAAGGTTTATAAAAGCATACGTTAGCAGttg                                                                                                                    >  2:331111/1‑148 (MQ=255)
                                   gCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTc                                                                                                                  <  1:1316518/148‑1 (MQ=255)
                                     tttCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCAc                                                                                                                <  2:2454193/148‑1 (MQ=255)
                                      ttCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCAcc                                                                                                               <  1:2178056/148‑1 (MQ=255)
                                         aaGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATg                                                                                                            <  2:1009990/148‑1 (MQ=255)
                                               aTGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGaaaaaa                                                                                                      >  2:478383/1‑148 (MQ=255)
                                                tGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAg                                                                                                     >  1:2578114/1‑148 (MQ=255)
                                                tGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCACACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAg                                                                                                     >  2:3232747/1‑148 (MQ=255)
                                                 ggCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGc                                                                                                    <  2:156344/148‑1 (MQ=255)
                                                 ggCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGc                                                                                                    <  1:1099120/148‑1 (MQ=255)
                                                     tAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTTTAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATa                                                                                                <  1:1791632/148‑1 (MQ=255)
                                                     tAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATa                                                                                                >  2:2122649/1‑148 (MQ=255)
                                                      aaGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAg                                                                                               >  1:2329633/1‑148 (MQ=255)
                                                      aaGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAg                                                                                               >  1:133123/1‑148 (MQ=255)
                                                      aaGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAg                                                                                               <  1:3243997/148‑1 (MQ=255)
                                                       aGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACAACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGc                                                                                              <  1:3089720/148‑1 (MQ=255)
                                                        gAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCa                                                                                             <  2:2286401/148‑1 (MQ=255)
                                                        gAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCa                                                                                             >  1:1996352/1‑148 (MQ=255)
                                                           aaTTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCatat                                                                                          <  2:2089371/148‑1 (MQ=255)
                                                            aTTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCatata                                                                                         <  1:2950918/148‑1 (MQ=255)
                                                              tGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATg                                                                                       >  2:1518292/1‑148 (MQ=255)
                                                                  aTACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAt                                                                                   >  1:3139594/1‑148 (MQ=255)
                                                                  aTACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAt                                                                                   <  1:3179775/148‑1 (MQ=255)
                                                                  aTACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAt                                                                                   >  2:1251689/1‑148 (MQ=255)
                                                                  aTACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACAACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAt                                                                                   >  2:1211021/1‑148 (MQ=255)
                                                                    aCGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAttt                                                                                 <  1:2898591/148‑1 (MQ=255)
                                                                     cGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGAtttt                                                                                <  1:285426/148‑1 (MQ=255)
                                                                      gATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTtc                                                                               <  2:1982903/148‑1 (MQ=255)
                                                                        ttattGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCCTAGCATATATGCGATTTtcat                                                                             >  2:3395587/1‑148 (MQ=255)
                                                                        ttattGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTtcat                                                                             >  2:1444169/1‑148 (MQ=255)
                                                                                gtATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGt                                                                     >  2:3358832/1‑148 (MQ=255)
                                                                                    atTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGa                                                                 <  2:549822/148‑1 (MQ=255)
                                                                                    atTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGa                                                                 <  2:2340400/148‑1 (MQ=255)
                                                                                      tCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATg                                                               >  2:1093311/1‑148 (MQ=255)
                                                                                      tCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATg                                                               <  2:3413609/148‑1 (MQ=255)
                                                                                          gATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTg                                                           <  1:1237411/148‑1 (MQ=255)
                                                                                          gATAATGTGACCTTAATTGTATTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTg                                                           <  1:2366175/148‑1 (MQ=255)
                                                                                            tAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACAACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGc                                                         >  1:301460/1‑148 (MQ=255)
                                                                                                  gACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCa                                                   >  1:3565756/1‑148 (MQ=255)
                                                                                                   aCCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCAt                                                  <  1:2433297/148‑1 (MQ=255)
                                                                                                     cTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTa                                                <  2:3079862/148‑1 (MQ=255)
                                                                                                     cTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTa                                                <  2:694073/148‑1 (MQ=255)
                                                                                                     cTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTa                                                >  2:667291/1‑148 (MQ=255)
                                                                                                       tAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTaaa                                              <  2:2539230/148‑1 (MQ=255)
                                                                                                          ttGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCg                                           >  2:299059/1‑148 (MQ=255)
                                                                                                          ttGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCg                                           >  2:2028960/1‑148 (MQ=255)
                                                                                                           tGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGa                                          <  1:1663219/148‑1 (MQ=255)
                                                                                                                ttttAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACagag                                     <  1:1186555/148‑1 (MQ=255)
                                                                                                                ttttAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTAATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACagag                                     >  2:1875484/1‑148 (MQ=255)
                                                                                                                 tttAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAgaga                                    >  2:90908/1‑148 (MQ=255)
                                                                                                                  ttAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGagaa                                   <  1:1853118/148‑1 (MQ=255)
                                                                                                                    aaataaaaaataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGagaaga                                 >  1:3416362/1‑148 (MQ=255)
                                                                                                                    aaataaaaaataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGagaaga                                 >  1:1925744/1‑148 (MQ=255)
                                                                                                                      ataaaaaataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACt                               >  1:1398251/1‑148 (MQ=255)
                                                                                                                        aaaaaataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTa                             <  1:2121981/148‑1 (MQ=255)
                                                                                                                          aaaataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTAtt                           >  2:2660811/1‑148 (MQ=255)
                                                                                                                             ataaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCt                        >  2:3169767/1‑148 (MQ=255)
                                                                                                                              taaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTg                       >  2:12933/1‑148 (MQ=255)
                                                                                                                              taaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTg                       <  2:1269510/148‑1 (MQ=255)
                                                                                                                               aaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGa                      >  2:1496920/1‑148 (MQ=255)
                                                                                                                               aaaCAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGa                      >  2:1605307/1‑148 (MQ=255)
                                                                                                                                aaCAAAAATTATATCCCACAACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGaa                     <  2:1708469/148‑1 (MQ=255)
                                                                                                                                  cAAAAATTATATCCCACCACTAAGGTTTATAAAATCATACGTTTGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGTTGTCTGGCGACTCATTAAACAGACATAGAAGACTTATTGCTGAATg                   <  1:55227/148‑1 (MQ=255)
                                                                                                                                  cAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATg                   <  1:649690/148‑1 (MQ=255)
                                                                                                                                   aaaaaTTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATgg                  <  1:822698/148‑1 (MQ=255)
                                                                                                                                    aaaaTTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGt                 >  2:2834435/1‑148 (MQ=255)
                                                                                                                                       aTTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGACGACTTATTGCTGAATGGTTaa              >  2:3361066/1‑148 (MQ=255)
                                                                                                                                       aTTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTaa              >  1:280525/1‑148 (MQ=255)
                                                                                                                                        ttATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTaaa             <  2:783717/148‑1 (MQ=255)
                                                                                                                                        ttATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTaaa             >  2:2746750/1‑148 (MQ=255)
                                                                                                                                         tataTCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAg            >  1:304259/1‑148 (MQ=255)
                                                                                                                                         tataTCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAg            >  1:2046505/1‑148 (MQ=255)
                                                                                                                                         tataTCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAg            >  2:581176/1‑148 (MQ=255)
                                                                                                                                         tataTCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAg            >  1:1723873/1‑148 (MQ=255)
                                                                                                                                          atatCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAgg           >  2:37208/1‑148 (MQ=255)
                                                                                                                                            atCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTa         >  1:2699720/1‑148 (MQ=255)
                                                                                                                                            atCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTa         <  1:1131403/148‑1 (MQ=255)
                                                                                                                                             tCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTaa        >  2:609451/1‑148 (MQ=255)
                                                                                                                                             tCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTaa        <  2:2181086/148‑1 (MQ=255)
                                                                                                                                              cccACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTaaa       >  2:1944322/1‑148 (MQ=255)
                                                                                                                                              cccACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTaaa       >  1:3473442/1‑148 (MQ=255)
                                                                                                                                                caccacTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTAAATa     <  1:3227071/148‑1 (MQ=255)
                                                                                                                                                caccacTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTAAATa     >  2:18830/1‑148 (MQ=255)
                                                                                                                                                  ccacTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTAAATAGt   >  1:1344062/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAAATAAACAAAAATTATATCCCACCACTAAGGTTTATAAAAGCATACGTTAGCAGGTGTCACCATGAAAAAAGCCATAGCATATATGCGATTTTCATCACCAGGTCAGATGTCTGGCGACTCATTAAACCGACAGAGAAGACTTATTGCTGAATGGTTAAAGGTAAATAGTG  >  CP009273/564172‑564466

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: