Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,888,274 C→A 30.1% P88Q (CCG→CAG)  chrR → chromate reductase, Class I, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,888,2740CA30.1% 408.3 / 138.7 216P88Q (CCG→CAG) chrRchromate reductase, Class I, flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (84/65);  new base A (42/24);  total (126/90)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

CGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCC  >  CP009273/3888127‑3888420
                                                                                                                                                   |                                                                                                                                                  
cGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                                                                                                                    <  1:200615/148‑1 (MQ=255)
cGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                                                                                                                    <  1:197747/148‑1 (MQ=255)
  ttACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCgg                                                                                                                                                  <  2:3211184/148‑1 (MQ=255)
    aCCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGt                                                                                                                                                >  2:1510325/1‑148 (MQ=255)
      catccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTgg                                                                                                                                              >  1:2326096/1‑148 (MQ=255)
       atccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                                                                                                                             <  2:718944/148‑1 (MQ=255)
        tccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGc                                                                                                                                            >  2:2730581/1‑148 (MQ=255)
               ccGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                                                                                                                     <  1:2626923/148‑1 (MQ=255)
               ccGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                                                                                                                     >  1:3168722/1‑148 (MQ=255)
               ccGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                                                                                                                     <  2:1294162/148‑1 (MQ=255)
               ccGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                                                                                                                     <  2:529154/148‑1 (MQ=255)
                 gACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                                                                                                                   >  1:782896/1‑148 (MQ=255)
                  aCATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGTTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGc                                                                                                                                  <  1:352862/148‑1 (MQ=255)
                  aCATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGc                                                                                                                                  <  1:2016286/148‑1 (MQ=255)
                    aTTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                                                                                                                <  1:3065370/148‑1 (MQ=255)
                     ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGACTATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCAt                                                                                                                               >  2:1012329/1‑148 (MQ=255)
                     ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCAt                                                                                                                               >  2:2788817/1‑148 (MQ=255)
                      tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                                                                                                                              >  1:2712311/1‑148 (MQ=255)
                      tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATc                                                                                                                              <  2:2054235/148‑1 (MQ=255)
                       cccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCg                                                                                                                             <  2:1944831/148‑1 (MQ=255)
                         cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGAc                                                                                                                           <  2:1160156/148‑1 (MQ=255)
                           tGTATGACGCTGACGTACAGCAGGAAGATGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACTGATCCGTCTGGCTGAGGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                                                                                                                         <  1:1023783/148‑1 (MQ=255)
                           tGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                                                                                                                         >  2:1372804/1‑148 (MQ=255)
                                gACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGACGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCttt                                                                                                                    <  2:2435382/148‑1 (MQ=255)
                                 aCGCTGACGTACAGCAGTAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                                                                                                                   <  2:1558396/148‑1 (MQ=255)
                                 aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                                                                                                                   <  2:2727445/148‑1 (MQ=255)
                                 aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCAGGCGGAACAGATCCGTCAGGCGGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                                                                                                                   <  2:1652570/148‑1 (MQ=255)
                                  cGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                                                                                                                  >  2:3507763/1‑148 (MQ=255)
                                   gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTccc                                                                                                                 >  2:2853151/1‑148 (MQ=255)
                                    cTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCg                                                                                                                >  1:1862353/1‑148 (MQ=255)
                                        cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTg                                                                                                            >  1:871711/1‑148 (MQ=255)
                                              gcagTAAGATGGTTTTCCAGCAACGTTTGAAGCTCTGGCGTAAAAGACCCGTCCGGCTGACGGTGTGGTGATCGCCACGCCGTAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGAt                                                                                                      <  1:3318266/148‑1 (MQ=255)
                                               cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc                                                                                                     >  1:2748338/1‑148 (MQ=255)
                                               cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc                                                                                                     >  2:1853235/1‑148 (MQ=255)
                                                  gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc                                                                                                  >  1:1365331/1‑148 (MQ=255)
                                                    agaagGTTTTCCAGCAACGGTTGAAGCTCGGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCg                                                                                                <  1:2586177/148‑1 (MQ=255)
                                                      aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCt                                                                                              >  1:723629/1‑148 (MQ=255)
                                                       agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTg                                                                                             >  1:891898/1‑148 (MQ=255)
                                                           tttCCACCAACGGTTGACGCTCGGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGACACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCg                                                                                         <  2:920003/148‑1 (MQ=255)
                                                             tCCAGCAACGGTTGAAGCTCAGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGACACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGt                                                                                       <  1:501673/148‑1 (MQ=255)
                                                               cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaa                                                                                     >  1:2858577/1‑148 (MQ=255)
                                                               cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaa                                                                                     >  1:6487/1‑148 (MQ=255)
                                                                accacgggTTGAAGCTCTGGGGTAACATTTCCGTCAGGCTGTGGGTGTGTTGATCGTCAGGCCGTATTTGAACTACTCGTTACCGGGTGGGCTGAAAATTGCACTCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaaa                                                                                    <  1:1012329/142‑1 (MQ=255)
                                                                aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaaa                                                                                    <  2:3300748/148‑1 (MQ=255)
                                                                aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaaa                                                                                    >  2:561733/1‑148 (MQ=255)
                                                                  cAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGACAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAAcc                                                                                  >  2:3188821/1‑148 (MQ=255)
                                                                   aaCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCg                                                                                 <  2:2020181/148‑1 (MQ=255)
                                                                    aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                                >  2:255143/1‑148 (MQ=255)
                                                                    aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                                >  2:437735/1‑148 (MQ=255)
                                                                     cGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGt                                                                               >  1:1731452/1‑148 (MQ=255)
                                                                       gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTAt                                                                             >  1:1223370/1‑148 (MQ=255)
                                                                        ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTACACCGGTAtt                                                                            >  2:794878/1‑148 (MQ=255)
                                                                          gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGa                                                                          <  2:64354/148‑1 (MQ=255)
                                                                          gAAGCTCTGGCGGAACAGATCCGTCAGGCGGACGGTGTGGTGATCGTCACGCCGTAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCGGGCCGGTAAACCGGTATTGa                                                                          <  2:3352337/148‑1 (MQ=255)
                                                                            aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAtt                                                                        <  1:2034723/148‑1 (MQ=255)
                                                                             gCTCGGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTc                                                                       <  2:1045307/148‑1 (MQ=255)
                                                                              ctctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCa                                                                      >  2:1845918/1‑148 (MQ=255)
                                                                               tctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAg                                                                     <  1:2853151/148‑1 (MQ=255)
                                                                               tctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAg                                                                     >  1:2849298/1‑148 (MQ=255)
                                                                                 tGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGAc                                                                   >  2:3172841/1‑148 (MQ=255)
                                                                                  ggCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGAcc                                                                  <  2:1084443/148‑1 (MQ=255)
                                                                                       aaCAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTc                                                             >  1:1648651/1‑148 (MQ=255)
                                                                                       aaCAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTc                                                             >  2:83261/1‑148 (MQ=255)
                                                                                        aCAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCa                                                            >  2:3184349/1‑148 (MQ=255)
                                                                                         cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                           >  2:2799080/1‑148 (MQ=255)
                                                                                            aTCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATgg                                                        <  2:1621623/148‑1 (MQ=255)
                                                                                             tCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATggg                                                       <  1:1467829/148‑1 (MQ=255)
                                                                                               ggTCTGGCTGAGGGTGTGGTGATCGTCACGCCGTAATATAACTACTCGGTACGGGGTGGGCTTAAAAATGCCATNGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCg                                                     <  1:1806564/147‑1 (MQ=255)
                                                                                                gTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGt                                                    <  2:3042291/148‑1 (MQ=255)
                                                                                                 tCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTg                                                   >  2:1890421/1‑148 (MQ=255)
                                                                                                 tCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTg                                                   >  1:1887788/1‑148 (MQ=255)
                                                                                                 tCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTg                                                   >  1:481456/1‑148 (MQ=255)
                                                                                                 tCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTg                                                   >  1:368925/1‑148 (MQ=255)
                                                                                                  cAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGa                                                  >  1:1409635/1‑148 (MQ=255)
                                                                                                     gCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTg                                               >  1:2885627/1‑148 (MQ=255)
                                                                                                     gCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTg                                               >  1:2965144/1‑148 (MQ=255)
                                                                                                       tGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTggc                                             <  1:2730581/148‑1 (MQ=255)
                                                                                                          cGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTggcggc                                          >  2:1785996/1‑148 (MQ=255)
                                                                                                           ggTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGgcggcg                                         >  1:2404317/1‑148 (MQ=255)
                                                                                                            gtgtGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCTGACCAGCTCAATGGGCGTGATTGGCGgcgc                                        <  2:1952940/148‑1 (MQ=255)
                                                                                                               tggtgATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCTGACCAGCTCAATGGGCGTGATTGGCGgcgcgcg                                     <  2:3152920/148‑1 (MQ=255)
                                                                                                               tggtgATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcgcg                                     <  2:670857/148‑1 (MQ=255)
                                                                                                               tggtgATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcgcg                                     <  1:367970/148‑1 (MQ=255)
                                                                                                                ggtgATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcgcgc                                    >  2:500645/1‑148 (MQ=255)
                                                                                                                  tgATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTg                                  >  1:3199708/1‑148 (MQ=255)
                                                                                                                  tgATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTg                                  >  2:765098/1‑148 (MQ=255)
                                                                                                                   gATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGt                                 >  1:2254751/1‑148 (MQ=255)
                                                                                                                     tCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCa                               >  1:1180895/1‑148 (MQ=255)
                                                                                                                       gTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGt                             <  2:2157483/148‑1 (MQ=255)
                                                                                                                             ccGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCAcc                       >  1:972157/1‑148 (MQ=255)
                                                                                                                                  atatAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgcgc                  <  1:1925634/148‑1 (MQ=255)
                                                                                                                                     taACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAg               >  1:1832140/1‑148 (MQ=255)
                                                                                                                                      aaCTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGa              >  2:4674/1‑148 (MQ=255)
                                                                                                                                        ctactCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAtt            <  1:847080/148‑1 (MQ=255)
                                                                                                                                          actCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCTGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCt          >  1:3429257/1‑148 (MQ=255)
                                                                                                                                             cGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGt       >  1:2181661/1‑148 (MQ=255)
                                                                                                                                                tACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGtttt    >  2:2843707/1‑148 (MQ=255)
                                                                                                                                                 aCCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTc   <  1:957111/148‑1 (MQ=255)
                                                                                                                                                 aCCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTc   <  2:934863/148‑1 (MQ=255)
                                                                                                                                                  ccGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc  >  2:1283799/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCC  >  CP009273/3888127‑3888420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: