Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,626,858 T→A 5.5% H191L (CAC→CTC)  guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,626,8580TA 5.5% 928.2 / 32.7 289H191L (CAC→CTC) guaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (141/132);  new base A (6/10);  total (147/142)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 295 total aligned reads displayed.

TGCGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGTCGGTAAC  >  CP009273/2626710‑2627006
                                                                                                                                                    |                                                                                                                                                    
tGCGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGt                                                                                                                                                      <  1:3713907/149‑1 (MQ=255)
 gCGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTg                                                                                                                                                     <  2:2725418/149‑1 (MQ=255)
 gCGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTg                                                                                                                                                     >  1:2431159/1‑149 (MQ=255)
  cGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTgg                                                                                                                                                    <  2:3270281/149‑1 (MQ=255)
  cGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTgg                                                                                                                                                    <  1:2551013/149‑1 (MQ=255)
    gCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGaa                                                                                                                                                  <  1:1978666/149‑1 (MQ=255)
       aCCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACTGTGATCATGCCGATCAGGTGGAATTc                                                                                                                                               >  2:2262223/1‑149 (MQ=255)
         cAGCGCGTCAACCCGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGt                                                                                                                                             <  2:3588789/149‑1 (MQ=255)
          aGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGCATTCGtc                                                                                                                                            >  2:2645347/1‑149 (MQ=255)
          aGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtc                                                                                                                                            >  1:2718422/1‑149 (MQ=255)
           gcgcGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtca                                                                                                                                           <  2:1049753/149‑1 (MQ=255)
           gcgcGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtca                                                                                                                                           <  2:1050218/149‑1 (MQ=255)
             gcgTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTACTGGAAGTCTTTCACGGTGATCAAGCCGATCAGGTGGAATTCGtcctc                                                                                                                                         >  1:2077517/1‑149 (MQ=255)
               gTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCaa                                                                                                                                       >  1:135045/1‑149 (MQ=255)
                tCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAac                                                                                                                                      <  1:391294/149‑1 (MQ=255)
                   acacGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAaccac                                                                                                                                   >  1:4712547/1‑149 (MQ=255)
                    cacGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAaccacc                                                                                                                                  <  1:4407760/149‑1 (MQ=255)
                     acGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAAccacca                                                                                                                                 <  1:3339332/149‑1 (MQ=255)
                      cGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAg                                                                                                                                >  1:1812618/1‑149 (MQ=255)
                        ctctTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcg                                                                                                                              >  2:423277/1‑149 (MQ=255)
                         tctTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGGTTACGGTCCGCGGTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCCTCAACCACCCgcgc                                                                                                                             >  1:3606037/1‑149 (MQ=255)
                          ctTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCt                                                                                                                            <  1:1421547/149‑1 (MQ=255)
                           ttCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtt                                                                                                                           >  2:3680913/1‑149 (MQ=255)
                           ttCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtt                                                                                                                           <  1:1066837/149‑1 (MQ=255)
                            tCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCttt                                                                                                                          >  2:94281/1‑149 (MQ=255)
                              gTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttttt                                                                                                                        >  2:2488697/1‑149 (MQ=255)
                              gTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttttt                                                                                                                        <  2:3298364/149‑1 (MQ=255)
                               ttACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTc                                                                                                                       <  2:3762350/149‑1 (MQ=255)
                               ttACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTc                                                                                                                       <  1:2940331/149‑1 (MQ=255)
                                  cccGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAAc                                                                                                                    <  2:4314948/149‑1 (MQ=255)
                                      cACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGt                                                                                                                <  2:4676399/149‑1 (MQ=255)
                                      cACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGt                                                                                                                <  2:518182/149‑1 (MQ=255)
                                         cTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGtttt                                                                                                             >  1:343787/1‑149 (MQ=255)
                                         cTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGtttt                                                                                                             <  1:721887/149‑1 (MQ=255)
                                          tGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGttttt                                                                                                            <  2:2471446/149‑1 (MQ=255)
                                          tGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGttttt                                                                                                            <  1:3162726/149‑1 (MQ=255)
                                                aaCCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCa                                                                                                      <  2:1738327/149‑1 (MQ=255)
                                                   cGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt                                                                                                   <  1:3081389/149‑1 (MQ=255)
                                                     cTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTGTTCAACGCGTTTTTCGTGCAttttt                                                                                                 <  2:58246/149‑1 (MQ=255)
                                                      tGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTAAGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTAAACGCGTTTTTCGGGCATTTTTg                                                                                                >  1:137653/1‑149 (MQ=255)
                                                        cACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCCACGCGTTTTTCGTGCATTTTTGcc                                                                                              <  2:3211092/149‑1 (MQ=255)
                                                        cACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGcc                                                                                              <  2:3871998/149‑1 (MQ=255)
                                                                  cAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACtt                                                                                    >  2:1174428/1‑149 (MQ=255)
                                                                    gACGGCCTTGCTCGTCTTTACAGGCGTTCGGTGTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCa                                                                                  <  2:3917516/149‑1 (MQ=255)
                                                                     aCGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCAc                                                                                 <  1:3030124/149‑1 (MQ=255)
                                                                      cGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGTAAGTATTTAACGGTGATCATGCCGATAAGGTGGAATTCGTCATCAACCAACCGCGCATTTTCAACACGTTTTTAGTGCATTTTTTCAAGAACCACAACACg                                                                                >  1:2209754/1‑149 (MQ=255)
                                                                      cGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACg                                                                                <  1:1970595/149‑1 (MQ=255)
                                                                      cGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACg                                                                                <  1:1732826/149‑1 (MQ=255)
                                                                       ggCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTTCCAGCACCACATCACgg                                                                               >  1:3818138/1‑149 (MQ=255)
                                                                       ggCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACgg                                                                               >  2:779023/1‑149 (MQ=255)
                                                                        gCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAAAGCGTTTTTCGTGCCTTTTTGCCAGCACCACATAACggg                                                                              >  1:2228405/1‑149 (MQ=255)
                                                                        gCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCAAGCCGATAAGGTGGAAATCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAATTTTTCCAGCACCACTTAACcgg                                                                              >  1:646714/1‑149 (MQ=255)
                                                                           ttGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGctt                                                                           >  2:2104265/1‑149 (MQ=255)
                                                                            tGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGcttc                                                                          >  1:4238971/1‑149 (MQ=255)
                                                                             gCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTGTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGcttca                                                                         >  2:4426285/1‑149 (MQ=255)
                                                                                cGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTAAAAGGGTTTTTAGTGCATTTTTTCCAGCACCACATCACGGGcttcacct                                                                      >  1:4653001/1‑149 (MQ=255)
                                                                                 gTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACATCACGGGcttcaccat                                                                     >  1:1006437/1‑147 (MQ=255)
                                                                                  tCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGcttcaccttc                                                                    <  2:2405883/149‑1 (MQ=255)
                                                                                       aCAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTcacgca                                                               <  2:3116327/149‑1 (MQ=255)
                                                                                             gTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGa                                                         <  2:293735/149‑1 (MQ=255)
                                                                                             gTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGa                                                         <  2:4570809/149‑1 (MQ=255)
                                                                                              ttCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGAc                                                        >  2:4721642/1‑149 (MQ=255)
                                                                                              ttCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCATCACGCACGGTGAc                                                        >  1:2115746/1‑149 (MQ=255)
                                                                                                  gTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGa                                                    <  2:1146083/149‑1 (MQ=255)
                                                                                                      aCGTTCCGCTTTCTGGAAGTCTTTAACGGTGATCATGCCGATAAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTTCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACCCACGGTGACCAGACGct                                                >  1:1750484/1‑149 (MQ=255)
                                                                                                        gTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCCACGCGTTTTTAGTGCATTTTTTCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGctct                                              >  1:589605/1‑149 (MQ=255)
                                                                                                        gTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGctct                                              >  1:1337766/1‑149 (MQ=255)
                                                                                                         ttCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCtt                                             >  2:2015197/1‑149 (MQ=255)
                                                                                                           ccGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTc                                           <  2:3340723/149‑1 (MQ=255)
                                                                                                             gCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCgg                                         <  1:1697231/149‑1 (MQ=255)
                                                                                                               tttCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCg                                       >  1:1615246/1‑149 (MQ=255)
                                                                                                               tttCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCg                                       >  1:1454353/1‑149 (MQ=255)
                                                                                                               tttCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCg                                       >  2:3059985/1‑149 (MQ=255)
                                                                                                                 tCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTc                                     <  2:2695409/149‑1 (MQ=255)
                                                                                                                 tCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTc                                     <  1:1792323/149‑1 (MQ=255)
                                                                                                                   tGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCAt                                   >  1:3258474/1‑149 (MQ=255)
                                                                                                                    ggAAGTCTTTCACGGTGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATg                                  >  1:2342761/1‑149 (MQ=255)
                                                                                                                      aaGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTa                                >  1:1275024/1‑149 (MQ=255)
                                                                                                                       aGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTaa                               >  2:3700776/1‑149 (MQ=255)
                                                                                                                        gTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTaaa                              >  2:1409544/1‑149 (MQ=255)
                                                                                                                           tttCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGc                           <  2:3430496/149‑1 (MQ=255)
                                                                                                                             tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTa                         >  1:4000432/1‑149 (MQ=255)
                                                                                                                              cACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTaa                        >  1:3002329/1‑149 (MQ=255)
                                                                                                                               aCGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAAc                       >  2:6704/1‑149 (MQ=255)
                                                                                                                               aCGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAAc                       >  2:4589060/1‑149 (MQ=255)
                                                                                                                                  gTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCgg                    >  2:3498042/1‑149 (MQ=255)
                                                                                                                                   tGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGc                   <  2:691473/149‑1 (MQ=255)
                                                                                                                                   tGATCATGCCGATCAGGAGGAATTCGTCATCAACCACCAGCGCTGTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGc                   <  2:726218/149‑1 (MQ=255)
                                                                                                                                    gATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCt                  >  2:1119024/1‑149 (MQ=255)
                                                                                                                                     aTCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTg                 >  1:3007930/1‑149 (MQ=255)
                                                                                                                                     aTCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTg                 >  1:3007911/1‑149 (MQ=255)
                                                                                                                                      tCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTgg                >  1:1561015/1‑149 (MQ=255)
                                                                                                                                       cATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGt               >  2:955584/1‑149 (MQ=255)
                                                                                                                                        aTGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGtt              >  1:619530/1‑149 (MQ=255)
                                                                                                                                          gCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCa            <  2:4507363/149‑1 (MQ=255)
                                                                                                                                            cGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAgg          >  1:1176353/1‑149 (MQ=255)
                                                                                                                                             gATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGt         <  2:1617902/149‑1 (MQ=255)
                                                                                                                                                 aGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGTCGGt     <  1:933689/149‑1 (MQ=255)
                                                                                                                                                  ggTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGTCGGTa    <  1:2857699/149‑1 (MQ=255)
                                                                                                                                                    tGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGTCGGTAAc  <  1:803537/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGGCAACCAGCGCGTCAACACGCTCTTCGTTACCCGCACCTGCGCCAACCGCTGCACCAACACGCAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCACCACTTCACGGGCTTCACCTTCACGCACGGTGACCAGACGCTCTTTCGGCGTCATGTAAACGCTAACCGGCTGGTTCAGGTCGGTAAC  >  CP009273/2626710‑2627006

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