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,8270CA98.5% 932.2 / 2.8 288D69E (GAC→GAAygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (152/130);  total (156/132)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 292 total aligned reads displayed.

CGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACG  >  CP009273/2999680‑2999972
                                                                                                                                                   |                                                                                                                                                 
  agaAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACt                                                                                                                                               >  1:3500757/1‑150 (MQ=255)
      aGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCACAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                           <  1:4079201/150‑1 (MQ=255)
        aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAccc                                                                                                                                         >  1:1401095/1‑150 (MQ=255)
            tttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                     <  2:4987514/150‑1 (MQ=255)
              tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAccc                                                                                                                                   <  2:2342690/150‑1 (MQ=255)
              tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAccc                                                                                                                                   <  2:695966/150‑1 (MQ=255)
               gTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                                  >  1:1872501/1‑150 (MQ=255)
               gTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                                  >  2:1991330/1‑150 (MQ=255)
               gTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                                  >  1:5197396/1‑150 (MQ=255)
               gTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                                  >  2:415337/1‑150 (MQ=255)
                    cGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                             <  1:3471763/150‑1 (MQ=255)
                     gTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtctc                                                                                                                            >  2:5166049/1‑150 (MQ=255)
                      tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTctct                                                                                                                           >  2:2819714/1‑150 (MQ=255)
                       gggAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                          <  2:1521946/150‑1 (MQ=255)
                        ggAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTa                                                                                                                         >  1:444572/1‑150 (MQ=255)
                         gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                        <  2:5025282/150‑1 (MQ=255)
                         gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCATATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                        <  2:4127350/150‑1 (MQ=255)
                          aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                       <  2:1656365/150‑1 (MQ=255)
                          aaCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                       >  1:2455811/1‑150 (MQ=255)
                           aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGc                                                                                                                      >  1:840843/1‑150 (MQ=255)
                             aGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                                    <  2:3012073/150‑1 (MQ=255)
                               aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACTTCTCAACCAAAGTCTTTATTAGTGGTTTAGGTATTTCCGTATTCCGCGTTATCTCATACCGTACCCGATTCTCTTATGCTTcc                                                                                                                  <  1:822280/150‑1 (MQ=255)
                               aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTcc                                                                                                                  <  1:578225/150‑1 (MQ=255)
                                cGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCcg                                                                                                                 <  1:1117102/150‑1 (MQ=255)
                                    aGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCt                                                                                                             >  1:4756606/1‑150 (MQ=255)
                                     gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCtt                                                                                                            <  2:3129070/150‑1 (MQ=255)
                                             gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt                                                                                                    <  1:3332837/150‑1 (MQ=255)
                                             gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt                                                                                                    <  2:578504/150‑1 (MQ=255)
                                              aaaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                                                                   <  2:1794587/150‑1 (MQ=255)
                                              aaaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                                                                   <  1:2854687/150‑1 (MQ=255)
                                              aaaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                                                                   <  2:3453259/150‑1 (MQ=255)
                                              aaaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                                                                   >  1:4115202/1‑150 (MQ=255)
                                               aaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCg                                                                                                  <  2:922278/150‑1 (MQ=255)
                                                  aaGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                               >  1:1767605/1‑150 (MQ=255)
                                                  aaGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                               >  2:4362810/1‑150 (MQ=255)
                                                     tACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTg                                                                                            <  2:2551516/150‑1 (MQ=255)
                                                      aCTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGc                                                                                           >  1:2221140/1‑150 (MQ=255)
                                                           aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCACAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaac                                                                                      <  2:3769024/150‑1 (MQ=255)
                                                                 ccGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                                <  1:3640955/150‑1 (MQ=255)
                                                                  cGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                               <  2:1377940/150‑1 (MQ=255)
                                                                  cGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                               >  1:2466237/1‑150 (MQ=255)
                                                                    cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcg                                                                             <  1:3177468/150‑1 (MQ=255)
                                                                    cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcg                                                                             <  1:499854/150‑1 (MQ=255)
                                                                       caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacg                                                                          >  1:4232329/1‑150 (MQ=255)
                                                                        aTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                                                                         >  2:1818993/1‑150 (MQ=255)
                                                                          tAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                       >  2:1478186/1‑150 (MQ=255)
                                                                           aaaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                      >  1:2783723/1‑150 (MQ=255)
                                                                            aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGc                                                                     <  1:3677770/150‑1 (MQ=255)
                                                                                cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                                 >  2:1313863/1‑150 (MQ=255)
                                                                                        gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATc                                                         <  2:1038987/150‑1 (MQ=255)
                                                                                          tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAATATCTCGACGAAGGCAAATCACAAATCGc                                                       >  1:2572533/1‑150 (MQ=255)
                                                                                           gAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCATGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCt                                                      >  1:633777/1‑150 (MQ=255)
                                                                                            aaaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACACATCGCTc                                                     >  2:300170/1‑150 (MQ=255)
                                                                                            aaaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGAATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTCTGCTTCCGCGCTACACCTGCTCGGCCTTGCACACCCAGAACTCGACGAAGGCCCATCACACAGCGCTc                                                     >  1:4547650/1‑150 (MQ=255)
                                                                                             aaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                                    >  2:372855/1‑150 (MQ=255)
                                                                                             aaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                                    >  1:15149/1‑150 (MQ=255)
                                                                                              anCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACACCACGATCTCGACGACGGCACATCACANANCGCTccc                                                   >  1:1862763/1‑148 (MQ=255)
                                                                                               aCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg                                                  <  1:1096976/150‑1 (MQ=255)
                                                                                                      aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAc                                           >  2:3799708/1‑150 (MQ=255)
                                                                                                      aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAc                                           <  2:3306594/150‑1 (MQ=255)
                                                                                                       aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                          <  2:4466943/150‑1 (MQ=255)
                                                                                                       aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                          <  2:4201323/150‑1 (MQ=255)
                                                                                                       aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                          <  1:5104849/150‑1 (MQ=255)
                                                                                                       aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                          <  2:5324984/150‑1 (MQ=255)
                                                                                                         cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgt                                        >  1:3821051/1‑150 (MQ=255)
                                                                                                          aaaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtg                                       >  2:4603964/1‑150 (MQ=255)
                                                                                                           aaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgc                                      <  1:933475/150‑1 (MQ=255)
                                                                                                             ggcTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCAGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACTAAGGCAAATCACAAATCGCTCACGTCGAAACagtgcgt                                    >  1:740834/3‑150 (MQ=255)
                                                                                                             gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                    <  2:3161220/150‑1 (MQ=255)
                                                                                                             gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                    <  2:3349963/150‑1 (MQ=255)
                                                                                                                tttAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTTCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaa                                 >  1:3538479/1‑150 (MQ=255)
                                                                                                                  tAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAAc                               <  1:1197305/150‑1 (MQ=255)
                                                                                                                    aTAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCg                             <  1:2154218/150‑1 (MQ=255)
                                                                                                                    aTAGTGGAGTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAAACg                             >  1:557003/1‑150 (MQ=255)
                                                                                                                        tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                         >  1:1646931/1‑150 (MQ=255)
                                                                                                                        tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                         <  1:894724/150‑1 (MQ=255)
                                                                                                                         ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAat                        >  1:3741307/1‑150 (MQ=255)
                                                                                                                           aTTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                      <  1:3128834/150‑1 (MQ=255)
                                                                                                                           aTTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                      <  2:4457027/150‑1 (MQ=255)
                                                                                                                            ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATg                     <  1:3364384/150‑1 (MQ=255)
                                                                                                                              aGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAAAATGAt                   >  2:2124541/1‑150 (MQ=255)
                                                                                                                              aGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAt                   <  2:2504642/150‑1 (MQ=255)
                                                                                                                               ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtc                  <  2:888420/150‑1 (MQ=255)
                                                                                                                               ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtc                  >  2:1436696/1‑150 (MQ=255)
                                                                                                                                 tATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc                >  2:2067676/1‑150 (MQ=255)
                                                                                                                                    ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt             <  2:3999643/150‑1 (MQ=255)
                                                                                                                                     tCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc            >  2:4743833/1‑150 (MQ=255)
                                                                                                                                     tCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc            <  2:79571/150‑1 (MQ=255)
                                                                                                                                      ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTTCGTGAAACCGCCAATATGATCTCCTTCt           >  2:2745511/1‑150 (MQ=255)
                                                                                                                                      ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCt           >  2:2968732/1‑150 (MQ=255)
                                                                                                                                       cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg          >  2:4863653/1‑150 (MQ=255)
                                                                                                                                       cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTAGCACAACAAGATCTCGACGAAGGCAAATCACAAAACGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg          >  1:65074/1‑150 (MQ=255)
                                                                                                                                        gTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgc         >  2:5018819/1‑150 (MQ=255)
                                                                                                                                        gTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgc         >  1:1317710/1‑150 (MQ=255)
                                                                                                                                         tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg        >  2:567686/1‑150 (MQ=255)
                                                                                                                                         tATTCCGCGAAAACTCCACCCGTACCCGCGTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg        >  1:395242/1‑150 (MQ=255)
                                                                                                                                           ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGcc      <  2:4534437/150‑1 (MQ=255)
                                                                                                                                           ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGcc      >  1:3781857/1‑150 (MQ=255)
                                                                                                                                           ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAAGATGATCTCCTTCTGCGcc      >  2:4702716/1‑150 (MQ=255)
                                                                                                                                               gcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACg  >  1:4697830/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACG  >  CP009273/2999680‑2999972

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: