Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,425,324 G→A 6.8% S353S (AGC→AGT cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,425,3240GA 6.8% 920.7 / 37.9 266S353S (AGC→AGTcpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (118/130);  new base A (11/7);  total (129/137)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

CACCGGCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCAT  >  CP009273/4425177‑4425472
                                                                                                                                                   |                                                                                                                                                    
cACCGGCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCt                                                                                                                                                    <  2:1121097/150‑1 (MQ=255)
  ccGGCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGa                                                                                                                                                  >  2:2828782/1‑150 (MQ=255)
   cGGCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGAc                                                                                                                                                 >  1:2868432/1‑150 (MQ=255)
    ggCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg                                                                                                                                                >  2:4342282/1‑150 (MQ=255)
      cAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGaa                                                                                                                                              >  2:3579448/1‑150 (MQ=255)
       aGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAAc                                                                                                                                             >  1:226040/1‑150 (MQ=255)
        gTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACt                                                                                                                                            <  2:441623/150‑1 (MQ=255)
           ttGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCGGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTACTGACGAACTGGc                                                                                                                                         >  1:2338851/1‑150 (MQ=255)
           ttGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCGGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTACTGACGAACTGGc                                                                                                                                         >  1:2338808/1‑150 (MQ=255)
           ttGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCGGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTACTGACGAACTGGc                                                                                                                                         >  1:2098626/1‑150 (MQ=255)
             gCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCgt                                                                                                                                       >  1:1706798/1‑150 (MQ=255)
              ccAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCgtg                                                                                                                                      <  2:183631/150‑1 (MQ=255)
               cAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGAATTACCGATTGGCTTGCTGACGAACTGGCgtgt                                                                                                                                     <  1:671961/150‑1 (MQ=255)
                aGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGtgtg                                                                                                                                    <  2:4194519/150‑1 (MQ=255)
                  aTCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGc                                                                                                                                  >  2:427741/1‑150 (MQ=255)
                      ggATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCg                                                                                                                              <  1:3360374/150‑1 (MQ=255)
                        aTCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTg                                                                                                                            <  2:1780919/150‑1 (MQ=255)
                          cGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGAt                                                                                                                          >  2:1709062/1‑150 (MQ=255)
                            cccTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCg                                                                                                                        >  1:2884648/1‑150 (MQ=255)
                            cccTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCg                                                                                                                        <  2:3709673/150‑1 (MQ=255)
                               tGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCt                                                                                                                     >  1:1168522/1‑150 (MQ=255)
                                gAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAGTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCtt                                                                                                                    >  2:253564/1‑150 (MQ=255)
                                 aaTGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCttt                                                                                                                   >  1:1597150/1‑150 (MQ=255)
                                 aaTGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCttt                                                                                                                   >  2:3049709/1‑150 (MQ=255)
                                      aaTGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAgtg                                                                                                              <  1:4324525/150‑1 (MQ=255)
                                      aaTGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAgtg                                                                                                              >  1:3622568/1‑150 (MQ=255)
                                       aTGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAgtgt                                                                                                             >  1:2755181/1‑150 (MQ=255)
                                        tGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGtt                                                                                                            <  1:938579/150‑1 (MQ=255)
                                            cGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTa                                                                                                        <  2:4061416/150‑1 (MQ=255)
                                             gACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTAc                                                                                                       >  1:1364507/1‑150 (MQ=255)
                                               cATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCa                                                                                                     <  1:3919621/150‑1 (MQ=255)
                                               cATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCa                                                                                                     >  1:1057711/1‑150 (MQ=255)
                                               cATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCa                                                                                                     <  1:3565608/150‑1 (MQ=255)
                                               cATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCa                                                                                                     <  1:3423167/150‑1 (MQ=255)
                                                  aCGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCt                                                                                                  >  2:3937324/1‑150 (MQ=255)
                                                  aCGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCt                                                                                                  >  1:67708/1‑150 (MQ=255)
                                                   cGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCtt                                                                                                 <  2:3133649/150‑1 (MQ=255)
                                                   cGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCtt                                                                                                 >  1:1435315/1‑150 (MQ=255)
                                                    gCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTg                                                                                                >  2:1668806/1‑150 (MQ=255)
                                                     cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGc                                                                                               <  2:2644605/150‑1 (MQ=255)
                                                     cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGc                                                                                               >  2:3712062/1‑150 (MQ=255)
                                                       ttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTg                                                                                             <  2:2401922/150‑1 (MQ=255)
                                                         ttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTc                                                                                           <  2:67647/150‑1 (MQ=255)
                                                            cgcgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTc                                                                                        <  1:1567538/150‑1 (MQ=255)
                                                              cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTccg                                                                                      <  1:2291841/150‑1 (MQ=255)
                                                                ttgttCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTGTGAGTGTTTCTACCAGCTTGCTGTCTTccgcc                                                                                    >  2:4444829/1‑150 (MQ=255)
                                                                  gttCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCcgccgc                                                                                  <  2:1981901/150‑1 (MQ=255)
                                                                  gttCACCACCTGCACGGTCGGATCGTCCGGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTACTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCcgccgc                                                                                  <  1:2279500/150‑1 (MQ=255)
                                                                     caccacCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAg                                                                               >  1:3189027/1‑150 (MQ=255)
                                                                        caccTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGa                                                                            >  2:2744794/1‑150 (MQ=255)
                                                                             gCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGAttttt                                                                       >  1:647527/1‑150 (MQ=255)
                                                                              cACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGAtttttt                                                                      <  2:970931/150‑1 (MQ=255)
                                                                               aCGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTtta                                                                     >  1:2447451/1‑150 (MQ=255)
                                                                               aCGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTtta                                                                     >  1:277784/1‑150 (MQ=255)
                                                                                 ggTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTttatt                                                                   >  2:2310115/1‑150 (MQ=255)
                                                                                  gTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTttatta                                                                  <  1:3594584/150‑1 (MQ=255)
                                                                                  gTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTttatta                                                                  <  1:1325010/150‑1 (MQ=255)
                                                                                   tCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAg                                                                 >  2:3395570/1‑150 (MQ=255)
                                                                                   tCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAg                                                                 >  2:1060910/1‑150 (MQ=255)
                                                                                      gATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGa                                                              >  2:2481140/1‑150 (MQ=255)
                                                                                      gATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGa                                                              >  1:4244423/1‑150 (MQ=255)
                                                                                      gATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGa                                                              >  2:792449/1‑150 (MQ=255)
                                                                                      gATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGa                                                              <  2:3837694/150‑1 (MQ=255)
                                                                                         cGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGt                                                           <  2:1466225/150‑1 (MQ=255)
                                                                                         cGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGt                                                           <  1:3668369/150‑1 (MQ=255)
                                                                                          gTCCTGCACCAGCGCCAGATTGCTTTACATATTGTCGGCAGATTTACCGATTGGCTTGCTGTCGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTc                                                          <  1:1077205/150‑1 (MQ=255)
                                                                                           tCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCg                                                         <  2:2334104/150‑1 (MQ=255)
                                                                                            ccTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGt                                                        <  2:1005308/150‑1 (MQ=255)
                                                                                             cTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTa                                                       <  1:1176600/150‑1 (MQ=255)
                                                                                              tGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTaa                                                      >  1:2309431/1‑150 (MQ=255)
                                                                                                   cAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCgg                                                 <  2:226009/150‑1 (MQ=255)
                                                                                                            aTAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTc                                        <  2:1534680/150‑1 (MQ=255)
                                                                                                             tAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTcc                                       >  1:1341367/1‑150 (MQ=255)
                                                                                                             tAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTACTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTcc                                       <  1:1974910/150‑1 (MQ=255)
                                                                                                                  ataCATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCttt                                  <  1:1682973/150‑1 (MQ=255)
                                                                                                                  ataCATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCttt                                  <  1:1683507/150‑1 (MQ=255)
                                                                                                                   taCATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTc                                 <  1:3962426/150‑1 (MQ=255)
                                                                                                                     cATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGc                               >  2:344239/1‑150 (MQ=255)
                                                                                                                       tatTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCt                             <  1:2509842/150‑1 (MQ=255)
                                                                                                                       tatTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCt                             >  2:1292246/1‑150 (MQ=255)
                                                                                                                       tatTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCt                             >  2:1292149/1‑150 (MQ=255)
                                                                                                                        atTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTg                            >  2:196629/1‑150 (MQ=255)
                                                                                                                        atTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTg                            >  2:4114921/1‑150 (MQ=255)
                                                                                                                          tGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCg                          <  2:2787508/150‑1 (MQ=255)
                                                                                                                           gTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGt                         >  2:2832477/1‑150 (MQ=255)
                                                                                                                           gTCGGCAGATTTACCGATTGGCGTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGt                         <  2:516370/150‑1 (MQ=255)
                                                                                                                            tCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTc                        >  1:818384/1‑150 (MQ=255)
                                                                                                                               gCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCAcc                     <  1:2675509/150‑1 (MQ=255)
                                                                                                                                cAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCt                    <  2:1737510/150‑1 (MQ=255)
                                                                                                                                    tttACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCa                >  2:2276571/1‑150 (MQ=255)
                                                                                                                                         cGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTAc           <  2:605989/150‑1 (MQ=255)
                                                                                                                                          gATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTAcc          <  1:2210351/150‑1 (MQ=255)
                                                                                                                                          gATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTAcc          <  1:3634594/150‑1 (MQ=255)
                                                                                                                                          gATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTAcc          >  2:2682955/1‑150 (MQ=255)
                                                                                                                                          gATTGGCTTGCAGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTAATAGCGATTTAGTAAATCGGTCGTGCTTGCGCTTTCGCATGCGTCTCCTGCCATGTAca          >  2:2912533/1‑149 (MQ=255)
                                                                                                                                           aTTGGCTTACTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCg         <  2:812131/150‑1 (MQ=255)
                                                                                                                                              ggCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTg      >  2:4509045/1‑150 (MQ=255)
                                                                                                                                              ggCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTg      >  2:2194414/1‑150 (MQ=255)
                                                                                                                                               gCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGt     <  1:654876/150‑1 (MQ=255)
                                                                                                                                                cTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATGTACCGCTGTc    >  2:3637273/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACCGGCAGTTTTGCCAGATCCGGATCGCCCTGAATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCAT  >  CP009273/4425177‑4425472

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: