Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,021,954 C→A 38.8% E298D (GAG→GAT yedI ← DUF808 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,021,9540CA38.8% 308.9 / 242.7 242E298D (GAG→GATyedIDUF808 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (86/62);  new base A (51/43);  total (137/105)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

GGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAATTCCAT  >  CP009273/2021806‑2022099
                                                                                                                                                    |                                                                                                                                                 
ggCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTtc                                                                                                                                                   <  2:2380242/149‑1 (MQ=255)
ggCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTtc                                                                                                                                                   >  2:308033/1‑149 (MQ=255)
  ccgccCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTtatc                                                                                                                                                 <  2:1432479/149‑1 (MQ=255)
     cccTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTAc                                                                                                                                              >  1:1184034/1‑149 (MQ=255)
     cccTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTAc                                                                                                                                              >  2:1099507/1‑149 (MQ=255)
      ccTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACg                                                                                                                                             >  2:3260334/1‑149 (MQ=255)
       ctctTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACgg                                                                                                                                            >  2:1150484/1‑149 (MQ=255)
       ctctTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACgg                                                                                                                                            >  1:2931224/1‑149 (MQ=255)
        tctTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGc                                                                                                                                           >  2:1611829/1‑149 (MQ=255)
           ttAACCTTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCttt                                                                                                                                        <  2:2332555/149‑1 (MQ=255)
           ttAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCttt                                                                                                                                        <  1:1395711/149‑1 (MQ=255)
                  tAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTccc                                                                                                                                 >  2:1130319/1‑149 (MQ=255)
                  tAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTccc                                                                                                                                 <  1:2252299/149‑1 (MQ=255)
                     cTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGc                                                                                                                              >  1:3205881/1‑149 (MQ=255)
                       tCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGcac                                                                                                                            >  1:955170/1‑149 (MQ=255)
                        ccTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTCCGGCTTTCACTCCCAGcacc                                                                                                                           >  2:2977299/1‑149 (MQ=255)
                          tCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGcaccac                                                                                                                         <  1:436394/149‑1 (MQ=255)
                            aCTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGa                                                                                                                       >  2:428463/1‑149 (MQ=255)
                             cTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGAt                                                                                                                      >  2:3159551/1‑149 (MQ=255)
                               aaCCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATgc                                                                                                                    >  1:3140726/1‑149 (MQ=255)
                                aCCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATgcc                                                                                                                   >  2:340756/1‑149 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 >  2:1995478/1‑149 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 >  2:1995798/1‑149 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 <  1:1998383/149‑1 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 >  1:1334089/1‑149 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 >  1:858531/1‑149 (MQ=255)
                                  cATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATgccgc                                                                                                                 >  1:24206/1‑149 (MQ=255)
                                   aTAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATgccgcc                                                                                                                >  1:768522/1‑149 (MQ=255)
                                   aTAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATgccgcc                                                                                                                >  2:3616326/1‑149 (MQ=255)
                                     aaaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCga                                                                                                              <  2:1764987/149‑1 (MQ=255)
                                     aaaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCga                                                                                                              <  1:2333230/149‑1 (MQ=255)
                                     aaaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCga                                                                                                              <  2:1765305/149‑1 (MQ=255)
                                      aaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgat                                                                                                             >  1:385446/1‑149 (MQ=255)
                                      aaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgat                                                                                                             <  1:193529/149‑1 (MQ=255)
                                      aaGCTCATACTAGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgat                                                                                                             <  1:1375135/149‑1 (MQ=255)
                                        gCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgatga                                                                                                           >  1:3584760/1‑149 (MQ=255)
                                          tCATACTCGCCTCCTTGTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATa                                                                                                         <  2:2588263/149‑1 (MQ=255)
                                            aTACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATaaa                                                                                                       <  1:2979147/149‑1 (MQ=255)
                                            aTACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATaaa                                                                                                       <  1:3048990/149‑1 (MQ=255)
                                            aTACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATaaa                                                                                                       >  2:1564494/1‑149 (MQ=255)
                                             tACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                                                                                      >  2:676052/1‑149 (MQ=255)
                                             tACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                                                                                      <  1:1613592/149‑1 (MQ=255)
                                             tACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                                                                                      <  1:3920072/149‑1 (MQ=255)
                                             tACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                                                                                      <  1:3008459/149‑1 (MQ=255)
                                             tACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTAGAAGATGGAGAATTTTGCATAATAGCCAGTGCGAGTATTAATGAGCCTGACCTCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                                                                                      >  2:1001125/1‑149 (MQ=255)
                                              aCTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATc                                                                                                     <  1:3712345/149‑1 (MQ=255)
                                               cTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATcc                                                                                                    >  2:858039/1‑149 (MQ=255)
                                               cTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATcc                                                                                                    <  1:639117/149‑1 (MQ=255)
                                                    ctcctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGaa                                                                                               <  2:3584760/149‑1 (MQ=255)
                                                     tcctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAAt                                                                                              >  2:167987/1‑149 (MQ=255)
                                                       ctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCa                                                                                            <  2:3512434/149‑1 (MQ=255)
                                                        ttttttCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAg                                                                                           >  1:794591/1‑149 (MQ=255)
                                                           tttCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGAtt                                                                                        <  1:3103874/149‑1 (MQ=255)
                                                                aGTGAAACTTGTTCAACTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTaaaa                                                                                   <  1:4021298/149‑1 (MQ=255)
                                                                 gTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAAc                                                                                  >  1:2522754/1‑149 (MQ=255)
                                                                      aCTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCg                                                                             >  1:1740232/1‑149 (MQ=255)
                                                                      aCTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCg                                                                             <  1:3472626/149‑1 (MQ=255)
                                                                        ttGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGt                                                                           <  1:3839661/149‑1 (MQ=255)
                                                                          gTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTaa                                                                         <  2:4000074/149‑1 (MQ=255)
                                                                          gTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTaa                                                                         >  1:1513972/1‑149 (MQ=255)
                                                                          gTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTaa                                                                         <  2:3932106/149‑1 (MQ=255)
                                                                            tCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAata                                                                       >  1:2206261/1‑149 (MQ=255)
                                                                                 ttAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATatcatc                                                                  >  1:259805/1‑149 (MQ=255)
                                                                                     ttGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCa                                                              <  2:2764004/149‑1 (MQ=255)
                                                                                     ttGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCa                                                              <  2:385555/149‑1 (MQ=255)
                                                                                     ttGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCa                                                              <  1:1554042/149‑1 (MQ=255)
                                                                                       gAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCact                                                            >  2:3727621/1‑149 (MQ=255)
                                                                                        aaGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCacta                                                           >  1:517007/1‑149 (MQ=255)
                                                                                          gATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCactact                                                         >  2:2215460/1‑149 (MQ=255)
                                                                                          gATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCactact                                                         >  2:230157/1‑149 (MQ=255)
                                                                                            tGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGc                                                       >  1:1022772/1‑149 (MQ=255)
                                                                                                gAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACtt                                                   <  1:1130319/149‑1 (MQ=255)
                                                                                                gAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACtt                                                   <  1:1127060/149‑1 (MQ=255)
                                                                                                 aaTTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACttt                                                  >  1:1667121/1‑149 (MQ=255)
                                                                                                   ttttGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTtgc                                                <  1:890706/149‑1 (MQ=255)
                                                                                                         ataataGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCg                                          >  1:3074617/1‑149 (MQ=255)
                                                                                                           aataGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGc                                        <  2:2832809/149‑1 (MQ=255)
                                                                                                               gCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGaaa                                    <  1:2790462/149‑1 (MQ=255)
                                                                                                               gCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGaaa                                    <  2:2931077/149‑1 (MQ=255)
                                                                                                               gCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGaaa                                    <  2:621075/149‑1 (MQ=255)
                                                                                                                   gTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGTTTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGtt                                >  1:3130408/1‑149 (MQ=255)
                                                                                                                   gTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGtt                                >  1:3949273/1‑149 (MQ=255)
                                                                                                                    tGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTc                               <  1:3250670/149‑1 (MQ=255)
                                                                                                                    tGCGAGTATTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTc                               >  1:2675095/1‑149 (MQ=255)
                                                                                                                       gAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAAt                            >  2:945673/1‑149 (MQ=255)
                                                                                                                       gAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAAt                            >  2:3329679/1‑149 (MQ=255)
                                                                                                                         gTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATgg                          >  2:1478642/1‑149 (MQ=255)
                                                                                                                           aTTAATGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCg                        >  1:3926748/1‑149 (MQ=255)
                                                                                                                             tAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTg                      <  1:3159558/149‑1 (MQ=255)
                                                                                                                              aaTGTGCCTGACCGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGCAATGTTCAATGGCGTGa                     >  2:3101943/1‑149 (MQ=255)
                                                                                                                                tgtgCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATg                   >  1:983430/1‑149 (MQ=255)
                                                                                                                                  tgCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCa                 >  2:2739003/1‑149 (MQ=255)
                                                                                                                                       gACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAgcggc            >  2:2160996/1‑149 (MQ=255)
                                                                                                                                         ccGCGCATTTTATCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGgcgc          <  1:1055127/149‑1 (MQ=255)
                                                                                                                                           gcgcATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCaa        <  1:3651608/149‑1 (MQ=255)
                                                                                                                                             gcATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAAtt      >  2:2875078/1‑149 (MQ=255)
                                                                                                                                             gcATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAAtt      >  1:3399906/1‑149 (MQ=255)
                                                                                                                                             gcATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAAtt      >  1:1350562/1‑149 (MQ=255)
                                                                                                                                                ttttCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAATTCCa   >  2:3171716/1‑149 (MQ=255)
                                                                                                                                                ttttCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATACAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAATTCCa   >  2:2828881/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACTGCACTTTGCTGCCCGGCGAAATGTTCAATGGCGTGATGCAGCGGCGCAATTCCAT  >  CP009273/2021806‑2022099

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: