Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 569,111 T→A 5.5% V95V (GTT→GTA rusA → DLP12 prophage; endonuclease RUS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273569,1110TA 5.5% 1150.3 / 39.5 325V95V (GTT→GTArusADLP12 prophage; endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base T (175/132);  new base A (13/5);  total (188/137)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 334 total aligned reads displayed.

ATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGT  >  CP009273/568962‑569260
                                                                                                                                                     |                                                                                                                                                     
aTGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGtt                                                                                                                                                       >  1:1105390/1‑150 (MQ=255)
aTGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGtt                                                                                                                                                       >  2:3136865/1‑150 (MQ=255)
  gCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGa                                                                                                                                                     >  2:4256092/1‑150 (MQ=255)
  gCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGGCGTTGa                                                                                                                                                     >  2:5222115/1‑150 (MQ=255)
     ggATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTAGATTa                                                                                                                                                  >  1:626264/1‑150 (MQ=255)
        tatCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcg                                                                                                                                               >  1:4653311/1‑150 (MQ=255)
         atCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcgc                                                                                                                                              >  1:1140917/1‑150 (MQ=255)
         atCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcgc                                                                                                                                              >  1:3981695/1‑150 (MQ=255)
            ggCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTAGATTACCGCGtt                                                                                                                                           >  2:3534916/1‑150 (MQ=255)
               cTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTtgtg                                                                                                                                        >  1:455951/1‑150 (MQ=255)
                tGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGa                                                                                                                                       >  1:3336028/1‑150 (MQ=255)
                  gCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAg                                                                                                                                     <  1:693442/150‑1 (MQ=255)
                  gCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAg                                                                                                                                     <  2:4897955/150‑1 (MQ=255)
                  gCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAg                                                                                                                                     <  2:3815483/150‑1 (MQ=255)
                    tATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGAt                                                                                                                                   >  1:1686572/1‑150 (MQ=255)
                    tATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGAt                                                                                                                                   >  2:3716731/1‑150 (MQ=255)
                      tGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGc                                                                                                                                 >  1:4274562/1‑150 (MQ=255)
                       gCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGcc                                                                                                                                >  2:3227577/1‑150 (MQ=255)
                                aaTCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCTTTGTGAAGATGCCTGTTACCaa                                                                                                                       >  2:3896394/1‑150 (MQ=255)
                                  tCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAg                                                                                                                     <  1:5028421/150‑1 (MQ=255)
                                     gCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTg                                                                                                                  <  2:4193678/150‑1 (MQ=255)
                                          gAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAgg                                                                                                             >  1:2246583/1‑150 (MQ=255)
                                           aGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGc                                                                                                            <  2:4945871/150‑1 (MQ=255)
                                           aGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGc                                                                                                            >  1:3742561/1‑150 (MQ=255)
                                              gCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTgg                                                                                                         >  1:3450431/1‑150 (MQ=255)
                                              gCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTgg                                                                                                         <  2:1731469/150‑1 (MQ=255)
                                              gCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTgg                                                                                                         >  2:2460201/1‑150 (MQ=255)
                                                 acaTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAAc                                                                                                      >  1:5626373/1‑150 (MQ=255)
                                                  caTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACt                                                                                                     <  2:1668284/150‑1 (MQ=255)
                                                    tGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGa                                                                                                   >  1:222714/1‑150 (MQ=255)
                                                       cGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCa                                                                                                >  2:1152638/1‑150 (MQ=255)
                                                           tcgccgtcgccgTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCAc                                                                                            >  2:2958252/1‑150 (MQ=255)
                                                            cgccgtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCAcc                                                                                           <  2:2333998/150‑1 (MQ=255)
                                                            cgccgtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCAcc                                                                                           <  2:2928176/150‑1 (MQ=255)
                                                             gccgtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTAGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCg                                                                                          >  1:4345236/1‑150 (MQ=255)
                                                               cgtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGaa                                                                                        <  1:2070079/150‑1 (MQ=255)
                                                                gtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGaaa                                                                                       >  1:2841959/1‑150 (MQ=255)
                                                                 tcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTCCCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAAt                                                                                      >  2:733947/1‑150 (MQ=255)
                                                                 tcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAAt                                                                                      <  2:2698464/150‑1 (MQ=255)
                                                                  cgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATg                                                                                     <  2:687428/150‑1 (MQ=255)
                                                                   gccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATgg                                                                                    >  2:4635457/1‑150 (MQ=255)
                                                                   gccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATgg                                                                                    >  1:5412327/1‑150 (MQ=255)
                                                                    ccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATggg                                                                                   <  1:2976766/150‑1 (MQ=255)
                                                                        gACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaat                                                                               >  1:879747/1‑150 (MQ=255)
                                                                         aCCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatg                                                                              >  1:3919390/1‑150 (MQ=255)
                                                                         aCCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatg                                                                              >  2:2597399/1‑150 (MQ=255)
                                                                          ccTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGGTGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatgt                                                                             >  1:4580165/1‑149 (MQ=255)
                                                                          ccTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACGAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatga                                                                             <  2:1117338/150‑1 (MQ=255)
                                                                           cTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatgaa                                                                            <  2:4026160/150‑1 (MQ=255)
                                                                              gATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAatgaatga                                                                         <  1:2681055/150‑1 (MQ=255)
                                                                               aTAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAatgat                                                                        >  2:537564/1‑150 (MQ=255)
                                                                                tAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAatgatg                                                                       >  2:5350860/1‑150 (MQ=255)
                                                                                 aaTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAACTGGGGAATGAATGATGt                                                                      >  1:4410638/1‑150 (MQ=255)
                                                                                 aaTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCGGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGt                                                                      >  2:83765/1‑150 (MQ=255)
                                                                                  aTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGtt                                                                     >  1:4151630/1‑150 (MQ=255)
                                                                                   tCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGttt                                                                    >  2:1323847/1‑150 (MQ=255)
                                                                                    cTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTg                                                                   <  2:2108461/150‑1 (MQ=255)
                                                                                     tGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGa                                                                  <  1:954220/150‑1 (MQ=255)
                                                                                      gCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAg                                                                 <  1:1868345/150‑1 (MQ=255)
                                                                                      gCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAg                                                                 <  1:1169588/150‑1 (MQ=255)
                                                                                         aaaaGCCGCTTTTGACGCACTCACTAAAGAAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACAGCTTTGTGAAGATGCATGTTACCAAAGTTGGGAGGCTGGAACTGACCATAACCGAAATGGGGAATGAATGATGTTTGAGttt                                                              >  1:4775894/1‑150 (MQ=255)
                                                                                          aaaGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTa                                                             >  2:1462038/1‑150 (MQ=255)
                                                                                             gCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAata                                                          <  1:5091718/150‑1 (MQ=255)
                                                                                              ccGCTTTTGACACACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAatat                                                         <  2:3616976/150‑1 (MQ=255)
                                                                                                 cTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAACTGGGGAATGAATGATGTTTGAGTTTAATATGGc                                                      >  2:5300926/1‑150 (MQ=255)
                                                                                                 cTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGc                                                      >  2:337394/1‑150 (MQ=255)
                                                                                                       aCGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAAct                                                <  1:3851314/150‑1 (MQ=255)
                                                                                                       aCGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAAct                                                <  1:1039705/150‑1 (MQ=255)
                                                                                                          cactcactAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTAGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAActtct                                             >  2:1815809/1‑150 (MQ=255)
                                                                                                               actAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGcc                                        >  2:899582/1‑150 (MQ=255)
                                                                                                                 tAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCAc                                      <  2:5104529/150‑1 (MQ=255)
                                                                                                                  aaaGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCAcc                                     >  1:1630159/1‑150 (MQ=255)
                                                                                                                     gCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCt                                  <  1:1368358/150‑1 (MQ=255)
                                                                                                                      cAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTg                                 >  1:40431/1‑150 (MQ=255)
                                                                                                                      cAGGTTTCTGGCTGGATGATGCTCAGGTCGTAGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTg                                 >  1:309964/1‑150 (MQ=255)
                                                                                                                         gTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTgggg                              >  2:476074/1‑150 (MQ=255)
                                                                                                                         gTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTgggg                              >  2:456401/1‑150 (MQ=255)
                                                                                                                           ttCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGc                            >  2:5066649/1‑150 (MQ=255)
                                                                                                                            tCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCg                           <  2:2905896/150‑1 (MQ=255)
                                                                                                                              tggctggATGATGCTCAGGTCGTTGATTACCGCGTTCTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTc                         <  2:5633237/150‑1 (MQ=255)
                                                                                                                               ggctggATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCt                        >  2:3212170/1‑150 (MQ=255)
                                                                                                                               ggctggATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCt                        <  2:2025636/150‑1 (MQ=255)
                                                                                                                                 ctggATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTgc                      >  1:821799/1‑150 (MQ=255)
                                                                                                                                     atgatgCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCtt                  >  1:1581742/1‑150 (MQ=255)
                                                                                                                                       gatgCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTtat                >  1:1932227/1‑150 (MQ=255)
                                                                                                                                        atgCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTtata               <  1:1508932/150‑1 (MQ=255)
                                                                                                                                            tCAGGTCGTAGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGt           >  1:5142387/1‑150 (MQ=255)
                                                                                                                                              aGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGttt         <  1:1323847/150‑1 (MQ=255)
                                                                                                                                              aGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGAGCTTATATCGttt         >  2:1513866/1‑150 (MQ=255)
                                                                                                                                               ggTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTc        >  1:323078/1‑150 (MQ=255)
                                                                                                                                               ggTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTc        <  1:2316566/150‑1 (MQ=255)
                                                                                                                                                gTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTcc       >  1:520427/1‑150 (MQ=255)
                                                                                                                                                 tCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTccc      <  2:2813053/150‑1 (MQ=255)
                                                                                                                                                 tCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTccc      >  2:4464841/1‑150 (MQ=255)
                                                                                                                                                  cGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTcccc     <  1:3700368/150‑1 (MQ=255)
                                                                                                                                                   gTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCg    >  1:4801348/1‑150 (MQ=255)
                                                                                                                                                    ttGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCgg   <  1:5279413/150‑1 (MQ=255)
                                                                                                                                                    ttGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCgg   >  1:5308075/1‑150 (MQ=255)
                                                                                                                                                     tGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGt  <  1:2947803/150‑1 (MQ=255)
                                                                                                                                                     tGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGt  <  1:74617/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGT  >  CP009273/568962‑569260

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: