Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,194,365 G→T 22.7% W72L (TGG→TTG)  lit → e14 prophage; cell death peptidase, inhibitor of T4 late gene expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,194,3650GT22.7% 740.3 / 196.8 317W72L (TGG→TTG) lite14 prophage; cell death peptidase, inhibitor of T4 late gene expression
Reads supporting (aligned to +/- strand):  ref base G (120/125);  new base T (40/32);  total (160/157)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 330 total aligned reads displayed.

GAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAAACAACAG  >  CP009273/1194216‑1194514
                                                                                                                                                     |                                                                                                                                                     
  aGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTggg                                                                                                                                                     >  1:958638/1‑150 (MQ=255)
  aGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATGTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTggg                                                                                                                                                     <  1:697758/150‑1 (MQ=255)
   ggAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGc                                                                                                                                                    >  1:2300801/1‑150 (MQ=255)
        aaGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCAtttt                                                                                                                                               >  2:4961600/1‑150 (MQ=255)
           aTCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCa                                                                                                                                            <  2:4851979/150‑1 (MQ=255)
            tCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCaa                                                                                                                                           >  2:149102/1‑150 (MQ=255)
              ttAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAAtt                                                                                                                                         <  2:3682841/150‑1 (MQ=255)
              ttAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAAtt                                                                                                                                         <  1:4048168/150‑1 (MQ=255)
                  aaCGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCtt                                                                                                                                     >  1:2814261/1‑150 (MQ=255)
                   aCGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCttt                                                                                                                                    >  1:1539161/1‑150 (MQ=255)
                   aCGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCttt                                                                                                                                    <  2:129470/150‑1 (MQ=255)
                    cGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCtttt                                                                                                                                   <  1:4231731/150‑1 (MQ=255)
                     gATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTg                                                                                                                                  >  1:2901196/1‑150 (MQ=255)
                       tAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTggg                                                                                                                                <  2:9752/150‑1 (MQ=255)
                       tAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTggg                                                                                                                                <  1:1050923/150‑1 (MQ=255)
                        aGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGt                                                                                                                               >  1:2677964/1‑150 (MQ=255)
                        aGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGt                                                                                                                               <  1:1600069/150‑1 (MQ=255)
                          ttGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGttt                                                                                                                             <  1:2150941/150‑1 (MQ=255)
                          ttGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGttt                                                                                                                             <  2:3372501/150‑1 (MQ=255)
                          ttGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGttt                                                                                                                             >  1:4708584/1‑150 (MQ=255)
                                    aaaaaTTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAaga                                                                                                                   <  2:72659/150‑1 (MQ=255)
                                    aaaaaTTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAaga                                                                                                                   <  2:232977/150‑1 (MQ=255)
                                      aaaTTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGt                                                                                                                 >  2:2513928/1‑150 (MQ=255)
                                      aaaTTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGt                                                                                                                 >  2:5105367/1‑150 (MQ=255)
                                       aaTTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTa                                                                                                                >  2:4254552/1‑150 (MQ=255)
                                        aTTTTAATTTCTGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGATTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTAc                                                                                                               <  1:1639544/150‑1 (MQ=255)
                                              aTTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTATTTGGTTTTTACGCAAGAGTACTCCaaa                                                                                                         >  2:1401480/1‑150 (MQ=255)
                                                   aGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGTGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCa                                                                                                    <  2:3830550/150‑1 (MQ=255)
                                                       gtgagtgaACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGACTATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGaaa                                                                                                >  2:2554609/1‑150 (MQ=255)
                                                         gagtgaACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGaaaaa                                                                                              >  1:5130789/1‑150 (MQ=255)
                                                          agtgaACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAaat                                                                                             >  1:1981732/1‑150 (MQ=255)
                                                             gaACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAaataat                                                                                          >  2:3063159/1‑150 (MQ=255)
                                                              aaCCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatg                                                                                         >  1:1442815/1‑150 (MQ=255)
                                                              aaCCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatg                                                                                         >  2:1066112/1‑150 (MQ=255)
                                                              aaCCTGGTTTGAATATCCGAGTCAGGAAGAATGAGAGTAATAAGTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatg                                                                                         <  2:5155061/150‑1 (MQ=255)
                                                               aCCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatga                                                                                        <  1:2182644/150‑1 (MQ=255)
                                                               aCCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatga                                                                                        >  2:2078627/1‑150 (MQ=255)
                                                                ccTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatgat                                                                                       <  2:1190878/150‑1 (MQ=255)
                                                                 cTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCCAGAGTACTCCAAATCTCAGAAAAATAattatg                                                                                      >  2:3702079/1‑150 (MQ=255)
                                                                 cTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAatgatg                                                                                      <  1:136694/150‑1 (MQ=255)
                                                                  tGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAAtgatga                                                                                     >  2:3935464/1‑150 (MQ=255)
                                                                  tGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAAtgatga                                                                                     <  2:1925893/150‑1 (MQ=255)
                                                                        tAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACtt                                                                               >  2:1751209/1‑150 (MQ=255)
                                                                         aaTATGCGAGTCATGAAGAATGAGAGTTATTATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACttt                                                                              <  2:4916772/150‑1 (MQ=255)
                                                                         aaTATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACttt                                                                              <  1:4144103/150‑1 (MQ=255)
                                                                         aaTATCCGAGTCAGGAAGAATGAGAGTAATAATTAACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACttt                                                                              <  2:330113/150‑1 (MQ=255)
                                                                          atatCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTg                                                                             >  2:2908940/1‑150 (MQ=255)
                                                                          atatCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTg                                                                             <  2:595568/150‑1 (MQ=255)
                                                                           tatCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGa                                                                            >  2:3750278/1‑150 (MQ=255)
                                                                                gAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTaa                                                                       <  2:2909985/150‑1 (MQ=255)
                                                                                 aGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAAc                                                                      <  1:983944/150‑1 (MQ=255)
                                                                                 aGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAAc                                                                      <  1:235761/150‑1 (MQ=255)
                                                                                 aGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAAc                                                                      >  2:2202675/1‑150 (MQ=255)
                                                                                 aGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAAc                                                                      <  2:2178711/150‑1 (MQ=255)
                                                                                  gTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACa                                                                     <  1:3069660/150‑1 (MQ=255)
                                                                                   tCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAg                                                                    <  2:5177997/150‑1 (MQ=255)
                                                                                    cAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAgg                                                                   >  2:4310538/1‑150 (MQ=255)
                                                                                         ataaTGAGAGTATTTATTCTCATGAAATATTTCTAACAGTAGCTTCACTTGAATATATTTTGGTATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGaaaaa                                                              <  2:4089287/148‑1 (MQ=255)
                                                                                         agaaTGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGaaaaa                                                              >  1:3602593/1‑150 (MQ=255)
                                                                                          gaaTGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGaaaaaa                                                             >  2:3806947/1‑150 (MQ=255)
                                                                                            aTGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATa                                                           <  2:648426/150‑1 (MQ=255)
                                                                                            aTGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATa                                                           <  2:313738/150‑1 (MQ=255)
                                                                                             tGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAg                                                          <  1:5042883/150‑1 (MQ=255)
                                                                                              gagaGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAgg                                                         >  2:1695218/1‑150 (MQ=255)
                                                                                               agagTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGc                                                        >  1:1801163/1‑150 (MQ=255)
                                                                                                gagTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCt                                                       <  1:2066089/150‑1 (MQ=255)
                                                                                                gagTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCt                                                       <  2:958638/150‑1 (MQ=255)
                                                                                                    aataatTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTaaa                                                   >  2:3045600/1‑150 (MQ=255)
                                                                                                    aataatTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTaaa                                                   >  2:564122/1‑150 (MQ=255)
                                                                                                     ataatTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTaaaa                                                  >  2:1937015/1‑150 (MQ=255)
                                                                                                      taatTCACATCAAATAGTTCTAACATTAGCTTCACTTGAATATATTTGGGCATTTTACAATTTCTTTTGTGTTTTTACGCAAGAGTACTCCAAATCTCATAAAAATAATGATGATCACTTTGATTTAACAGGAAAAAATAGGCTTaaaaa                                                 <  2:3242702/150‑1 (MQ=255)
                                                                                                       aatTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAg                                                <  1:2222016/150‑1 (MQ=255)
                                                                                                        atTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGt                                               >  1:5213311/1‑150 (MQ=255)
                                                                                                         ttCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTc                                              <  1:1300587/150‑1 (MQ=255)
                                                                                                          tCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCt                                             >  1:4781314/1‑150 (MQ=255)
                                                                                                           cacaTGAAATAGTTCTAACAGTAGCTTCACTTGAATATCTTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGCAAAAATAGGCTTAAAAAGTCtg                                            >  2:4483122/1‑150 (MQ=255)
                                                                                                           cacaTGAAATAGTTCTAACAGTAGCTTCAATTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCtg                                            <  2:4703372/150‑1 (MQ=255)
                                                                                                           cacaTGAAATAGGTCTAACAGTAGCTTCACTTGCATATATTTGGGCATTTTCCAATTTCATTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCtg                                            >  2:4553557/1‑150 (MQ=255)
                                                                                                               tGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCtgatga                                        >  1:2506328/1‑150 (MQ=255)
                                                                                                                     aGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTACAATTTCTTTTGGGTTTTTACGCAAGAGTACTCAAAATCTCAGAAAAATAATGATGAGAACTTTGATTTAACAGGAAAAAATAGGATTAAAAAGTCTGAAGAActtct                                  >  1:4589758/1‑150 (MQ=255)
                                                                                                                      gTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAActtctt                                 >  2:3488594/1‑150 (MQ=255)
                                                                                                                            aCAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTGTGAGTTAACAGGCAAAAATAGGCTTAACAAGTCTGCTGCACTACTTACATgg                           >  2:2742101/1‑150 (MQ=255)
                                                                                                                             cattAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATggg                          <  2:4207965/147‑1 (MQ=255)
                                                                                                                             cAGTAGCTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATggg                          <  1:1212941/150‑1 (MQ=255)
                                                                                                                             cAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATggg                          <  2:1367361/150‑1 (MQ=255)
                                                                                                                              aGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGc                         <  1:3800048/150‑1 (MQ=255)
                                                                                                                               gTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCa                        <  1:891340/150‑1 (MQ=255)
                                                                                                                                   cTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGa                    <  1:1949217/150‑1 (MQ=255)
                                                                                                                                   cTTCACTTGAATATATTTTGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGa                    >  1:1464476/1‑150 (MQ=255)
                                                                                                                                    ttCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGaa                   >  1:1186101/1‑150 (MQ=255)
                                                                                                                                          tGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACtt             <  1:1228586/150‑1 (MQ=255)
                                                                                                                                          tGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACtt             >  2:1417943/1‑150 (MQ=255)
                                                                                                                                           gAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTg            >  2:2601908/1‑150 (MQ=255)
                                                                                                                                           gAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTg            <  2:883187/150‑1 (MQ=255)
                                                                                                                                           gAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTg            <  1:3935464/150‑1 (MQ=255)
                                                                                                                                               atatTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCaaa        >  1:4855384/1‑150 (MQ=255)
                                                                                                                                                tatTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAaac       <  1:3515333/150‑1 (MQ=255)
                                                                                                                                                  tttGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAaacaa     >  1:1747436/1‑150 (MQ=255)
                                                                                                                                                    tGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCACGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAAacaaca   >  1:1338622/1‑150 (MQ=255)
                                                                                                                                                     gggCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAAACAACAg  >  1:667403/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACAGTAGCTTCACTTGAATATATTTGGGCATTTTCCAATTTCTTTTGGGTTTTTACGCAAGAGTACTCCAAATCTCAGAAAAATAATGATGAGCACTTTGATTTAACAGGAAAAAATAGGCTTAAAAAGTCTGATGAACTTCTTAAATGGGCAAGGAAAAACTTGCAAACAACAG  >  CP009273/1194216‑1194514

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: