Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,629,460 T→G 10.2% A215A (GCA→GCC stfQ ← Qin prophage; putative side tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,629,4600TG10.2% 760.8 / 82.7 285A215A (GCA→GCCstfQQin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base T (136/120);  new base G (16/13);  total (152/133)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 286 total aligned reads displayed.

GACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGG  >  CP009273/1629311‑1629608
                                                                                                                                                     |                                                                                                                                                    
    gCGAGTGGGTATGATCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTAGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGAAc                                                                                                                                                  >  2:3016986/1‑150 (MQ=255)
         tGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCAGTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGaa                                                                                                                                             >  1:2870167/1‑150 (MQ=255)
          gggTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAAc                                                                                                                                            <  2:2686722/150‑1 (MQ=255)
          gggTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAAc                                                                                                                                            <  2:2686436/150‑1 (MQ=255)
            gTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAActc                                                                                                                                          >  2:2264374/1‑150 (MQ=255)
             tATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAActct                                                                                                                                         <  1:907690/150‑1 (MQ=255)
                 aGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTtgtg                                                                                                                                     >  1:3640994/1‑150 (MQ=255)
                   cACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGtt                                                                                                                                   <  1:1603396/150‑1 (MQ=255)
                     ccATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTgg                                                                                                                                 <  1:3923476/150‑1 (MQ=255)
                         cTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTgcg                                                                                                                             <  2:4433111/150‑1 (MQ=255)
                         cTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTgcg                                                                                                                             >  1:3217776/1‑150 (MQ=255)
                           gatgatGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCa                                                                                                                           <  1:2251352/150‑1 (MQ=255)
                            atgatgTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCAc                                                                                                                          >  1:1874281/1‑150 (MQ=255)
                              gatgTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCg                                                                                                                        >  1:2517651/1‑150 (MQ=255)
                                     tCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTAt                                                                                                                 >  2:865374/1‑150 (MQ=255)
                                     tCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTAt                                                                                                                 <  1:2213217/150‑1 (MQ=255)
                                      cccTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTAtt                                                                                                                <  2:1620325/150‑1 (MQ=255)
                                       ccTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCTACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTAttt                                                                                                               <  1:1213460/150‑1 (MQ=255)
                                           cATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCgg                                                                                                           <  1:2878916/150‑1 (MQ=255)
                                           cATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCgg                                                                                                           >  2:4304290/1‑150 (MQ=255)
                                            aTTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGt                                                                                                          >  2:971231/1‑150 (MQ=255)
                                             ttACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTc                                                                                                         >  1:927747/1‑150 (MQ=255)
                                               aCGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTccc                                                                                                       <  1:2859016/150‑1 (MQ=255)
                                               aCGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTccc                                                                                                       >  1:3058723/1‑150 (MQ=255)
                                               aCGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTccc                                                                                                       >  1:2086589/1‑150 (MQ=255)
                                                   gTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTa                                                                                                   <  1:306156/150‑1 (MQ=255)
                                                     cTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTAAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAAt                                                                                                 <  2:258280/150‑1 (MQ=255)
                                                      tGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATa                                                                                                <  2:2505078/150‑1 (MQ=255)
                                                      tGGCCACTACCGCTTGTTGTGCTCATAATCCCCACGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATa                                                                                                <  1:4481185/150‑1 (MQ=255)
                                                       ggCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATAc                                                                                               >  2:2130548/1‑150 (MQ=255)
                                                          cACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTAt                                                                                            <  2:4378604/150‑1 (MQ=255)
                                                           aCTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATg                                                                                           >  1:15128/1‑150 (MQ=255)
                                                               ccGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTa                                                                                       <  1:1171877/150‑1 (MQ=255)
                                                               ccGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTa                                                                                       >  2:1067528/1‑150 (MQ=255)
                                                                cGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTAt                                                                                      <  1:2510724/150‑1 (MQ=255)
                                                                 gCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATg                                                                                     <  1:3397419/150‑1 (MQ=255)
                                                                        gtgCTCATAAGCCCCGCGCTGAGATTTGAAATCGCGGTCTAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCg                                                                              <  2:3751528/150‑1 (MQ=255)
                                                                           cTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGgtg                                                                           >  2:3319797/1‑150 (MQ=255)
                                                                             cATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGtt                                                                         <  1:1741794/150‑1 (MQ=255)
                                                                             cATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGtt                                                                         <  1:843647/150‑1 (MQ=255)
                                                                             cATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGtt                                                                         >  1:4217051/1‑150 (MQ=255)
                                                                             cATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGtt                                                                         >  1:1973170/1‑150 (MQ=255)
                                                                              aTAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTa                                                                        >  2:4028263/1‑150 (MQ=255)
                                                                               tAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTAt                                                                       <  1:4047693/150‑1 (MQ=255)
                                                                                aaTCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTAtt                                                                      <  2:2793956/150‑1 (MQ=255)
                                                                                aaTCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTAtt                                                                      >  1:2162589/1‑150 (MQ=255)
                                                                                  tCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCg                                                                    <  1:3906680/150‑1 (MQ=255)
                                                                                    cccGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTg                                                                  >  1:1500107/1‑150 (MQ=255)
                                                                                      cgcgCTTAGATTTGAAATCGCGGTATAACCATAAGGGAAAATGCTCGTGTACGTGCCACCAAATGCACCGGAACTCGTGGGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGAGTTATTCGTGGa                                                                >  2:3155636/1‑150 (MQ=255)
                                                                                         gcTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGAttt                                                             >  1:3240076/1‑150 (MQ=255)
                                                                                          cTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTg                                                            >  1:682255/1‑150 (MQ=255)
                                                                                          cTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTGATTCGTGGATTTg                                                            >  1:2372863/1‑150 (MQ=255)
                                                                                              gATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTc                                                        >  1:376496/1‑150 (MQ=255)
                                                                                              gATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTc                                                        >  1:403074/1‑150 (MQ=255)
                                                                                               aTTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTcc                                                       <  1:2692355/150‑1 (MQ=255)
                                                                                                     aaTCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATc                                                 >  2:2980118/1‑150 (MQ=255)
                                                                                                      aTCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATa                                                                                                >  2:1591557/1‑102 (MQ=255)
                                                                                                      aTCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATa                                                                                                <  1:1590794/102‑1 (MQ=255)
                                                                                                      aTCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGGGTTATTCGTGGATTTGTTTCCGTAATCa                                                >  1:3321181/1‑150 (MQ=255)
                                                                                                        cgcgGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCaaa                                              <  2:2768723/150‑1 (MQ=255)
                                                                                                        cgcgGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCaaa                                              <  1:458739/150‑1 (MQ=255)
                                                                                                        cgcgGTATAACCATTAGGGAAAAGGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCaaa                                              <  2:3250570/150‑1 (MQ=255)
                                                                                                         gcgGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAac                                             <  2:2513431/150‑1 (MQ=255)
                                                                                                           ggTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacga                                           >  2:1493972/1‑150 (MQ=255)
                                                                                                            gTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgac                                          <  2:2625243/150‑1 (MQ=255)
                                                                                                             tataACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                         >  2:700440/1‑150 (MQ=255)
                                                                                                             tataACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                         >  2:2264838/1‑150 (MQ=255)
                                                                                                              ataACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAAcgacga                                        >  2:4485170/1‑150 (MQ=255)
                                                                                                                acccATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtg                                      >  2:134803/3‑150 (MQ=255)
                                                                                                                aaCCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGTTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtg                                      <  1:2926835/150‑1 (MQ=255)
                                                                                                                aaCCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtg                                      >  2:3375093/1‑150 (MQ=255)
                                                                                                                 aCCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtgt                                     >  1:3454545/1‑150 (MQ=255)
                                                                                                                  ccATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtgtg                                    >  1:877993/1‑150 (MQ=255)
                                                                                                                  ccATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAtgtg                                    >  1:3547595/1‑150 (MQ=255)
                                                                                                                    aTTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGt                                  >  1:567347/1‑150 (MQ=255)
                                                                                                                     ttAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGtt                                 <  2:3486779/150‑1 (MQ=255)
                                                                                                                     ttAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGtt                                 <  2:643808/150‑1 (MQ=255)
                                                                                                                      tAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGttt                                >  2:3794641/1‑150 (MQ=255)
                                                                                                                      tAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGttt                                <  2:3211113/150‑1 (MQ=255)
                                                                                                                      tAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGttt                                >  2:1358141/1‑150 (MQ=255)
                                                                                                                           aaaaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTc                           >  2:3392166/1‑150 (MQ=255)
                                                                                                                           aaaaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTc                           <  2:2058743/150‑1 (MQ=255)
                                                                                                                           aaaaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGAGGTGGTTTCCGTc                           >  2:2998860/1‑150 (MQ=255)
                                                                                                                            aaaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTcc                          <  2:1102592/150‑1 (MQ=255)
                                                                                                                            aaaTGCACGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTcc                          <  2:1807986/150‑1 (MQ=255)
                                                                                                                             aaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTccc                         >  2:580051/1‑150 (MQ=255)
                                                                                                                              aTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTcccc                        <  2:34293/150‑1 (MQ=255)
                                                                                                                                gCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCaa                      <  1:4022817/150‑1 (MQ=255)
                                                                                                                                 cTCGTGTTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCaaa                     >  1:1283802/1‑150 (MQ=255)
                                                                                                                                    gtgtTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATcc                  >  1:1933112/1‑150 (MQ=255)
                                                                                                                                     tgtTCGTGCCACCAAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCg                 <  2:2517651/150‑1 (MQ=255)
                                                                                                                                       ttCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTa               >  1:974491/1‑150 (MQ=255)
                                                                                                                                         cGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACt             <  1:2980259/150‑1 (MQ=255)
                                                                                                                                         cGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACt             <  1:4304290/150‑1 (MQ=255)
                                                                                                                                          gTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTg            <  2:1562055/150‑1 (MQ=255)
                                                                                                                                             ccaccaAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGAt         >  2:1570082/1‑150 (MQ=255)
                                                                                                                                              caccaAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATg        <  2:3843542/150‑1 (MQ=255)
                                                                                                                                              caccaAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATg        >  2:1553856/1‑150 (MQ=255)
                                                                                                                                                 caAAGGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgc     >  2:2754269/1‑150 (MQ=255)
                                                                                                                                                   aaTGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgctg   <  1:1121749/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGG  >  CP009273/1629311‑1629608

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: