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,8270CA96.1% 552.6 / 12.1 181D69E (GAC→GAAygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (78/96);  total (80/101)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 182 total aligned reads displayed.

AGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCG  >  CP009273/2999682‑2999969
                                                                                                                                                 |                                                                                                                                              
agaAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGaaaa                                                                                                                                              <  1:1377824/148‑1 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                            >  1:959540/1‑148 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                            <  1:2186667/148‑1 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                            >  1:1678571/1‑148 (MQ=255)
   aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                           <  1:666654/148‑1 (MQ=255)
   aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                           <  1:666653/148‑1 (MQ=255)
     gATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCa                                                                                                                                         >  1:1721110/1‑148 (MQ=255)
          tttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGt                                                                                                                                    <  2:2650323/148‑1 (MQ=255)
            tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                  <  2:1677522/148‑1 (MQ=255)
            tGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTAc                                                                                                                                  <  1:1861437/148‑1 (MQ=255)
                 aCGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCt                                                                                                                             <  1:179646/148‑1 (MQ=255)
                 aCGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCt                                                                                                                             <  1:179618/148‑1 (MQ=255)
                    tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                          <  2:3445501/148‑1 (MQ=255)
                       gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                       <  1:1740983/148‑1 (MQ=255)
                       gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                       <  1:1195011/148‑1 (MQ=255)
                         aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                     <  2:1890245/148‑1 (MQ=255)
                          cAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                    <  1:2484632/148‑1 (MQ=255)
                             aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                                 >  1:1275998/1‑148 (MQ=255)
                             aCGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCtt                                                                                                                 >  2:1870504/1‑148 (MQ=255)
                              cGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTc                                                                                                                <  1:40740/148‑1 (MQ=255)
                                   gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                           <  2:1722233/148‑1 (MQ=255)
                                   gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                           <  2:288970/148‑1 (MQ=255)
                                    aTGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCt                                                                                                          <  2:1131309/148‑1 (MQ=255)
                                    aTGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCt                                                                                                          <  1:87977/148‑1 (MQ=255)
                                      gAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTa                                                                                                        >  2:2623146/1‑148 (MQ=255)
                                       aaCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTaa                                                                                                       >  2:729775/1‑148 (MQ=255)
                                       aaCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTaa                                                                                                       <  1:3050526/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                   <  2:3180640/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                   <  2:1220322/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                   <  2:2871500/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                   <  1:3422065/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                   <  2:2869661/148‑1 (MQ=255)
                                             aaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt                                                                                                 <  2:2224836/148‑1 (MQ=255)
                                                  gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                            <  2:3310263/148‑1 (MQ=255)
                                                  gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                            >  2:3448400/1‑148 (MQ=255)
                                                      tGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGc                                                                                        <  1:2881160/148‑1 (MQ=255)
                                                         aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                     <  1:1521945/148‑1 (MQ=255)
                                                            ttGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaaca                                                                                  <  1:1550990/148‑1 (MQ=255)
                                                             tGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCAcaacaa                                                                                 >  2:623487/1‑148 (MQ=255)
                                                              gCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                                                                <  2:2519210/148‑1 (MQ=255)
                                                                cGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAt                                                                              <  1:2238345/148‑1 (MQ=255)
                                                                 gcatcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGTTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                             <  1:2837737/148‑1 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            >  1:1941140/1‑148 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            <  2:2551065/148‑1 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtct                                                                            <  2:3181183/148‑1 (MQ=255)
                                                                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTGAACCAGCTCGGCCTTGCACACCACGAtct                                                                            >  1:642987/1‑148 (MQ=255)
                                                                     caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                                                                         <  1:3361871/148‑1 (MQ=255)
                                                                     caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                                                                         >  1:2309767/1‑148 (MQ=255)
                                                                        tAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCGTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                                                                      >  2:2948061/1‑148 (MQ=255)
                                                                         aaaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                     <  1:606830/148‑1 (MQ=255)
                                                                         aaaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                     >  2:2044201/1‑148 (MQ=255)
                                                                         aaaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                     <  1:2126766/148‑1 (MQ=255)
                                                                          aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                    <  2:1557371/148‑1 (MQ=255)
                                                                              cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaa                                                                <  2:2452351/148‑1 (MQ=255)
                                                                              cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaa                                                                <  2:3358295/148‑1 (MQ=255)
                                                                                 tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATc                                                             >  2:607279/1‑148 (MQ=255)
                                                                                      gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaaa                                                        >  1:3054551/1‑148 (MQ=255)
                                                                                         gAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                     >  1:1649556/1‑148 (MQ=255)
                                                                                                 ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGc                                             >  2:1024414/1‑148 (MQ=255)
                                                                                                 ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGc                                             >  1:835425/1‑148 (MQ=255)
                                                                                                   cAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGa                                           >  1:1638163/1‑148 (MQ=255)
                                                                                                    aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGaa                                          <  1:1001201/148‑1 (MQ=255)
                                                                                                      acAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAc                                        <  2:2364001/147‑1 (MQ=255)
                                                                                                       cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAAcc                                       <  2:427516/148‑1 (MQ=255)
                                                                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtg                                    >  2:1959023/1‑148 (MQ=255)
                                                                                                            tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaa                                                         >  2:1259216/1‑125 (MQ=255)
                                                                                                             cTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                 <  1:1939132/148‑1 (MQ=255)
                                                                                                                tAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaa                              <  2:1720868/148‑1 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                        >  2:1759317/1‑148 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                        >  1:130295/1‑148 (MQ=255)
                                                                                                                       ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                       <  2:3098864/148‑1 (MQ=255)
                                                                                                                        gATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                      >  1:1218364/1‑148 (MQ=255)
                                                                                                                        gATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                      >  1:2378544/1‑148 (MQ=255)
                                                                                                                        gATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCaa                      >  2:3244185/1‑148 (MQ=255)
                                                                                                                         aTTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAat                     <  1:3549342/148‑1 (MQ=255)
                                                                                                                          ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTAGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAata                    <  1:2008061/148‑1 (MQ=255)
                                                                                                                          ttAGGTATTTCCGTATTCCGCGAAAACTCCACCCGGACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTAGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAata                    <  1:2008057/148‑1 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                   >  2:520238/1‑148 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                   <  1:97278/148‑1 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                   >  2:3093693/1‑148 (MQ=255)
                                                                                                                              gTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAt                <  1:326864/148‑1 (MQ=255)
                                                                                                                               tATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtc               >  2:3126334/1‑148 (MQ=255)
                                                                                                                               tATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtc               >  2:1344142/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            >  1:1861314/1‑148 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaa                              >  2:1437768/1‑130 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaa                              <  1:1437768/130‑1 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            <  2:882393/148‑1 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc            <  2:1937266/148‑1 (MQ=255)
                                                                                                                                   tCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCt           <  1:752447/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATAACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt          <  1:1011265/148‑1 (MQ=255)
                                                                                                                                     cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc         <  2:959540/148‑1 (MQ=255)
                                                                                                                                       tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg       >  1:3500199/1‑148 (MQ=255)
                                                                                                                                       tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg       >  1:3496069/1‑148 (MQ=255)
                                                                                                                                         ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg     >  2:1533052/1‑148 (MQ=255)
                                                                                                                                         ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg     >  2:677824/1‑148 (MQ=255)
                                                                                                                                          tCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcgc    >  1:1321914/1‑148 (MQ=255)
                                                                                                                                           ccGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCGCAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGcc   >  1:963834/1‑148 (MQ=255)
                                                                                                                                           ccGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCCATATGATCTCCTTCTGCGcc   >  1:2698577/1‑148 (MQ=255)
                                                                                                                                            cgcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCTCCAATATGATCTCCTTCTGCGCCg  >  1:793628/1‑148 (MQ=255)
                                                                                                                                            cgcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCg  <  2:3521284/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: