Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,556,744 T→C 7.9% intergenic (+268/‑207) yffP → / → yffQ CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,556,7440TC 7.9% 1066.6 / 58.8 343intergenic (+268/‑207)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (157/158);  new base C (13/14);  total (171/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

GCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGT  >  CP009273/2556595‑2556893
                                                                                                                                                     |                                                                                                                                                     
gcgcCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAAt                                                                                                                                                       <  2:2768950/150‑1 (MQ=255)
  gcCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCa                                                                                                                                                     <  1:3269762/150‑1 (MQ=255)
    cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAAc                                                                                                                                                   >  1:2091892/1‑150 (MQ=255)
    cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCaaa                                                                                                                                                   >  2:4450286/1‑149 (MQ=255)
    cAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAAc                                                                                                                                                   >  2:2553635/1‑150 (MQ=255)
      aaaaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTc                                                                                                                                                 >  2:202477/1‑150 (MQ=255)
       aaaTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTcc                                                                                                                                                >  2:4934045/1‑150 (MQ=255)
          tAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCccac                                                                                                                                             >  1:370004/1‑150 (MQ=255)
             aGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaa                                                                                                                                          >  2:3967770/1‑150 (MQ=255)
             aGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaa                                                                                                                                          >  2:3076908/1‑150 (MQ=255)
              ggCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCaaa                                                                                                                                         >  2:5099227/1‑150 (MQ=255)
                cccAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGc                                                                                                                                       >  2:3732582/1‑150 (MQ=255)
                 ccAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGcc                                                                                                                                      <  2:2127188/150‑1 (MQ=255)
                  cAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCca                                                                                                                                     <  2:3926688/150‑1 (MQ=255)
                   aGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCcaa                                                                                                                                    >  2:1922622/1‑150 (MQ=255)
                       aGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAg                                                                                                                                >  1:4589387/1‑150 (MQ=255)
                       aGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTGTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAg                                                                                                                                <  1:976179/150‑1 (MQ=255)
                         ttAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGAt                                                                                                                              <  2:1107856/150‑1 (MQ=255)
                          tAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATg                                                                                                                             >  1:1713870/1‑150 (MQ=255)
                          tAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATg                                                                                                                             <  1:3592224/150‑1 (MQ=255)
                           aGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATgg                                                                                                                            >  1:3967678/1‑150 (MQ=255)
                                 cTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTg                                                                                                                      <  1:4010957/150‑1 (MQ=255)
                                 cTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTg                                                                                                                      >  1:4336491/1‑150 (MQ=255)
                                  tACCCGGGCCTTTCTCACACTCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGa                                                                                                                     <  1:2372620/150‑1 (MQ=255)
                                  tACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGa                                                                                                                     <  1:3920600/150‑1 (MQ=255)
                                    cccGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAAt                                                                                                                   <  1:626432/150‑1 (MQ=255)
                                         gCCTTTCTCACCATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTtgt                                                                                                              <  2:4000282/150‑1 (MQ=255)
                                          ccTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTtgtt                                                                                                             >  2:5192182/1‑149 (MQ=255)
                                            tttCTCACACTCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTAt                                                                                                           <  1:1214032/150‑1 (MQ=255)
                                            tttCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTAt                                                                                                           <  1:4213456/150‑1 (MQ=255)
                                             ttCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATg                                                                                                          <  1:1019966/150‑1 (MQ=255)
                                             ttCTCACAATCCCCGCTGACATTATATCACCCATCTTTACGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTACTTTCTGAACATTCAAATCCACTCCCACCAAAGCCAACACGAGGGGAGGGGACTGTTGTTATg                                                                                                          >  2:3109249/1‑150 (MQ=255)
                                              tctcACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGc                                                                                                         >  1:4679318/1‑150 (MQ=255)
                                                tcACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAg                                                                                                       <  2:244409/150‑1 (MQ=255)
                                                tcACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTGGTTATGCAg                                                                                                       >  2:1333309/1‑150 (MQ=255)
                                                 cacaATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGa                                                                                                      <  2:1727049/150‑1 (MQ=255)
                                                  acaATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGAt                                                                                                     >  2:126521/1‑150 (MQ=255)
                                                   caATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATc                                                                                                    <  1:2282365/150‑1 (MQ=255)
                                                    aaTCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAGGCCAACAAGATGGGAGGTGAATTTGGTTATGCAGATCg                                                                                                   >  1:3639130/1‑150 (MQ=255)
                                                        cccGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCttt                                                                                               >  1:1772862/1‑150 (MQ=255)
                                                         ccGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCtttt                                                                                              <  2:4275916/150‑1 (MQ=255)
                                                         ccGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCtttt                                                                                              <  2:1849040/150‑1 (MQ=255)
                                                            cTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAAGTGAATTTTGTTATGCAGATCGCTTTTTcc                                                                                           >  1:405920/1‑150 (MQ=255)
                                                              gACATTATATCACCCATCTGTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCt                                                                                         <  2:4304847/150‑1 (MQ=255)
                                                               aCATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCtt                                                                                        <  1:2647052/150‑1 (MQ=255)
                                                               aCATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCtt                                                                                        <  2:4061912/150‑1 (MQ=255)
                                                               aCATTAGATCACCCATCTTTCCGTGTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAGTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCACCAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCtt                                                                                        <  2:234240/150‑1 (MQ=255)
                                                                cATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTc                                                                                       >  1:976921/1‑150 (MQ=255)
                                                                cATTATATCACCCATCTGTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTc                                                                                       <  2:2606332/150‑1 (MQ=255)
                                                                 aTTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCa                                                                                      >  1:1545385/1‑150 (MQ=255)
                                                                  ttATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAg                                                                                     >  1:1480031/1‑150 (MQ=255)
                                                                   tataTCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAgg                                                                                    <  1:4423932/150‑1 (MQ=255)
                                                                        cACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCAttt                                                                               >  1:578213/1‑150 (MQ=255)
                                                                          cccATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCAttttt                                                                             <  1:3075066/150‑1 (MQ=255)
                                                                           ccATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTg                                                                            >  1:942753/1‑150 (MQ=255)
                                                                           ccATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTg                                                                            <  2:4204723/150‑1 (MQ=255)
                                                                             aTCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGtt                                                                          <  1:309094/150‑1 (MQ=255)
                                                                                 ttCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAAt                                                                      >  2:3020653/1‑150 (MQ=255)
                                                                                  tCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAAtt                                                                     <  2:354449/150‑1 (MQ=255)
                                                                                  tCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAAtt                                                                     >  1:5009051/1‑150 (MQ=255)
                                                                                  tCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAAtt                                                                     >  1:5215814/1‑150 (MQ=255)
                                                                                   ccGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAATGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATtc                                                                    >  1:3864531/1‑150 (MQ=255)
                                                                                      ttttCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATtcaat                                                                 <  2:3225286/150‑1 (MQ=255)
                                                                                         tCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATtcaatcaa                                                              <  2:5287080/150‑1 (MQ=255)
                                                                                              tGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGaaa                                                         <  1:3548890/150‑1 (MQ=255)
                                                                                               gggAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAg                                                        >  2:1656624/1‑150 (MQ=255)
                                                                                               gggAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAg                                                        >  2:3804122/1‑150 (MQ=255)
                                                                                                ggAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAgg                                                       >  2:2699108/1‑150 (MQ=255)
                                                                                                 gagaCTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGc                                                      >  2:3431457/1‑150 (MQ=255)
                                                                                                  agaCTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCa                                                     <  1:2736840/150‑1 (MQ=255)
                                                                                                        aTGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTccc                                               <  1:195751/150‑1 (MQ=255)
                                                                                                         tGGCCGCCGACGTTGCACAATTGATTAATTCCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCca                                              <  1:4782328/150‑1 (MQ=255)
                                                                                                                cgACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCcaacaacaa                                       >  2:4578836/1‑150 (MQ=255)
                                                                                                                cgACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCATTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGAAGATCCcaaaaaaaa                                       >  2:4537514/1‑150 (MQ=255)
                                                                                                                     ttGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATaa                                  >  2:1423529/1‑150 (MQ=255)
                                                                                                                      tGCACAATTGATTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATaaa                                 <  2:4336491/150‑1 (MQ=255)
                                                                                                                         acaATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGAt                              <  2:480628/150‑1 (MQ=255)
                                                                                                                           aaTTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATgt                            >  2:786978/1‑149 (MQ=255)
                                                                                                                            aTTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATggg                           >  1:107570/1‑150 (MQ=255)
                                                                                                                             ttGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAGTTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGa                          <  1:1341820/150‑1 (MQ=255)
                                                                                                                             ttGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGa                          >  2:919104/1‑150 (MQ=255)
                                                                                                                             ttGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGa                          <  1:832424/150‑1 (MQ=255)
                                                                                                                                aTTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACACCAAAATAAAGATGGGAGGt                       >  1:3218322/1‑150 (MQ=255)
                                                                                                                                aTTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                       >  2:2167744/1‑150 (MQ=255)
                                                                                                                                aTTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                       >  1:2963528/1‑150 (MQ=255)
                                                                                                                                aTTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAAGAAAGATGGGAGGt                       >  1:1838769/1‑150 (MQ=255)
                                                                                                                                aTTAATTTCTGAACATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                       >  2:3582139/1‑150 (MQ=255)
                                                                                                                                   aaTTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                    >  1:4365866/1‑150 (MQ=255)
                                                                                                                                     tttCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                  <  1:4519970/150‑1 (MQ=255)
                                                                                                                                         tGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaa              >  1:4936175/1‑150 (MQ=255)
                                                                                                                                          gAAAATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaaa             <  2:5040421/150‑1 (MQ=255)
                                                                                                                                           aaCATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc            <  1:5242043/150‑1 (MQ=255)
                                                                                                                                            aCATTCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACg           >  2:4378569/1‑150 (MQ=255)
                                                                                                                                             cATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACgg          <  1:4919015/150‑1 (MQ=255)
                                                                                                                                             cATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACgg          >  2:3645855/1‑150 (MQ=255)
                                                                                                                                               atCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTg        <  2:2918390/149‑1 (MQ=255)
                                                                                                                                                tCAAACCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTgg       <  2:1232412/150‑1 (MQ=255)
                                                                                                                                                  aaaTCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGTAACGGTGGGa     >  2:3247609/1‑150 (MQ=255)
                                                                                                                                                  aaaTCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa     >  2:574544/1‑150 (MQ=255)
                                                                                                                                                   aaTCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTTATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGaa    >  2:197861/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGT  >  CP009273/2556595‑2556893

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: