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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,999,8270CA89.9% 511.9 / 53.5 199D69E (GAC→GAAygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (6/14);  new base A (86/93);  total (92/107)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 203 total aligned reads displayed.

AGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGA  >  CP009273/2999682‑2999970
                                                                                                                                                 |                                                                                                                                               
agaAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGaaaa                                                                                                                                               >  2:2723325/1‑148 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                             >  1:1002323/1‑148 (MQ=255)
  aaaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACt                                                                                                                                             >  1:1443407/1‑148 (MQ=255)
   aaGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTc                                                                                                                                            <  1:1942290/148‑1 (MQ=255)
    aGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTcc                                                                                                                                           >  1:3276422/1‑148 (MQ=255)
      aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                         <  1:916720/148‑1 (MQ=255)
      aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                         <  2:1361483/148‑1 (MQ=255)
      aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                         >  2:1927613/1‑148 (MQ=255)
      aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                         >  1:2175157/1‑148 (MQ=255)
      aTTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTACTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAc                                                                                                                                         <  2:1361641/148‑1 (MQ=255)
       ttttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAcc                                                                                                                                        <  1:1271507/148‑1 (MQ=255)
       ttttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAcc                                                                                                                                        <  2:2036148/148‑1 (MQ=255)
        tttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAACACTCCAccc                                                                                                                                       >  1:2633729/1‑148 (MQ=255)
        tttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCAccc                                                                                                                                       >  2:2848074/1‑148 (MQ=255)
         ttttGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCg                                                                                                                                      >  2:2159005/1‑148 (MQ=255)
               tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCg                                                                                                                                >  2:1250070/1‑148 (MQ=255)
               tAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCg                                                                                                                                >  1:690932/1‑148 (MQ=255)
                 aCGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCt                                                                                                                              <  1:3022977/148‑1 (MQ=255)
                    tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                           <  2:3021843/148‑1 (MQ=255)
                    tGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTtct                                                                                                                           >  1:2827162/1‑148 (MQ=255)
                       gAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCtt                                                                                                                        >  2:1816365/1‑148 (MQ=255)
                         aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                      <  1:1855130/148‑1 (MQ=255)
                         aCAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTAt                                                                                                                      >  1:750203/1‑148 (MQ=255)
                          cAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                     >  2:487299/1‑148 (MQ=255)
                          cAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATg                                                                                                                     <  1:607365/148‑1 (MQ=255)
                            gACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCt                                                                                                                   >  1:367381/1‑148 (MQ=255)
                                   gATGAACTGACACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAACGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                            <  2:54943/148‑1 (MQ=255)
                                   gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                            >  1:732832/1‑148 (MQ=255)
                                   gATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCgcgc                                                                                                            >  1:2392064/1‑148 (MQ=255)
                                    aTGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCt                                                                                                           >  1:2261323/1‑148 (MQ=255)
                                      gAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTa                                                                                                         <  1:1391506/148‑1 (MQ=255)
                                      gAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTa                                                                                                         >  1:1161468/1‑148 (MQ=255)
                                       aaCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTaa                                                                                                        <  2:1509203/148‑1 (MQ=255)
                                        aCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAc                                                                                                       <  1:2203596/148‑1 (MQ=255)
                                        aCTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAc                                                                                                       <  1:1780114/148‑1 (MQ=255)
                                         cTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAAcc                                                                                                      <  1:1573824/148‑1 (MQ=255)
                                           gAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                                                                    <  1:3128991/148‑1 (MQ=255)
                                            aaaCAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGc                                                                                                   <  2:1479452/148‑1 (MQ=255)
                                                 aGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGGACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGc                                                                                              <  2:907832/148‑1 (MQ=255)
                                                  gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                             <  2:1405696/148‑1 (MQ=255)
                                                  gTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGcc                                                                                             >  1:1647316/1‑148 (MQ=255)
                                                   tACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCt                                                                                            >  1:92764/1‑148 (MQ=255)
                                                         aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                      >  2:288132/1‑148 (MQ=255)
                                                         aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                      >  2:288047/1‑148 (MQ=255)
                                                          cGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaa                                                                                     >  1:3223802/1‑148 (MQ=255)
                                                          cGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaa                                                                                     >  1:2559149/1‑148 (MQ=255)
                                                            ttGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaaca                                                                                   >  2:792727/1‑148 (MQ=255)
                                                              gCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                                                                 <  2:3202793/148‑1 (MQ=255)
                                                              gCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                                                                 <  1:521367/148‑1 (MQ=255)
                                                              gCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                                                                 <  1:2941645/148‑1 (MQ=255)
                                                                cGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAt                                                                               <  1:2313194/148‑1 (MQ=255)
                                                                     caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                                                                          <  2:2374667/148‑1 (MQ=255)
                                                                     caTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                                                                          <  2:712347/148‑1 (MQ=255)
                                                                        tAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                                                                       >  2:255596/1‑148 (MQ=255)
                                                                         aaaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                      <  2:2869890/148‑1 (MQ=255)
                                                                          aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                     <  1:3070538/148‑1 (MQ=255)
                                                                          aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                                                                     <  2:1002323/148‑1 (MQ=255)
                                                                              cACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaa                                                                 >  1:2362042/1‑148 (MQ=255)
                                                                                cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                               >  1:102832/1‑148 (MQ=255)
                                                                                cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAt                                                               >  1:301102/1‑148 (MQ=255)
                                                                                 tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATc                                                              <  1:158155/148‑1 (MQ=255)
                                                                                  gCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATca                                                             <  1:585801/148‑1 (MQ=255)
                                                                                      gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCACATCACaaa                                                         >  2:399438/1‑148 (MQ=255)
                                                                                      gCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaaa                                                         >  2:399278/1‑148 (MQ=255)
                                                                                                 ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGc                                              >  1:381007/1‑148 (MQ=255)
                                                                                                 ctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGc                                              >  1:2073045/1‑148 (MQ=255)
                                                                                                      ccAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAAc                                         >  2:2943397/1‑148 (MQ=255)
                                                                                                        aaaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcg                                       >  1:2487383/1‑148 (MQ=255)
                                                                                                                 aaTAGTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGaaa                              <  1:1212921/148‑1 (MQ=255)
                                                                                                                     gTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGc                          <  1:949757/148‑1 (MQ=255)
                                                                                                                     gTGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGc                          <  2:661964/148‑1 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                         >  1:1844792/1‑148 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGcc                         >  1:1338412/1‑148 (MQ=255)
                                                                                                                      tGGATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGcc                         >  1:985719/1‑148 (MQ=255)
                                                                                                                       ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAACTCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                        <  2:469847/148‑1 (MQ=255)
                                                                                                                       ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                        <  2:450422/148‑1 (MQ=255)
                                                                                                                       ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                        <  2:1680665/148‑1 (MQ=255)
                                                                                                                       ggATTAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                        <  2:2673326/148‑1 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                    >  1:2027349/1‑148 (MQ=255)
                                                                                                                           tAGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAatat                    >  1:2027290/1‑148 (MQ=255)
                                                                                                                            aGGTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATg                   <  2:259307/148‑1 (MQ=255)
                                                                                                                             ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGa                  >  2:2217965/1‑148 (MQ=255)
                                                                                                                             ggTATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGa                  <  1:2150690/148‑1 (MQ=255)
                                                                                                                               tATTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCTAAACCGTGAGTGAAACAGCCAATATGAtc                >  2:1409651/1‑148 (MQ=255)
                                                                                                                                aTTTCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtct               >  1:2754953/1‑148 (MQ=255)
                                                                                                                                 tttCCGTATTCCTCGAATACTCCACCCGTACCCTCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGTCGAATCCGTGCGTGAAACCGCCAATATGAtctc              <  2:2231162/148‑1 (MQ=255)
                                                                                                                                 tttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGAtctc              <  2:1258112/148‑1 (MQ=255)
                                                                                                                                  ttCCGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTcc             <  2:3065908/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt           <  1:3153504/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt           <  2:530675/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt           <  1:3163653/148‑1 (MQ=255)
                                                                                                                                    ccGTATTCCGCGAAAACTCCACCCGTACCCGCGTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt           <  2:1634112/148‑1 (MQ=255)
                                                                                                                                     cGTATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc          <  1:1872736/148‑1 (MQ=255)
                                                                                                                                     cGTAGTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTc          <  2:1920887/148‑1 (MQ=255)
                                                                                                                                       tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg        >  2:412936/1‑148 (MQ=255)
                                                                                                                                       tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg        >  2:1613490/1‑148 (MQ=255)
                                                                                                                                       tATTCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTg        >  1:3077370/1‑148 (MQ=255)
                                                                                                                                         ttCCGCGAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgcg      >  2:3298278/1‑148 (MQ=255)
                                                                                                                                             gcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGa  <  2:3223802/148‑1 (MQ=255)
                                                                                                                                             gcgAAAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATAGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGa  <  1:1938232/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGA  >  CP009273/2999682‑2999970

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: