Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,194,194 C→A 8.0% S15* (TCA→TAA)  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,1940CA 8.0% 253.4 / 5.2 70S15* (TCA→TAA) lite14 prophage; cell death peptidase, inhibitor of T4 late gene expression
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base A (2/4);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACA  >  CP009273/1194052‑1194342
                                                                                                                                              |                                                                                                                                                    
aTAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCAttt                                                                                                                                               >  2:997940/1‑150 (MQ=255)
  aaTTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTaa                                                                                                                                             >  2:1513337/1‑150 (MQ=255)
  aaTTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTaa                                                                                                                                             >  1:2058590/1‑150 (MQ=255)
           cGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCAcc                                                                                                                                    <  2:2193985/150‑1 (MQ=255)
                ataGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAgaga                                                                                                                               >  1:4412508/1‑150 (MQ=255)
                      tCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCCCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGc                                                                                                                         >  2:3646143/1‑150 (MQ=255)
                       cGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCa                                                                                                                        <  2:1690656/150‑1 (MQ=255)
                        gAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                                                                                       <  1:1740245/150‑1 (MQ=255)
                        gAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                                                                                       <  1:1741049/150‑1 (MQ=255)
                        gAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCaa                                                                                                                       <  1:204740/150‑1 (MQ=255)
                             tAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCt                                                                                                                  >  2:4462419/1‑150 (MQ=255)
                              aGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCtt                                                                                                                 <  1:378411/150‑1 (MQ=255)
                                  aaaaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTaaaa                                                                                                             >  2:1825393/1‑150 (MQ=255)
                                  aaaaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTaaaa                                                                                                             >  2:1825398/1‑150 (MQ=255)
                                   aaaCTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAAc                                                                                                            >  1:1556358/1‑150 (MQ=255)
                                      cTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGAt                                                                                                         >  1:875502/1‑150 (MQ=255)
                                        gAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAg                                                                                                       >  2:3823179/1‑150 (MQ=255)
                                          aGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGtt                                                                                                     <  1:889697/150‑1 (MQ=255)
                                          aGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGtt                                                                                                     <  1:4309421/150‑1 (MQ=255)
                                               cATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTgacga                                                                                                <  2:1689708/150‑1 (MQ=255)
                                                  cTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACaa                                                                                             <  2:768661/150‑1 (MQ=255)
                                                  cTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACaa                                                                                             <  2:768638/150‑1 (MQ=255)
                                                  cTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACaa                                                                                             <  2:768330/150‑1 (MQ=255)
                                                    tATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACaaaa                                                                                           <  2:2818788/150‑1 (MQ=255)
                                                          gTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATta                                                                                     >  1:2063341/1‑150 (MQ=255)
                                                          gTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATta                                                                                     >  1:4309158/1‑150 (MQ=255)
                                                             tataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATaa                                                                                  >  2:3438519/1‑150 (MQ=255)
                                                             tataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATaa                                                                                  >  2:3438530/1‑150 (MQ=255)
                                                             tataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATaa                                                                                  >  2:3439229/1‑150 (MQ=255)
                                                              ataCTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAAt                                                                                 <  1:2898280/150‑1 (MQ=255)
                                                                actacCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAAttt                                                                               >  2:2490768/1‑150 (MQ=255)
                                                                  tacCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCa                                                                             <  2:271296/150‑1 (MQ=255)
                                                                              caaTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGaa                                                                 <  1:1449800/150‑1 (MQ=255)
                                                                               aaTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAAc                                                                <  1:1339262/150‑1 (MQ=255)
                                                                                aTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAAcc                                                               <  1:1131306/150‑1 (MQ=255)
                                                                                aTGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAAcc                                                               <  1:1128184/150‑1 (MQ=255)
                                                                                 tGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCt                                                              >  1:1008413/1‑150 (MQ=255)
                                                                                 tGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCt                                                              >  1:1002628/1‑150 (MQ=255)
                                                                                 tGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCt                                                              >  1:613010/1‑150 (MQ=255)
                                                                                      agGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGttt                                                         >  2:3680095/1‑150 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:1046632/150‑1 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:3322900/150‑1 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:4194861/150‑1 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:4194859/150‑1 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:4194824/150‑1 (MQ=255)
                                                                                       ggTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGtttt                                                        <  1:3322922/150‑1 (MQ=255)
                                                                                        gTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTa                                                       >  2:3363493/1‑150 (MQ=255)
                                                                                               aaaaaTGCGCTCACAAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCg                                                >  2:3206994/1‑150 (MQ=255)
                                                                                                aaaaTGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGa                                               <  1:508461/150‑1 (MQ=255)
                                                                                                aaaaTGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGa                                               <  1:2342821/150‑1 (MQ=255)
                                                                                                    tGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCa                                           <  1:1964168/150‑1 (MQ=255)
                                                                                                       gcTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGga                                        >  2:2485051/1‑150 (MQ=255)
                                                                                                                 ttGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGGGTCAGGAAGAATGagag                              <  1:785498/150‑1 (MQ=255)
                                                                                                                 ttGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGagag                              <  1:2372769/150‑1 (MQ=255)
                                                                                                                 ttGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGagag                              <  1:2371860/150‑1 (MQ=255)
                                                                                                                  tGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGt                             >  1:1807849/1‑150 (MQ=255)
                                                                                                                     cATCTTTTCTCAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGtaat                          <  2:521828/150‑1 (MQ=255)
                                                                                                                     cATCTTTTCTCAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGtaat                          <  2:514884/150‑1 (MQ=255)
                                                                                                                      aTCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGtaata                         >  2:1677432/1‑150 (MQ=255)
                                                                                                                      aTCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGtaata                         >  1:2671747/1‑150 (MQ=255)
                                                                                                                         ttttCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAAtt                      <  1:2259643/150‑1 (MQ=255)
                                                                                                                         ttttCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAAtt                      <  1:2259297/150‑1 (MQ=255)
                                                                                                                            tctcAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTcac                   <  2:4412508/150‑1 (MQ=255)
                                                                                                                             ctcAGCAATTAATTCATAACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTcaca                  >  2:393724/1‑150 (MQ=255)
                                                                                                                               cAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGATCCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATg                <  1:3721279/150‑1 (MQ=255)
                                                                                                                                 gCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGaa              >  2:4187420/1‑150 (MQ=255)
                                                                                                                                   aattaattCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAAt            <  1:2101172/150‑1 (MQ=255)
                                                                                                                                   aattaattCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAAt            <  1:4349988/150‑1 (MQ=255)
                                                                                                                                       aattCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGtt        >  2:4176320/1‑150 (MQ=255)
                                                                                                                                         ttCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCt      <  2:2891208/150‑1 (MQ=255)
                                                                                                                                          tcatcaCCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTa     <  2:2671747/150‑1 (MQ=255)
                                                                                                                                             tcaCCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACa  <  1:2040571/150‑1 (MQ=255)
                                                                                                                                             tcaCCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACa  <  1:495151/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATAATTCATTACGCATATAGTATCGAACATAGAAAAAACTGAAGATTCATCTTATTTTGTATATACTACCTAGCCCAACAATGTAGAGGTTAACGAAAAATGCGCTCACCAATTTGTCATCTTTTCTCAGCAATTAATTCATCACCATTTAAGATTGCACCAGAGAAGGAGCAAGATCTTAAAACGATAGTTGACGACAAAAAAATTATAATTTCAGTTGTGAGTGAACCTGGTTTTAATATCCGAGTCAGGAAGAATGAGAGTAATAATTCACATGAAATAGTTCTAACA  >  CP009273/1194052‑1194342

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: