Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,088,235 C→A 5.1% A69S (GCC→TCC)  yggR ← putative PilT family AAA+ ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,088,2350CA 5.1% 1035.1 / 21.8 271A69S (GCC→TCC) yggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (125/132);  new base A (6/8);  total (131/140)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 280 total aligned reads displayed.

GCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGT  >  CP009273/3088086‑3088384
                                                                                                                                                     |                                                                                                                                                     
gCTCTTGAGTGATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGc                                                                                                                                                       <  2:102121/150‑1 (MQ=255)
gCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGc                                                                                                                                                       <  2:3090987/150‑1 (MQ=255)
gCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGc                                                                                                                                                       >  1:1006618/1‑150 (MQ=255)
 ctctTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCa                                                                                                                                                      >  2:3611811/1‑150 (MQ=255)
 ctctTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCa                                                                                                                                                      <  2:1903667/150‑1 (MQ=255)
 ctctTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCa                                                                                                                                                      >  2:1343513/1‑150 (MQ=255)
  tctTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCac                                                                                                                                                     >  2:4101968/1‑149 (MQ=255)
    ttGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCaaaa                                                                                                                                                   >  1:1007309/1‑150 (MQ=255)
    ttGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCaaaa                                                                                                                                                   >  2:205446/1‑150 (MQ=255)
       aGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACTCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATcc                                                                                                                                                <  1:2500333/150‑1 (MQ=255)
       aGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATcc                                                                                                                                                >  2:1519669/1‑150 (MQ=255)
       aGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATcc                                                                                                                                                <  1:1828792/150‑1 (MQ=255)
             tCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGa                                                                                                                                          >  2:2297806/1‑150 (MQ=255)
             tCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGa                                                                                                                                          <  1:4210054/150‑1 (MQ=255)
              ccGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGAc                                                                                                                                         <  2:1481891/150‑1 (MQ=255)
               cGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGAcc                                                                                                                                        <  2:3758964/150‑1 (MQ=255)
               cGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGAcc                                                                                                                                        <  1:2844929/150‑1 (MQ=255)
               cGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGAcc                                                                                                                                        <  1:3564053/150‑1 (MQ=255)
                tgCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCa                                                                                                                                       <  2:118596/149‑1 (MQ=255)
                ggCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCa                                                                                                                                       <  2:727593/150‑1 (MQ=255)
                ggCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCa                                                                                                                                       >  2:4318010/1‑150 (MQ=255)
                ggCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCa                                                                                                                                       >  2:1837653/1‑150 (MQ=255)
                 gCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCAt                                                                                                                                      <  2:912195/150‑1 (MQ=255)
                  cAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGATCAtt                                                                                                                                     >  2:2246573/1‑150 (MQ=255)
                  cAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCAtt                                                                                                                                     <  1:234422/150‑1 (MQ=255)
                   aaTACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATtc                                                                                                                                    <  1:961626/150‑1 (MQ=255)
                   aaTACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATtc                                                                                                                                    <  1:4178051/150‑1 (MQ=255)
                     tACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATtctc                                                                                                                                  <  1:1360935/150‑1 (MQ=255)
                     tACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATtctc                                                                                                                                  <  1:638932/150‑1 (MQ=255)
                      aCCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTcc                                                                                                                                 >  2:2650346/1‑150 (MQ=255)
                      aCCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCCTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTcc                                                                                                                                 >  2:2651441/1‑150 (MQ=255)
                       ccGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCa                                                                                                                                <  2:4488375/150‑1 (MQ=255)
                            ggTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAat                                                                                                                           >  2:438288/1‑150 (MQ=255)
                             gTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAata                                                                                                                          >  2:1080617/1‑150 (MQ=255)
                               gcgcCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATAtt                                                                                                                        <  1:2918856/150‑1 (MQ=255)
                                 gcCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGc                                                                                                                      >  1:2590134/1‑150 (MQ=255)
                                   cAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGccc                                                                                                                    >  2:3308163/1‑150 (MQ=255)
                                    aaGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGACAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCg                                                                                                                   >  2:1802824/1‑150 (MQ=255)
                                     aGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGTTGACCATTCTCCAGCAATATTGCCCGc                                                                                                                  <  1:1146865/150‑1 (MQ=255)
                                     aGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGc                                                                                                                  >  2:3769794/1‑150 (MQ=255)
                                      gCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCt                                                                                                                 <  1:346750/150‑1 (MQ=255)
                                      gCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCAGTAACGCCAGCGTAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCt                                                                                                                 <  1:3990805/150‑1 (MQ=255)
                                       cTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGACATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTg                                                                                                                >  1:4466305/1‑150 (MQ=255)
                                       cTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTg                                                                                                                >  1:2229632/1‑150 (MQ=255)
                                       cTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTg                                                                                                                >  1:962377/1‑150 (MQ=255)
                                        tGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGa                                                                                                               <  2:1488528/150‑1 (MQ=255)
                                        tGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGa                                                                                                               <  2:1571795/150‑1 (MQ=255)
                                        tGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGa                                                                                                               <  2:3198835/150‑1 (MQ=255)
                                          ttCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGAAAAATCCAGCTGACCATTCTCCCGCAATATTGCCCGCTGATc                                                                                                             <  2:4234383/150‑1 (MQ=255)
                                          ttCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGAAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATc                                                                                                             <  2:613116/150‑1 (MQ=255)
                                                 tGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATcc                                                                                                      <  2:3297385/150‑1 (MQ=255)
                                                    gggCAGTTCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGc                                                                                                   >  1:3590502/1‑150 (MQ=255)
                                                    gggCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGc                                                                                                   <  2:2875723/150‑1 (MQ=255)
                                                             gAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGc                                                                                          >  1:4363014/1‑150 (MQ=255)
                                                              aaGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCa                                                                                         >  2:1699534/1‑150 (MQ=255)
                                                               aGGTAACCGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAg                                                                                        <  2:3150288/150‑1 (MQ=255)
                                                                     cAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGctc                                                                                  >  1:2437522/1‑150 (MQ=255)
                                                                          gTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGa                                                                             >  2:339515/1‑150 (MQ=255)
                                                                              cGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTc                                                                         <  2:3472658/150‑1 (MQ=255)
                                                                              cGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTc                                                                         <  2:3050050/150‑1 (MQ=255)
                                                                              cGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTc                                                                         <  1:4425574/150‑1 (MQ=255)
                                                                               gCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTcc                                                                        <  1:3096821/150‑1 (MQ=255)
                                                                               gCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTcc                                                                        <  2:4516059/150‑1 (MQ=255)
                                                                                     gAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTa                                                                  >  2:2567250/1‑150 (MQ=255)
                                                                                     gAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTa                                                                  <  1:3300903/150‑1 (MQ=255)
                                                                                     gAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGAAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTa                                                                  <  1:1328797/150‑1 (MQ=255)
                                                                                        atgccatgccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCa                                                               >  2:2512516/1‑150 (MQ=255)
                                                                                         tgccatgccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCaa                                                              >  1:1235048/1‑150 (MQ=255)
                                                                                            catgccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTCGCTCTTCGACGTCCGGCGGATCAAACg                                                           >  2:4342307/1‑150 (MQ=255)
                                                                                            catgccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACg                                                           >  2:2370711/1‑150 (MQ=255)
                                                                                                      cgcgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTc                                                 <  1:3394064/150‑1 (MQ=255)
                                                                                                        cgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCa                                               <  1:802258/150‑1 (MQ=255)
                                                                                                          aaTGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCAtt                                             >  2:2919185/1‑150 (MQ=255)
                                                                                                            tGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCt                                           >  2:1620577/1‑150 (MQ=255)
                                                                                                                  gCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCcgcg                                     >  2:3278755/1‑150 (MQ=255)
                                                                                                                     aCGCAATCGCTGGTTTTCCGCCAGCGACACGGAAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAAt                                  <  1:183717/150‑1 (MQ=255)
                                                                                                                      cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATg                                 >  1:3358399/1‑150 (MQ=255)
                                                                                                                      cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATg                                 >  1:896849/1‑150 (MQ=255)
                                                                                                                      cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATg                                 >  1:322034/1‑150 (MQ=255)
                                                                                                                      cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATg                                 <  1:4232816/150‑1 (MQ=255)
                                                                                                                       gCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgc                                <  2:814576/150‑1 (MQ=255)
                                                                                                                       gCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgc                                >  2:4204123/1‑150 (MQ=255)
                                                                                                                       gCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgc                                <  2:3103291/150‑1 (MQ=255)
                                                                                                                        cAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgcg                               <  1:3884364/150‑1 (MQ=255)
                                                                                                                        cAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgcg                               <  2:1852277/150‑1 (MQ=255)
                                                                                                                         aaTCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATgcgc                              <  1:4112922/150‑1 (MQ=255)
                                                                                                                                ggTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTg                       >  2:1063955/1‑150 (MQ=255)
                                                                                                                                  ttttCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATGAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCg                     <  2:1116032/150‑1 (MQ=255)
                                                                                                                                  ttttCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCg                     <  2:4029822/150‑1 (MQ=255)
                                                                                                                                    ttCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCggg                   >  2:2807215/1‑150 (MQ=255)
                                                                                                                                          cAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCGAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGgc             <  2:3902666/150‑1 (MQ=255)
                                                                                                                                          cAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGgc             <  1:1491455/150‑1 (MQ=255)
                                                                                                                                          cAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGgc             <  2:3509298/150‑1 (MQ=255)
                                                                                                                                           aGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGgcg            >  2:4363991/1‑150 (MQ=255)
                                                                                                                                              gACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTg         >  1:1351087/1‑150 (MQ=255)
                                                                                                                                              gACACGGCAAAATCCAGCTGACCATTATCAAGCAATATTGCACGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTACGGCGTATCAAACGGCGCAGCTTACATTCTGCCGCGAATGAGCCATCGTGCGGGCCAGGCGCTg         >  1:3641915/1‑150 (MQ=255)
                                                                                                                                               acacGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGc        >  1:2085210/1‑150 (MQ=255)
                                                                                                                                                 acGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGcac      >  2:2090197/1‑150 (MQ=255)
                                                                                                                                                  cGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGcaca     >  1:2801252/1‑150 (MQ=255)
                                                                                                                                                     cAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt  >  2:749973/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGT  >  CP009273/3088086‑3088384

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: