Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,999,827 C→A 100% D69E (GAC→GAA ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,999,8270CA97.2% 548.2 / 6.2 177D69E (GAC→GAAygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (86/86);  total (89/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 179 total aligned reads displayed.

AGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCG  >  CP009273/2999682‑2999969
                                                                                                                                                 |                                                                                                                                              
agaAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGaaaa                                                                                                                                              >  1:1140939/1‑148 (MQ=255)
 gaAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAAc                                                                                                                                             >  2:2787512/1‑148 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                            <  2:1358137/148‑1 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                            >  2:1826030/1‑148 (MQ=255)
   aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTTAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                           <  2:2049392/148‑1 (MQ=255)
   aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTTAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                           <  2:2869694/148‑1 (MQ=255)
        tttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAccc                                                                                                                                      <  1:158641/148‑1 (MQ=255)
            tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                  >  1:3521098/1‑148 (MQ=255)
            tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                  >  2:1009394/1‑148 (MQ=255)
            tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                  <  1:308361/148‑1 (MQ=255)
             gTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAcc                                                                                                                                 <  1:296420/148‑1 (MQ=255)
              ttAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAccc                                                                                                                                >  1:2850997/1‑148 (MQ=255)
              ttAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAccc                                                                                                                                >  2:2355657/1‑148 (MQ=255)
              ttAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAccc                                                                                                                                >  2:2354441/1‑148 (MQ=255)
                  cGTGGGAACAGACGCCAGATGAACTGAAACAAGTGCTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                            <  1:1578918/148‑1 (MQ=255)
                  cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                            <  2:257321/148‑1 (MQ=255)
                  cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                            <  2:2184529/148‑1 (MQ=255)
                  cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCAGCtt                                                                                                                            <  1:1182633/148‑1 (MQ=255)
                     gggAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtctc                                                                                                                         >  1:1450509/1‑148 (MQ=255)
                        aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                      <  1:1466291/148‑1 (MQ=255)
                        aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                      <  1:3181923/148‑1 (MQ=255)
                         aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                     >  1:1923732/1‑148 (MQ=255)
                           aGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGc                                                                                                                   >  2:515062/1‑148 (MQ=255)
                             aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                                 <  1:2463821/148‑1 (MQ=255)
                                  aGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCcgcg                                                                                                            >  2:313115/1‑148 (MQ=255)
                                    aTGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCt                                                                                                          <  2:2408677/148‑1 (MQ=255)
                                     tGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCtt                                                                                                         <  1:1049905/148‑1 (MQ=255)
                                      gAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTa                                                                                                        >  2:156348/1‑148 (MQ=255)
                                      gAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTa                                                                                                        >  2:156356/1‑148 (MQ=255)
                                       aaCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTaa                                                                                                       <  1:2355657/148‑1 (MQ=255)
                                        aCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAc                                                                                                      >  1:997035/1‑148 (MQ=255)
                                          tGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCt                                                                                                    >  2:2761776/1‑148 (MQ=255)
                                             aaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt                                                                                                 <  1:3023662/148‑1 (MQ=255)
                                              aCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                                                                >  1:318499/1‑148 (MQ=255)
                                                aaGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCgg                                                                                              <  2:3068090/148‑1 (MQ=255)
                                                  gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                            <  1:664835/148‑1 (MQ=255)
                                                   tACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCt                                                                                           >  1:2414281/1‑148 (MQ=255)
                                                     cTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTg                                                                                         >  2:3300985/1‑148 (MQ=255)
                                                              gCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                                                                <  1:851147/148‑1 (MQ=255)
                                                                 gcagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                             <  2:1687671/148‑1 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            >  1:379278/1‑148 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            >  1:379339/1‑148 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            <  2:1219274/148‑1 (MQ=255)
                                                                    gcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcg                                                                          <  1:2432833/148‑1 (MQ=255)
                                                                     caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                                                                         <  1:3252592/148‑1 (MQ=255)
                                                                      aTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgac                                                                        >  2:1691238/1‑148 (MQ=255)
                                                                       ttAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacg                                                                       >  1:2684601/1‑148 (MQ=255)
                                                                       ttAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacg                                                                       >  1:899028/1‑148 (MQ=255)
                                                                        tAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                                                                      <  1:1709335/148‑1 (MQ=255)
                                                                           aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                   >  1:683466/1‑148 (MQ=255)
                                                                           aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                   >  2:492904/1‑148 (MQ=255)
                                                                           aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                   <  1:2327868/148‑1 (MQ=255)
                                                                           aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                   >  1:34066/1‑148 (MQ=255)
                                                                              cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaa                                                                >  2:402898/1‑148 (MQ=255)
                                                                               aCTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaaa                                                               <  1:3028916/148‑1 (MQ=255)
                                                                                cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                              >  2:2897846/1‑148 (MQ=255)
                                                                                cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                              >  1:104156/1‑148 (MQ=255)
                                                                                      gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaaa                                                        <  2:1714175/148‑1 (MQ=255)
                                                                                       cTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAAt                                                       <  1:2012046/148‑1 (MQ=255)
                                                                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATc                                                      >  1:3038316/1‑148 (MQ=255)
                                                                                         gAAAACATCTCAACCAAAGTCTTTAAGAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                     <  2:2193196/148‑1 (MQ=255)
                                                                                           aaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCt                                                   >  2:1883811/1‑148 (MQ=255)
                                                                                               aTCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCAct                                               >  2:2742175/1‑147 (MQ=255)
                                                                                               aTCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg                                               >  1:2959226/1‑148 (MQ=255)
                                                                                                tctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACgg                                              <  2:1934321/148‑1 (MQ=255)
                                                                                                  tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCg                                            >  2:2464955/1‑148 (MQ=255)
                                                                                                     aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGaaa                                         >  2:3522428/1‑148 (MQ=255)
                                                                                                       cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                       >  2:1882841/1‑148 (MQ=255)
                                                                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtg                                    >  2:633341/1‑148 (MQ=255)
                                                                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtg                                    >  2:1723902/1‑148 (MQ=255)
                                                                                                           gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACcgtgc                                   <  2:1276108/148‑1 (MQ=255)
                                                                                                            tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcg                                  <  2:354653/148‑1 (MQ=255)
                                                                                                             cTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                 <  2:2561101/148‑1 (MQ=255)
                                                                                                              tttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgtg                                >  2:1991244/1‑148 (MQ=255)
                                                                                                              tttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgtg                                <  2:3412385/148‑1 (MQ=255)
                                                                                                              tttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACcgtgcgtg                                >  1:1610346/1‑148 (MQ=255)
                                                                                                               ttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGa                               <  2:567469/148‑1 (MQ=255)
                                                                                                                    aGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCg                          >  1:624056/1‑148 (MQ=255)
                                                                                                                    aGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCg                          >  1:729715/1‑148 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                        <  1:3051582/148‑1 (MQ=255)
                                                                                                                        gATTAGGCATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                      >  2:1109329/1‑148 (MQ=255)
                                                                                                                          ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCTCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAata                    <  2:951199/148‑1 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                   >  1:2414093/1‑148 (MQ=255)
                                                                                                                            aGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATg                  <  1:2297622/148‑1 (MQ=255)
                                                                                                                             ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGa                 <  1:1619136/148‑1 (MQ=255)
                                                                                                                              gTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAt                <  2:3499201/148‑1 (MQ=255)
                                                                                                                                aTTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtct              <  2:2414281/148‑1 (MQ=255)
                                                                                                                                 tttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc             >  1:280911/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            >  1:2783409/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            >  1:162973/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            <  2:1428742/148‑1 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            >  1:571492/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            <  2:1429685/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt          >  2:3246626/1‑148 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt          <  2:356965/148‑1 (MQ=255)
                                                                                                                                     cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc         <  1:515062/148‑1 (MQ=255)
                                                                                                                                        aTTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgc      >  1:557852/1‑148 (MQ=255)
                                                                                                                                          tCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcgc    <  1:3076364/148‑1 (MQ=255)
                                                                                                                                           ccGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGcc   <  2:3518267/148‑1 (MQ=255)
                                                                                                                                            cgcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCg  <  2:44504/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCG  >  CP009273/2999682‑2999969

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: