Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,874,100 C→A 9.9% R161L (CGC→CTC)  recF ← gap repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,874,1000CA 9.9% 1275.4 / 75.4 350R161L (CGC→CTC) recFgap repair protein
Reads supporting (aligned to +/- strand):  ref base C (169/146);  new base A (23/12);  total (192/158)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 362 total aligned reads displayed.

GAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATC  >  CP009273/3873951‑3874249
                                                                                                                                                     |                                                                                                                                                     
gaAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgc                                                                                                                                                       >  1:214569/1‑150 (MQ=255)
gaAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgc                                                                                                                                                       >  2:2737955/1‑150 (MQ=255)
 aaaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgcg                                                                                                                                                      >  2:3196700/1‑150 (MQ=255)
 aaaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgcg                                                                                                                                                      >  1:3407492/1‑150 (MQ=255)
 aaaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCAAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgcg                                                                                                                                                      <  1:1480514/150‑1 (MQ=255)
    ttGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCgctg                                                                                                                                                   >  1:1202533/1‑150 (MQ=255)
    ttGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTAGTAACGTGTCACCTGGCGCAGCGCCGCATTGCgctg                                                                                                                                                   >  2:448803/1‑150 (MQ=255)
        ttACCGGGATCAGCCATCTCGGCCGCGAGACCGGCGCTAGACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTGATCCCACGGGCGTAGCTGTTCGTAACGTGTCCCCTGGCGCAGCGCCGCATTGCGCTGCTTg                                                                                                                                               <  1:812733/150‑1 (MQ=255)
         tACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGa                                                                                                                                              >  2:1577570/1‑150 (MQ=255)
            aGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCa                                                                                                                                           >  1:468061/1‑150 (MQ=255)
             ggTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCaa                                                                                                                                          >  2:2129231/1‑150 (MQ=255)
             ggTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCaa                                                                                                                                          >  2:1153654/1‑150 (MQ=255)
               tATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATc                                                                                                                                        <  1:4843144/150‑1 (MQ=255)
                aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCg                                                                                                                                       <  1:2762836/150‑1 (MQ=255)
                  cAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCt                                                                                                                                     >  2:301918/1‑150 (MQ=255)
                   aGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGcgg                                                                                                                                    >  1:3229947/1‑148 (MQ=255)
                   aGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATGTGCTCCGCCAGCGGGATCAGCTCTTTCTCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCtt                                                                                                                                    <  2:4281725/150‑1 (MQ=255)
                      cATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTGgag                                                                                                                                 >  1:1301283/1‑150 (MQ=255)
                       atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTgaga                                                                                                                                <  2:4413092/150‑1 (MQ=255)
                         atCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGAGCAGCGCCGCATTGCGCTGCTTGAGCAAACGCTTGAGAtt                                                                                                                              <  1:2039894/150‑1 (MQ=255)
                          tCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTg                                                                                                                             >  1:776849/1‑150 (MQ=255)
                          tCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTg                                                                                                                             <  2:1089921/150‑1 (MQ=255)
                           cGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCACTCGCTGGAGAGTGc                                                                                                                            >  1:2948215/1‑150 (MQ=255)
                           cGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGc                                                                                                                            >  2:925136/1‑150 (MQ=255)
                            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCt                                                                                                                           <  2:4211422/150‑1 (MQ=255)
                            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCt                                                                                                                           >  2:4746714/1‑150 (MQ=255)
                            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTGGAGATTGCt                                                                                                                           >  1:4538685/1‑150 (MQ=255)
                               cgcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCa                                                                                                                        <  1:3754335/150‑1 (MQ=255)
                                  gATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGc                                                                                                                     >  1:3631824/1‑150 (MQ=255)
                                  gATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGATCTTTATCCCACGGGCGTAGATGTTCGTAACGTGTCACCTGGAGCATCGCCGCATTGCGCTGCTTGAGCAATCGCTGGAGATTGCTCCAGGc                                                                                                                     >  1:2473097/1‑150 (MQ=255)
                                    tACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCgg                                                                                                                   >  2:381657/1‑150 (MQ=255)
                                     aCCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCCGCGCCGCATTGAGCTGCTTGAGCACTCGCTTGAGATTGCGCCAGGCGGt                                                                                                                  >  2:2673274/1‑150 (MQ=255)
                                     aCCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTGGAGATTGCTCCAGGCGGt                                                                                                                  >  1:1063763/1‑150 (MQ=255)
                                      ccGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTg                                                                                                                 <  2:2819996/150‑1 (MQ=255)
                                       cGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGa                                                                                                                >  1:3800173/1‑150 (MQ=255)
                                          cgcTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGaaaa                                                                                                             <  1:547425/150‑1 (MQ=255)
                                             tataCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATc                                                                                                          >  1:1928006/1‑150 (MQ=255)
                                               taCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCg                                                                                                        <  2:3501533/150‑1 (MQ=255)
                                               taCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCg                                                                                                        >  2:1751341/1‑150 (MQ=255)
                                                aCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCgg                                                                                                       >  2:2568875/1‑150 (MQ=255)
                                                aCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCgg                                                                                                       >  1:1452717/1‑150 (MQ=255)
                                                 cTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCggg                                                                                                      <  2:2462730/150‑1 (MQ=255)
                                                 cTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCggg                                                                                                      <  2:4857452/150‑1 (MQ=255)
                                                    cgcgcgCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTc                                                                                                   <  2:868694/150‑1 (MQ=255)
                                                    cgcgcgCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTc                                                                                                   <  2:1539192/150‑1 (MQ=255)
                                                       gcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGtt                                                                                                <  2:2432823/150‑1 (MQ=255)
                                                       gcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGtt                                                                                                <  1:4168780/150‑1 (MQ=255)
                                                           cAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGa                                                                                            <  1:3966391/150‑1 (MQ=255)
                                                              gTGCTGATGTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAg                                                                                         <  2:4667225/150‑1 (MQ=255)
                                                               tGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGc                                                                                        >  1:2286672/1‑150 (MQ=255)
                                                                gCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCa                                                                                       >  2:3608606/1‑150 (MQ=255)
                                                                  tGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTAGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATc                                                                                     <  2:2250663/150‑1 (MQ=255)
                                                                   gATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGGTGTGAAAGCATcc                                                                                    >  2:624213/1‑150 (MQ=255)
                                                                        gCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGt                                                                               >  1:1594358/1‑150 (MQ=255)
                                                                         cTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTc                                                                              >  1:4312820/1‑150 (MQ=255)
                                                                          tCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCg                                                                             <  1:2502602/150‑1 (MQ=255)
                                                                                aGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAAt                                                                       <  2:3700088/150‑1 (MQ=255)
                                                                                  cGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGc                                                                     >  2:1216643/1‑150 (MQ=255)
                                                                                     gATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGctct                                                                  >  1:3059592/1‑150 (MQ=255)
                                                                                     gATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGctct                                                                  >  1:1003228/1‑150 (MQ=255)
                                                                                      aTCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCtctt                                                                 >  1:2417454/1‑150 (MQ=255)
                                                                                         aGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTg                                                              <  2:3348102/150‑1 (MQ=255)
                                                                                           ctctTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTa                                                            >  2:3938423/1‑150 (MQ=255)
                                                                                            tctTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTAt                                                           <  1:729306/150‑1 (MQ=255)
                                                                                               ttATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGAGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTg                                                        >  2:1846636/1‑150 (MQ=255)
                                                                                                tATCCCACGGGCGTCGCTGTTCGTAACGTGGCACGTGGCGCTGCGCCGCATTGCGCTGCTTGAGCCATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGGTCGTTGTGAATGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTgg                                                       <  1:2343698/150‑1 (MQ=255)
                                                                                                     cACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGgccg                                                  >  2:4960216/1‑150 (MQ=255)
                                                                                                        gggCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATGGCGCTGCTTGAGCAATCGCTGGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGACTGCTCTTCTGTATTTGGGGccgccg                                               >  2:3272773/1‑150 (MQ=255)
                                                                                                          gCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGtt                                             <  2:3882686/150‑1 (MQ=255)
                                                                                                          gCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGtt                                             <  1:1384601/150‑1 (MQ=255)
                                                                                                          gCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGtt                                             <  2:3368478/150‑1 (MQ=255)
                                                                                                            gTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGa                                           >  1:3466754/1‑150 (MQ=255)
                                                                                                            gTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGa                                           >  1:1541121/1‑150 (MQ=255)
                                                                                                            gTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCGTCTGTATTTGGGGCCGCCGTTGa                                           >  2:1939211/1‑150 (MQ=255)
                                                                                                             tAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGag                                          >  1:260500/1‑150 (MQ=255)
                                                                                                             tAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGATGCTTGAGCAATCGCTGGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTGGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGACGGTGag                                          >  1:940211/1‑150 (MQ=255)
                                                                                                              aGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagt                                         <  2:1329557/150‑1 (MQ=255)
                                                                                                               gCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagta                                        >  1:3103681/1‑150 (MQ=255)
                                                                                                                 tGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagtaaa                                      >  1:1883050/1‑150 (MQ=255)
                                                                                                                 tGTTCGTAACGTGTCACCTGGCGCAGCGCCGCAGTGAGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGGTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCGGTCTTTGGGGCCGCCGTTGagtaaa                                      >  2:2975087/1‑150 (MQ=255)
                                                                                                                    tCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagtaaagta                                   <  2:573729/150‑1 (MQ=255)
                                                                                                                      gTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAgtaaagtaaa                                 >  1:512410/1‑150 (MQ=255)
                                                                                                                       tAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAc                                <  1:4492758/150‑1 (MQ=255)
                                                                                                                        aaCGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCATCTGTAGTTGGGGCCGCCGTTGAGTAAAGTAAAcc                               >  2:5083205/1‑150 (MQ=255)
                                                                                                                            tgtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTc                           >  2:1771257/1‑150 (MQ=255)
                                                                                                                            tgtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTc                           <  1:4985165/150‑1 (MQ=255)
                                                                                                                             gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCt                          >  1:2547194/1‑150 (MQ=255)
                                                                                                                              tCACCTGGCGCAGCGCCGCCTTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTg                         <  2:1619142/150‑1 (MQ=255)
                                                                                                                                aCCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGc                       <  1:2827098/150‑1 (MQ=255)
                                                                                                                                   tGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGtt                    >  1:3586002/1‑150 (MQ=255)
                                                                                                                                   tGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGtt                    >  1:1015073/1‑150 (MQ=255)
                                                                                                                                         aGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAGAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAAc              >  1:80298/1‑150 (MQ=255)
                                                                                                                                         aGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAAc              >  1:1663621/1‑150 (MQ=255)
                                                                                                                                          gcgcCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACt             >  2:651508/1‑150 (MQ=255)
                                                                                                                                             ccgcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCa          >  1:3193687/1‑150 (MQ=255)
                                                                                                                                             ccgcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCa          >  1:48894/1‑150 (MQ=255)
                                                                                                                                               gcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCAtt        <  2:2355416/150‑1 (MQ=255)
                                                                                                                                                  ttGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGGTATCAACTGCATTGGc     >  2:4069875/1‑150 (MQ=255)
                                                                                                                                                     cgcTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGTATc  <  2:781899/150‑1 (MQ=255)
                                                                                                                                                     cgcTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATc  <  2:4320370/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATC  >  CP009273/3873951‑3874249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑