Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 334,250 G→T 19.5% R72L (CGT→CTT)  yahG → DUF1116 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273334,2500GT19.5% 962.9 / 168.1 332R72L (CGT→CTT) yahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base G (120/146);  new base T (33/33);  total (153/179)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

ACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAGCG  >  CP009273/334103‑334398
                                                                                                                                                   |                                                                                                                                                    
 ccTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGta                                                                                                                                                    <  2:523372/149‑1 (MQ=255)
 ccTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGta                                                                                                                                                    >  1:1649881/1‑149 (MQ=255)
  cTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGtat                                                                                                                                                   <  1:3389184/149‑1 (MQ=255)
   tGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGtatt                                                                                                                                                  <  1:1957110/149‑1 (MQ=255)
    gAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGtatta                                                                                                                                                 >  1:3553516/1‑149 (MQ=255)
     aaaaaGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGtattat                                                                                                                                                >  1:345788/1‑149 (MQ=255)
      aaaaGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATc                                                                                                                                               <  1:1774640/149‑1 (MQ=255)
        aaGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCa                                                                                                                                             <  2:1080154/149‑1 (MQ=255)
        aaGCAGCATGTAGAAGGGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCa                                                                                                                                             <  1:4337984/149‑1 (MQ=255)
          gcagcaTGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGt                                                                                                                                           <  2:1285922/149‑1 (MQ=255)
          gcagcaTGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGt                                                                                                                                           >  1:111180/1‑149 (MQ=255)
               aTGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCAt                                                                                                                                      >  2:207350/1‑149 (MQ=255)
                tGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTGATTATCCAGTCGCATc                                                                                                                                     >  1:2532541/1‑149 (MQ=255)
                 gTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATcc                                                                                                                                    <  2:169661/149‑1 (MQ=255)
                   aGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCgg                                                                                                                                  <  1:1691404/149‑1 (MQ=255)
                   aGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCgg                                                                                                                                  <  1:4551524/149‑1 (MQ=255)
                      aGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGc                                                                                                                               >  1:1650856/1‑149 (MQ=255)
                      aGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGc                                                                                                                               >  1:1190503/1‑149 (MQ=255)
                       ttGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGTCAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCt                                                                                                                              <  1:2171256/148‑1 (MQ=255)
                         gACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGa                                                                                                                            >  1:935452/1‑149 (MQ=255)
                         gACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGa                                                                                                                            <  1:4477655/149‑1 (MQ=255)
                          aCTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGAt                                                                                                                           <  1:2684875/149‑1 (MQ=255)
                          aCTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGAt                                                                                                                           >  2:3561094/1‑149 (MQ=255)
                           cTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGAtt                                                                                                                          <  1:4200424/149‑1 (MQ=255)
                            tCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTg                                                                                                                         >  1:644303/1‑149 (MQ=255)
                             cAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTgg                                                                                                                        >  1:1526506/1‑149 (MQ=255)
                               aCTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGtt                                                                                                                      >  2:361504/1‑149 (MQ=255)
                                cTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGttt                                                                                                                     >  1:409070/1‑149 (MQ=255)
                                      tGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCa                                                                                                               >  2:1301394/1‑149 (MQ=255)
                                        gACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAgg                                                                                                             <  1:3782866/149‑1 (MQ=255)
                                         aCGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGc                                                                                                            >  1:3583156/1‑149 (MQ=255)
                                         aCGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGc                                                                                                            <  2:139152/149‑1 (MQ=255)
                                           gccgccGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGa                                                                                                          <  2:3969012/149‑1 (MQ=255)
                                              tccgGGGCAGGGCAATATGCAGGTGGTTCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCAGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTa                                                                                                       <  1:4157415/148‑1 (MQ=255)
                                              gccgGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTa                                                                                                       >  2:13762/1‑149 (MQ=255)
                                                 ggggCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACg                                                                                                    <  1:2302300/149‑1 (MQ=255)
                                                  gggcagggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgt                                                                                                   >  2:1848711/1‑149 (MQ=255)
                                                   ggcagggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgtg                                                                                                  <  2:3466939/149‑1 (MQ=255)
                                                    gcagggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCCGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgtgg                                                                                                 <  2:2181891/149‑1 (MQ=255)
                                                     cagggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgtggt                                                                                                <  1:1294312/149‑1 (MQ=255)
                                                      agggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgtggtg                                                                                               >  1:861387/1‑149 (MQ=255)
                                                        ggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGcc                                                                                             >  1:3033559/1‑149 (MQ=255)
                                                         gcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCAGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGACg                                                                                            >  2:3933135/1‑149 (MQ=255)
                                                              atGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCAt                                                                                       >  2:4540055/1‑149 (MQ=255)
                                                                gCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGa                                                                                     >  1:4670803/1‑149 (MQ=255)
                                                                 cAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGAc                                                                                    <  2:1611383/149‑1 (MQ=255)
                                                                  aGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACg                                                                                   <  2:4100724/149‑1 (MQ=255)
                                                                  aGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACg                                                                                   <  2:1151784/149‑1 (MQ=255)
                                                                       gtgCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaa                                                                              >  2:515647/1‑149 (MQ=255)
                                                                       gtgCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaa                                                                              >  2:3605530/1‑149 (MQ=255)
                                                                         gCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaaaa                                                                            <  1:3281289/149‑1 (MQ=255)
                                                                         gCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaaaa                                                                            <  1:1601897/149‑1 (MQ=255)
                                                                         gCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTGTGATCAGGCGATTAACGGGGTGCCGGGCAGGACGGCGaaaa                                                                            >  1:617294/1‑149 (MQ=255)
                                                                         gCAGGCGCTGGATAACAGTGCCGAGTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaaaa                                                                            <  1:4425440/149‑1 (MQ=255)
                                                                           aGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAAcc                                                                          <  1:1207620/149‑1 (MQ=255)
                                                                           aGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAAcc                                                                          <  2:3713806/149‑1 (MQ=255)
                                                                            ggCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCa                                                                         <  1:201964/149‑1 (MQ=255)
                                                                            ggCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCa                                                                         <  1:1267540/149‑1 (MQ=255)
                                                                             gcgcTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCAt                                                                        <  2:1892153/149‑1 (MQ=255)
                                                                               gcTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCAttc                                                                      <  1:2865643/149‑1 (MQ=255)
                                                                                 tGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCAttctt                                                                    <  1:1823248/149‑1 (MQ=255)
                                                                                   gATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                  <  1:4664550/149‑1 (MQ=255)
                                                                                   gATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                  >  1:3129347/1‑149 (MQ=255)
                                                                                    aTAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCAc                                                                 <  2:2405030/149‑1 (MQ=255)
                                                                                    aTAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCAc                                                                 <  1:1473682/149‑1 (MQ=255)
                                                                                      aaCATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGc                                                               >  2:4256763/1‑149 (MQ=255)
                                                                                       aCATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGcc                                                              >  2:813901/1‑149 (MQ=255)
                                                                                          ttGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCggg                                                           <  2:1137287/149‑1 (MQ=255)
                                                                                            gCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgcc                                                         >  1:3844461/1‑149 (MQ=255)
                                                                                            gCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgcc                                                         >  2:1522933/1‑149 (MQ=255)
                                                                                              cGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgccgc                                                       <  2:636205/149‑1 (MQ=255)
                                                                                              cGATGCGCCGCTGGCGGACAAACTCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGAGTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgccgc                                                       <  1:791439/149‑1 (MQ=255)
                                                                                               gATTCGCCGCTGGCGGACAAACTCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgccgcc                                                      <  1:2773991/149‑1 (MQ=255)
                                                                                               gATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGgccgcc                                                      <  1:4178305/149‑1 (MQ=255)
                                                                                                    gccgcTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCa                                                 >  1:2566987/1‑149 (MQ=255)
                                                                                                     ccgcTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCAc                                                >  1:3935996/1‑149 (MQ=255)
                                                                                                     ccgcTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGAGTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCAc                                                >  2:3831949/1‑149 (MQ=255)
                                                                                                      cgcTGGCGGACAAACTCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCTGTCATGACGGCGATTACCATTCTTCACGCCGGGCCGCCGATCAcc                                               <  1:3015173/149‑1 (MQ=255)
                                                                                                      cgcTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCAcc                                               <  1:1740001/149‑1 (MQ=255)
                                                                                                                 aaaaTCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAAt                                    >  2:4440362/1‑149 (MQ=255)
                                                                                                                  aaaTCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAAtg                                   <  2:1550535/149‑1 (MQ=255)
                                                                                                                       gccgccgcTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTgcgg                              <  2:1497637/149‑1 (MQ=255)
                                                                                                                       gccgccgcTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTgcgg                              <  1:1417323/149‑1 (MQ=255)
                                                                                                                        ccgccgcTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTgcggc                             >  2:1190022/1‑149 (MQ=255)
                                                                                                                        ccgccgcTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTgcggc                             >  1:3765516/1‑149 (MQ=255)
                                                                                                                          gccgcTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGgcgc                           >  2:3214995/1‑149 (MQ=255)
                                                                                                                            cgcTAACCAGCAGGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGa                         <  2:700872/149‑1 (MQ=255)
                                                                                                                                aaCCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGaa                     >  2:3308615/1‑149 (MQ=255)
                                                                                                                                 aCCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGaaa                    <  1:531906/149‑1 (MQ=255)
                                                                                                                                  ccAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAg                   >  2:2859975/1‑149 (MQ=255)
                                                                                                                                  ccAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAg                   <  2:2532037/149‑1 (MQ=255)
                                                                                                                                  ccAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAg                   >  2:4590647/1‑149 (MQ=255)
                                                                                                                                   cagcagGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAgg                  >  1:1134441/1‑149 (MQ=255)
                                                                                                                                    agcagGCGCTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGgc                 >  1:2905205/1‑149 (MQ=255)
                                                                                                                                           gcTGGAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCa          <  2:2727862/149‑1 (MQ=255)
                                                                                                                                            cTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCAc         <  2:2560098/149‑1 (MQ=255)
                                                                                                                                            cTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCAc         <  1:2506513/149‑1 (MQ=255)
                                                                                                                                             tGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCAcc        >  2:4329304/1‑149 (MQ=255)
                                                                                                                                             tGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCAcc        >  1:4415906/1‑149 (MQ=255)
                                                                                                                                              ggAGCTTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCg       >  2:1988892/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAGCG  >  CP009273/334103‑334398

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: