Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 566,230 C→T 5.5% intergenic (+346/‑119) ybcK → / → ybcL DLP12 prophage; putative recombinase/inactive polymorphonuclear leukocyte migration suppressor; DLP12 prophage; secreted protein, UPF0098 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273566,2300CT 5.5% 1091.6 / 34.2 293intergenic (+346/‑119)ybcK/ybcLDLP12 prophage; putative recombinase/inactive polymorphonuclear leukocyte migration suppressor; DLP12 prophage; secreted protein, UPF0098 family
Reads supporting (aligned to +/- strand):  ref base C (149/128);  new base T (12/4);  total (161/132)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 295 total aligned reads displayed.

TGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGG  >  CP009273/566081‑566379
                                                                                                                                                     |                                                                                                                                                     
  aaaaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCaa                                                                                                                                                     <  1:3884642/150‑1 (MQ=255)
   aaaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCaaa                                                                                                                                                    <  1:4957014/150‑1 (MQ=255)
   aaaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCaaa                                                                                                                                                    >  1:617088/1‑150 (MQ=255)
    aaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTAGTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGGTTTGTGTATATTGTGGTaatt                                                                                                                                                   >  1:3532875/1‑148 (MQ=255)
    aaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaat                                                                                                                                                   <  1:1096552/150‑1 (MQ=255)
    aaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaat                                                                                                                                                   <  2:3451832/150‑1 (MQ=255)
       aCGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaataat                                                                                                                                                <  1:1170038/150‑1 (MQ=255)
       aCGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaataat                                                                                                                                                <  1:1251134/150‑1 (MQ=255)
       aCGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaataat                                                                                                                                                >  2:4110863/1‑150 (MQ=255)
        cGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAtt                                                                                                                                               <  2:4818499/150‑1 (MQ=255)
         gTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttg                                                                                                                                              >  2:3995303/1‑150 (MQ=255)
         gTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttg                                                                                                                                              >  1:156120/1‑150 (MQ=255)
         gTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttg                                                                                                                                              <  1:3922975/150‑1 (MQ=255)
          tGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttgg                                                                                                                                             >  2:304253/1‑150 (MQ=255)
            aCAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttgggt                                                                                                                                           >  1:701303/1‑150 (MQ=255)
            aCAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAAttgggt                                                                                                                                           >  1:2692467/1‑150 (MQ=255)
                     tatGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTAc                                                                                                                                  <  2:1970631/150‑1 (MQ=255)
                       tGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACaa                                                                                                                                >  2:760156/1‑150 (MQ=255)
                        gAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAAt                                                                                                                               >  2:96891/1‑150 (MQ=255)
                        gAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAAt                                                                                                                               >  1:3028583/1‑150 (MQ=255)
                         aGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAAtt                                                                                                                              <  1:536121/150‑1 (MQ=255)
                           cAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATtgt                                                                                                                            <  1:331804/150‑1 (MQ=255)
                            aaaaTGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATtgtg                                                                                                                           <  2:4335747/150‑1 (MQ=255)
                             ataTGGTGTATGTCACTTATTTTGAATGGTAGGTTAAATAATAACACCGACTGTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGa                                                                                                                          <  2:1747143/148‑1 (MQ=255)
                               aTGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGAtt                                                                                                                        >  1:64527/1‑150 (MQ=255)
                               aTGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGAtt                                                                                                                        <  1:2177343/150‑1 (MQ=255)
                                tGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTg                                                                                                                       <  2:2190375/150‑1 (MQ=255)
                                tGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTg                                                                                                                       <  2:917480/150‑1 (MQ=255)
                                 ggTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGc                                                                                                                      >  2:705292/1‑150 (MQ=255)
                                   tgtATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCtt                                                                                                                    <  2:1003115/150‑1 (MQ=255)
                                    gtATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCttt                                                                                                                   >  1:809594/1‑150 (MQ=255)
                                            cGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGaa                                                                                                           >  2:3294845/1‑150 (MQ=255)
                                             gTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAAc                                                                                                          >  2:4932868/1‑150 (MQ=255)
                                              tATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACa                                                                                                         >  2:5110216/1‑150 (MQ=255)
                                               aTTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACAt                                                                                                        <  2:2817302/150‑1 (MQ=255)
                                                 tttGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTg                                                                                                      >  1:3685029/1‑150 (MQ=255)
                                                    gAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATAAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTgcgc                                                                                                   >  2:3652301/1‑150 (MQ=255)
                                                    gAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATAAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTGTTGCATGAACAGTgcgc                                                                                                   >  2:3652406/1‑150 (MQ=255)
                                                     aaTGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGcc                                                                                                  <  2:1224190/150‑1 (MQ=255)
                                                       tGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGTAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCtt                                                                                                >  2:1518931/1‑150 (MQ=255)
                                                          tAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTAt                                                                                             >  1:2351232/1‑150 (MQ=255)
                                                          tAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTAt                                                                                             >  1:4895357/1‑150 (MQ=255)
                                                           aGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATg                                                                                            <  2:4658747/150‑1 (MQ=255)
                                                               tAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGTAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATa                                                                                        <  1:2551595/150‑1 (MQ=255)
                                                                  aaaaTAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATg                                                                                     <  1:224394/150‑1 (MQ=255)
                                                                    aaTAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATgag                                                                                   >  1:2157136/1‑150 (MQ=255)
                                                                     aTAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATgaga                                                                                  >  1:2655245/1‑150 (MQ=255)
                                                                      tAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGAt                                                                                 >  1:4484737/1‑150 (MQ=255)
                                                                       aaCACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATa                                                                                <  2:4972720/150‑1 (MQ=255)
                                                                        acacCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATaa                                                                               <  2:4588969/150‑1 (MQ=255)
                                                                          acCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAg                                                                             <  2:2960056/150‑1 (MQ=255)
                                                                           ccGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAgg                                                                            <  1:1546682/150‑1 (MQ=255)
                                                                               cTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAata                                                                        <  2:4809275/150‑1 (MQ=255)
                                                                                tttCGTAGGTGTTACTAATAATAAAGCAGCGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAatat                                                                       <  2:4127585/150‑1 (MQ=255)
                                                                                tttCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGTAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGAAatat                                                                       >  2:3888363/1‑150 (MQ=255)
                                                                                tttCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAatat                                                                       <  2:3818936/150‑1 (MQ=255)
                                                                                tttCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAatat                                                                       >  1:1764589/1‑150 (MQ=255)
                                                                                tttCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAatat                                                                       >  2:1434514/1‑150 (MQ=255)
                                                                                  tCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCa                                                                     >  2:1586746/1‑150 (MQ=255)
                                                                                    gTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaa                                                                   >  1:4381121/1‑150 (MQ=255)
                                                                                    gTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaa                                                                   <  2:155285/150‑1 (MQ=255)
                                                                                    gTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaa                                                                   >  1:390198/1‑150 (MQ=255)
                                                                                      aGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaata                                                                 >  1:4250782/1‑150 (MQ=255)
                                                                                        gtgtTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaataaa                                                               >  2:373437/1‑150 (MQ=255)
                                                                                         tgtTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaataaaa                                                              >  1:327726/1‑150 (MQ=255)
                                                                                           ttACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCaaataaaata                                                            <  1:1405251/150‑1 (MQ=255)
                                                                                              cTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGTAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGAAATATCAAATAAAATAACg                                                         >  1:2111279/1‑150 (MQ=255)
                                                                                                  taataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGTATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAg                                                     >  2:4538446/1‑150 (MQ=255)
                                                                                                  taataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAg                                                     <  1:4471162/150‑1 (MQ=255)
                                                                                                   aataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAgg                                                    >  1:2818223/1‑150 (MQ=255)
                                                                                                    ataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGt                                                   >  2:4262825/1‑150 (MQ=255)
                                                                                                    ataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGt                                                   >  2:815339/1‑150 (MQ=255)
                                                                                                     taaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTc                                                  <  1:3292655/150‑1 (MQ=255)
                                                                                                      aaaGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCa                                                 >  2:1476714/1‑150 (MQ=255)
                                                                                                          cAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAAc                                             >  1:4429149/1‑150 (MQ=255)
                                                                                                           agagTTTTTAGATAGTAACAATGTTCTTTGTGTATATTGTGGCAAAAAATTGGGTTGGGGGTACAATTGTGATAGCTTTTGCATGAACATAGCGCCTTTATTCATAATGAGATAAAGGAATATCAAATAAAATAAAGAAACGACATAACa                                            >  2:4406715/1‑150 (MQ=255)
                                                                                                                 tttAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTAGGGGGAACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAgg                                      <  1:5138562/150‑1 (MQ=255)
                                                                                                                     gATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGtttt                                  <  1:3162005/150‑1 (MQ=255)
                                                                                                                     gATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGtttt                                  >  2:3176479/1‑150 (MQ=255)
                                                                                                                      aTAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGttttt                                 <  1:4449650/150‑1 (MQ=255)
                                                                                                                       tAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTa                                <  1:1306745/150‑1 (MQ=255)
                                                                                                                       tAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTa                                >  1:3432929/1‑150 (MQ=255)
                                                                                                                           aTCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGaa                            >  1:4307608/1‑150 (MQ=255)
                                                                                                                             cAATGTGCTTTGTGTATATTGTGGTAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGAAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGaaaa                          >  1:2049641/1‑150 (MQ=255)
                                                                                                                             cAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGaaaa                          >  2:2104734/1‑150 (MQ=255)
                                                                                                                                tgtgCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAacac                       >  1:241187/1‑150 (MQ=255)
                                                                                                                                 gtgCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACt                      >  1:1598958/1‑150 (MQ=255)
                                                                                                                                    cTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTAt                   >  2:2685363/1‑150 (MQ=255)
                                                                                                                                     tttGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATc                  >  1:3716361/1‑150 (MQ=255)
                                                                                                                                     tttGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATc                  <  1:383485/150‑1 (MQ=255)
                                                                                                                                      ttGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCg                 >  1:1464280/1‑150 (MQ=255)
                                                                                                                                       tgtgTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGt                >  2:4044090/1‑150 (MQ=255)
                                                                                                                                        gtgtATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGtt               >  2:3535131/1‑150 (MQ=255)
                                                                                                                                         tgtATATTGTGGCAAATAATTGTGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGttt              <  2:1061071/150‑1 (MQ=255)
                                                                                                                                          gtATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTc             >  2:2326479/1‑150 (MQ=255)
                                                                                                                                           tataTTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCa            <  2:4882523/150‑1 (MQ=255)
                                                                                                                                            atatTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCaa           >  2:620105/1‑150 (MQ=255)
                                                                                                                                            atatTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCaa           >  2:1848479/1‑150 (MQ=255)
                                                                                                                                                tgtgGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGt       >  1:1343356/1‑150 (MQ=255)
                                                                                                                                                tgtgGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGt       <  1:1881573/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGG  >  CP009273/566081‑566379

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: