Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,193,516 T→G 52.4% N60H (AAT→CAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,193,5160TG52.4% 43.0 / 126.7 90N60H (AAT→CAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/22);  new base G (19/28);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATTAT  >  CP009273/1193372‑1193662
                                                                                                                                                |                                                                                                                                                  
aCGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATgaataa                                                                                                                                               <  2:967572/150‑1 (MQ=255)
 cGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATaaa                                                                                                                                              <  1:1967886/150‑1 (MQ=255)
 cGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATaaa                                                                                                                                              <  1:2039881/150‑1 (MQ=255)
 cGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATaaa                                                                                                                                              <  1:2834224/150‑1 (MQ=255)
  ggTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAg                                                                                                                                             >  2:706676/1‑150 (MQ=255)
   gTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGc                                                                                                                                            >  1:1623303/1‑150 (MQ=255)
   gTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGc                                                                                                                                            >  1:3811952/1‑150 (MQ=255)
    ttttATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGct                                                                                                                                           <  2:2319571/150‑1 (MQ=255)
      ttATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGctct                                                                                                                                         >  2:1904524/1‑150 (MQ=255)
           aTTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGa                                                                                                                                    >  2:3669767/1‑150 (MQ=255)
           aTTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGa                                                                                                                                    >  2:1478422/1‑150 (MQ=255)
             tAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATc                                                                                                                                  <  1:3855692/150‑1 (MQ=255)
              aGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCt                                                                                                                                 <  1:92676/150‑1 (MQ=255)
              aGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCt                                                                                                                                 <  1:92688/150‑1 (MQ=255)
              aGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCt                                                                                                                                 >  1:2524333/1‑150 (MQ=255)
               gATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCtt                                                                                                                                <  1:964663/150‑1 (MQ=255)
                   cTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAAt                                                                                                                            <  1:2512390/150‑1 (MQ=255)
                    ttATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATa                                                                                                                           <  2:3937136/150‑1 (MQ=255)
                     tATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAg                                                                                                                          >  1:2561051/1‑150 (MQ=255)
                     tATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAg                                                                                                                          <  1:4277313/150‑1 (MQ=255)
                         ccAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTgg                                                                                                                      >  1:1875852/1‑150 (MQ=255)
                              gATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACAt                                                                                                                 >  1:3841275/1‑150 (MQ=255)
                                  tCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACtt                                                                                                             <  2:4065826/150‑1 (MQ=255)
                                  tCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACtt                                                                                                             <  2:2135337/150‑1 (MQ=255)
                                  tCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACtt                                                                                                             <  2:2135368/150‑1 (MQ=255)
                                  tCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACtt                                                                                                             <  2:2135363/150‑1 (MQ=255)
                                   cGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACttg                                                                                                            <  2:3704010/150‑1 (MQ=255)
                                    gCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACttgg                                                                                                           <  2:2711396/150‑1 (MQ=255)
                                    gCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACttgg                                                                                                           >  2:3748291/1‑150 (MQ=255)
                                       cccGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACttggttg                                                                                                        <  1:2776298/150‑1 (MQ=255)
                                          gAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTa                                                                                                     >  2:3049142/1‑150 (MQ=255)
                                            agtagtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAAt                                                                                                   <  1:3885188/150‑1 (MQ=255)
                                            agtagtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAAt                                                                                                   <  1:3885260/150‑1 (MQ=255)
                                               agtTTGGCTGCATTGTTATTTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGtt                                                                                                <  1:412810/150‑1 (MQ=255)
                                               agtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGtt                                                                                                <  1:2588514/150‑1 (MQ=255)
                                               agtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGtt                                                                                                >  1:4214313/1‑150 (MQ=255)
                                               agtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGtt                                                                                                >  1:4214496/1‑150 (MQ=255)
                                                gtTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGtta                                                                                               <  1:3599314/150‑1 (MQ=255)
                                                     gCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTaaa                                                                                          <  1:3748291/150‑1 (MQ=255)
                                                      cTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAg                                                                                         >  2:3313833/1‑150 (MQ=255)
                                                         cATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGa                                                                                      <  1:725356/150‑1 (MQ=255)
                                                          aTTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGAt                                                                                     >  2:2118382/1‑150 (MQ=255)
                                                             gTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTa                                                                                  <  1:2300739/150‑1 (MQ=255)
                                                                aTGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGAGAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTc                                                                               <  2:2053155/150‑1 (MQ=255)
                                                                   tAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTa                                                                            >  1:51430/1‑150 (MQ=255)
                                                                      gTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAAt                                                                         <  2:1505637/150‑1 (MQ=255)
                                                                      gTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAAt                                                                         <  2:3923143/150‑1 (MQ=255)
                                                                           aCGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTa                                                                    <  2:1323705/150‑1 (MQ=255)
                                                                            cGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAg                                                                   <  1:4326369/150‑1 (MQ=255)
                                                                                  gCTAATTAAATATGCAAAAATCGCAAAAGCACACGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACGTGGTTGGTAATGTGATTAAAGGGATGTATTCGTAAATAACTAGAAGTAt                                                             >  2:62482/1‑150 (MQ=255)
                                                                                   cTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATg                                                            <  2:3625915/150‑1 (MQ=255)
                                                                                   cTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATg                                                            <  2:3625860/150‑1 (MQ=255)
                                                                                     aaTTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCa                                                          <  2:2524333/150‑1 (MQ=255)
                                                                                     aaTTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCa                                                          <  2:2487757/150‑1 (MQ=255)
                                                                                       ttAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAat                                                        >  1:1507458/1‑150 (MQ=255)
                                                                                            tatGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATAtt                                                   <  1:706676/150‑1 (MQ=255)
                                                                                             atGCAAAAATCGCAAAGGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTa                                                  <  1:1763015/150‑1 (MQ=255)
                                                                                                cAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCg                                               >  2:3834199/1‑150 (MQ=255)
                                                                                                cAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCg                                               >  1:691395/1‑150 (MQ=255)
                                                                                                    aaTCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTAtt                                           <  1:3513424/150‑1 (MQ=255)
                                                                                                          aaaaGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTATTGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGaa                                     >  2:2302185/1‑150 (MQ=255)
                                                                                                          aaaaGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGaa                                     >  2:3960027/1‑150 (MQ=255)
                                                                                                           aaaGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGaaa                                    >  1:2902297/1‑150 (MQ=255)
                                                                                                           aaaGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGaaa                                    <  1:2746900/150‑1 (MQ=255)
                                                                                                           aaaGCACAAGGTAAGATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGaaa                                    <  2:2751849/150‑1 (MQ=255)
                                                                                                              gCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGtt                                 <  1:3669767/150‑1 (MQ=255)
                                                                                                              gCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGtt                                 <  1:1478422/150‑1 (MQ=255)
                                                                                                              gCACAAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGtt                                 <  2:691395/150‑1 (MQ=255)
                                                                                                               cacaAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCACTATATTAGCGTATTTAAGAAAGTTg                                >  1:1220462/1‑150 (MQ=255)
                                                                                                               cacaAGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTg                                >  1:2000360/1‑150 (MQ=255)
                                                                                                                 caAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCt                              >  2:3257642/1‑150 (MQ=255)
                                                                                                                  aaGGTAATATACCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCtt                             >  2:2978216/1‑150 (MQ=255)
                                                                                                                     gTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTAt                          <  2:3628997/150‑1 (MQ=255)
                                                                                                                     gTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTAt                          <  2:3628931/150‑1 (MQ=255)
                                                                                                                            aCCAAAAGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGttt                   >  1:1435371/1‑150 (MQ=255)
                                                                                                                            aCCAAAAGATAGTGCCCAAAAAATAAAGCTCTCATGATCTTAAATAGGTGGGACATACTTGGATGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATAAATTAGCGTATTTAAGAAAGTTGCTTTATGAAGttt                   >  2:3263305/1‑150 (MQ=255)
                                                                                                                                 aaGATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCa              >  1:1382947/1‑150 (MQ=255)
                                                                                                                                  aGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCaa             >  1:1944224/1‑150 (MQ=255)
                                                                                                                                   gATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAAc            >  1:3331371/1‑150 (MQ=255)
                                                                                                                                   gATAGTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAAc            >  1:3942454/1‑150 (MQ=255)
                                                                                                                                    aTAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACa           >  2:2999280/1‑150 (MQ=255)
                                                                                                                                    aTAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACa           >  2:2999250/1‑150 (MQ=255)
                                                                                                                                    aTAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACa           >  1:148124/1‑150 (MQ=255)
                                                                                                                                       gTGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACAttt        <  1:2107085/150‑1 (MQ=255)
                                                                                                                                        tGCCCAATGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTc       >  2:2545464/1‑150 (MQ=255)
                                                                                                                                          cccAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCAt     >  1:3903350/1‑150 (MQ=255)
                                                                                                                                            cAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCAtga   >  2:756933/1‑148 (MQ=255)
                                                                                                                                             aaTTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATtat  <  1:614037/150‑1 (MQ=255)
                                                                                                                                             aaTGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATtat  <  2:1798297/150‑1 (MQ=255)
                                                                                                                                             aaTGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATtat  <  2:1798128/150‑1 (MQ=255)
                                                                                                                                             aaTGAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATtat  <  2:1714935/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACGGTTTTATAATTAGATGCTTATCCCAAAGATATCGCACCCGAAGTAGTTTGGCTGCATTGTTATGTAGGTCTAACGCACCGCTAATTAAATATGCAAAAATCGCAAAAGCACAAGGTAATATACCAAAAGATAGTGCCCAATTAATAAAGCTCTCATGATCTTTAATAGGTGGGACATACTTGGTTGGTAATGTTATTAAAGGGATGTATTCGTAAATAACTAGAAGTATGCAATATATTAGCGTATTTAAGAAAGTTGCTTTATGAAGTTTAGGCAACATTTCATTAT  >  CP009273/1193372‑1193662

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: