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.7% 682.5 / 6.1 215D69E (GAC→GAAygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (103/107);  total (106/109)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 215 total aligned reads displayed.

AAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCG  >  CP009273/2999684‑2999966
                                                                                                                                               |                                                                                                                                           
aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                         >  1:2410890/1‑148 (MQ=255)
aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                         <  2:708693/148‑1 (MQ=255)
 aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                        >  1:1363617/1‑148 (MQ=255)
 aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                        >  1:238071/1‑148 (MQ=255)
  aGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAGGGTCTGTAATAGTGGCTTAGGTCTTTGCGGCTTGCGCGAAATCTcc                                                                                                                                       >  2:3237007/1‑148 (MQ=255)
  aGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTcc                                                                                                                                       >  1:3042238/1‑148 (MQ=255)
  aGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTcc                                                                                                                                       >  2:3405553/1‑148 (MQ=255)
             tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTCATAGTGGATTAGGTATTTCCGTCTTCCGCGAAAACTCCACCCGGACCCg                                                                                                                            >  1:3483662/1‑148 (MQ=255)
             tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                            <  1:675940/148‑1 (MQ=255)
             tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                            <  2:1094939/148‑1 (MQ=255)
             tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                            >  1:1814082/1‑148 (MQ=255)
               aCGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCt                                                                                                                          >  2:632219/1‑148 (MQ=255)
                cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                         <  2:3042238/148‑1 (MQ=255)
                cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                         <  1:521358/148‑1 (MQ=255)
                cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCtt                                                                                                                         <  1:554160/148‑1 (MQ=255)
                 gTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTtc                                                                                                                        <  2:2433903/148‑1 (MQ=255)
                  tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                       <  2:1925657/148‑1 (MQ=255)
                  tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                       <  1:2717892/148‑1 (MQ=255)
                   gggAACAGACGCCAGATGAACTGAATCAAGTACTGGACGTTGCAGCAGCATTAAAAGCACTGCGTGCTGATAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTACGCGAAAACTCCACCCGTACCCGCTtctc                                                                                                                      <  1:2480193/148‑1 (MQ=255)
                    ggAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTctct                                                                                                                     <  1:214775/148‑1 (MQ=255)
                     gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                    <  1:1547910/148‑1 (MQ=255)
                     gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                    <  2:2652436/148‑1 (MQ=255)
                     gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                    <  2:1534694/148‑1 (MQ=255)
                     gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                    >  2:1170876/1‑148 (MQ=255)
                      aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                   >  2:748622/1‑148 (MQ=255)
                      aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                   <  1:3594139/148‑1 (MQ=255)
                      aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                   <  2:1280955/148‑1 (MQ=255)
                       aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                  >  2:321680/1‑148 (MQ=255)
                       aCAGACGCCAGATGAACTGAAACAAGTACTGGACATTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                  <  2:644118/148‑1 (MQ=255)
                        cAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                 <  2:615418/148‑1 (MQ=255)
                         aGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGc                                                                                                                <  1:3589069/148‑1 (MQ=255)
                           aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                              <  1:1277563/148‑1 (MQ=255)
                           aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                              <  1:1788034/148‑1 (MQ=255)
                           aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                              <  1:1787821/148‑1 (MQ=255)
                            cGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTc                                                                                                             >  2:388317/1‑148 (MQ=255)
                                 gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                        <  1:441778/148‑1 (MQ=255)
                                      aCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAc                                                                                                   <  2:2381418/148‑1 (MQ=255)
                                       cTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAcc                                                                                                  >  1:2551712/1‑148 (MQ=255)
                                         gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                <  2:187963/148‑1 (MQ=255)
                                               aGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGc                                                                                          <  1:388317/148‑1 (MQ=255)
                                               aGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGc                                                                                          >  2:3767230/1‑148 (MQ=255)
                                                gTACTGGACGTTGCCGCATCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                         <  1:3326146/148‑1 (MQ=255)
                                                gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                         >  2:2599023/1‑148 (MQ=255)
                                                gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                         <  1:3541395/148‑1 (MQ=255)
                                                 tACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCt                                                                                        <  2:2796624/148‑1 (MQ=255)
                                                   cTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTg                                                                                      <  2:1009921/148‑1 (MQ=255)
                                                      gACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcac                                                                                   >  2:479495/1‑148 (MQ=255)
                                                       aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                  <  1:2430612/148‑1 (MQ=255)
                                                       aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                  <  1:2347487/148‑1 (MQ=255)
                                                               gcagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                          <  1:3293679/148‑1 (MQ=255)
                                                                        aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                 >  1:857789/1‑148 (MQ=255)
                                                                        aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                 <  2:2028297/148‑1 (MQ=255)
                                                                         aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                >  1:1014339/1‑148 (MQ=255)
                                                                          aGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGc                                                               >  1:2036876/1‑148 (MQ=255)
                                                                          aGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGc                                                               >  2:1381515/1‑148 (MQ=255)
                                                                            cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaa                                                             >  1:619996/1‑148 (MQ=255)
                                                                             aCTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaaa                                                            <  1:1477170/148‑1 (MQ=255)
                                                                              cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                           >  1:2790032/1‑148 (MQ=255)
                                                                                gCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATca                                                         <  1:1738513/148‑1 (MQ=255)
                                                                                   tGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaa                                                      >  2:1672857/1‑148 (MQ=255)
                                                                                    gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaaa                                                     >  2:2053561/1‑148 (MQ=255)
                                                                                     cTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAAt                                                    <  2:148778/148‑1 (MQ=255)
                                                                                        aaaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGc                                                 >  2:2012391/1‑148 (MQ=255)
                                                                                          aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTc                                               >  2:3604958/1‑148 (MQ=255)
                                                                                          aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTc                                               >  1:221449/1‑148 (MQ=255)
                                                                                          aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTc                                               <  2:778609/148‑1 (MQ=255)
                                                                                           aCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                              <  2:1538144/148‑1 (MQ=255)
                                                                                             aTCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg                                            >  1:155480/1‑148 (MQ=255)
                                                                                              tctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACgg                                           >  2:2206900/1‑148 (MQ=255)
                                                                                               ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGc                                          >  2:315561/1‑148 (MQ=255)
                                                                                               ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGc                                          <  1:3070209/148‑1 (MQ=255)
                                                                                                tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCg                                         <  2:1300326/148‑1 (MQ=255)
                                                                                                tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCg                                         <  2:3681292/148‑1 (MQ=255)
                                                                                                       aaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgt                                  >  2:1103252/1‑148 (MQ=255)
                                                                                                        aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtg                                 <  2:1802983/148‑1 (MQ=255)
                                                                                                            tttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgtg                             >  1:1376938/1‑148 (MQ=255)
                                                                                                              tAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaa                           <  2:202013/148‑1 (MQ=255)
                                                                                                               aaTAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaaa                          <  2:3721989/148‑1 (MQ=255)
                                                                                                                    tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                     >  1:2903680/1‑148 (MQ=255)
                                                                                                                      gATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                   >  1:1390819/1‑148 (MQ=255)
                                                                                                                        ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAata                 <  2:1373753/148‑1 (MQ=255)
                                                                                                                        ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAata                 <  1:1257886/148‑1 (MQ=255)
                                                                                                                          aGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATg               <  1:1170876/148‑1 (MQ=255)
                                                                                                                           ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCCATATGa              >  2:2719747/1‑148 (MQ=255)
                                                                                                                            gTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAGGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAt             >  1:835284/1‑148 (MQ=255)
                                                                                                                               tttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc          <  2:1986818/148‑1 (MQ=255)
                                                                                                                               tttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc          <  1:882497/148‑1 (MQ=255)
                                                                                                                               tttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc          <  1:1923513/148‑1 (MQ=255)
                                                                                                                                ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCNCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc         >  2:16032/1‑148 (MQ=255)
                                                                                                                                ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc         >  1:1692449/1‑148 (MQ=255)
                                                                                                                                ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc         <  1:479495/148‑1 (MQ=255)
                                                                                                                                ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc         >  1:3417516/1‑148 (MQ=255)
                                                                                                                                ttCCGTAGTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc         <  2:1019909/148‑1 (MQ=255)
                                                                                                                                 tCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCt        >  1:719718/1‑148 (MQ=255)
                                                                                                                                  ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt       <  2:1698781/148‑1 (MQ=255)
                                                                                                                                  ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt       <  1:562701/148‑1 (MQ=255)
                                                                                                                                   cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc      >  1:1507354/1‑148 (MQ=255)
                                                                                                                                   cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc      <  2:1737432/148‑1 (MQ=255)
                                                                                                                                   cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc      <  2:2110305/148‑1 (MQ=255)
                                                                                                                                       ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg  >  1:3235026/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                           
AAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCG  >  CP009273/2999684‑2999966

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: