Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,888,274 C→A 43.8% 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,2740CA43.8% 142.9 / 185.3 173P88Q (CCG→CAG) chrRchromate reductase, Class I, flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (53/44);  new base A (49/25);  total (102/71)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 177 total aligned reads displayed.

CGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCC  >  CP009273/3888127‑3888420
                                                                                                                                                   |                                                                                                                                                  
cGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                                                                                                                    <  2:2036668/148‑1 (MQ=255)
   tACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAggg                                                                                                                                                 >  2:1485728/1‑148 (MQ=255)
     ccatccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTg                                                                                                                                               >  1:603065/1‑148 (MQ=255)
     ccatccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTg                                                                                                                                               <  2:1125338/148‑1 (MQ=255)
       atccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTggg                                                                                                                                             <  2:2559544/148‑1 (MQ=255)
        tccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                                                                                                                            <  2:3018002/148‑1 (MQ=255)
         ccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                                                                                                                           >  2:3295857/1‑148 (MQ=255)
         ccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                                                                                                                           <  2:2243688/148‑1 (MQ=255)
         ccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCt                                                                                                                                           >  1:2282926/1‑148 (MQ=255)
            ttGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaa                                                                                                                                        >  1:2747100/1‑148 (MQ=255)
            ttGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGaa                                                                                                                                        <  1:818313/148‑1 (MQ=255)
             tGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaa                                                                                                                                       <  1:543688/148‑1 (MQ=255)
                 gACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATg                                                                                                                                   <  1:304357/148‑1 (MQ=255)
                  aCATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTATTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGc                                                                                                                                  >  2:1737360/1‑148 (MQ=255)
                   cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                                                                                                                 >  1:1202149/1‑148 (MQ=255)
                   cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                                                                                                                 >  1:2852277/1‑148 (MQ=255)
                    aTTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                                                                                                                >  1:2562957/1‑148 (MQ=255)
                       cccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCACAGGGTGAATCGCGGGAGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCg                                                                                                                             >  1:3249396/1‑148 (MQ=255)
                       cccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCg                                                                                                                             >  1:817617/1‑148 (MQ=255)
                         cTTGTATGACGCTGACGTAAAGCTGGAAGAAGGTTTTCCTGAAACGGTTGAAGCTCTGGCGTAACAGTTCGGTCTGGCTGACGGTGTGGTGATGGTCCCGCCGGAATATACCTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGAc                                                                                                                           <  1:1102278/148‑1 (MQ=255)
                          ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                                                                                                                          >  2:414751/1‑148 (MQ=255)
                          ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGGGTGGTGATCGTCACGCCGGCCTATAACTACTCGGTACCGGGTGGGCTGAACAATGCCATCGACt                                                                                                                          >  2:412603/1‑148 (MQ=255)
                            gTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCTGGCTGAGGGTGTGGTGATCGTCACGCCGTAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTgg                                                                                                                        <  1:1597087/148‑1 (MQ=255)
                            gTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTgg                                                                                                                        <  2:162794/148‑1 (MQ=255)
                             tATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGc                                                                                                                       >  2:2657317/1‑148 (MQ=255)
                              aTGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCTGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCt                                                                                                                      <  2:967979/148‑1 (MQ=255)
                                gACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCttt                                                                                                                    >  1:1094844/1‑148 (MQ=255)
                                     tGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGc                                                                                                               >  2:615053/1‑148 (MQ=255)
                                            cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCgg                                                                                                        >  2:2841731/1‑148 (MQ=255)
                                             agcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGa                                                                                                       >  1:2266870/1‑148 (MQ=255)
                                               cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc                                                                                                     >  1:1849541/1‑148 (MQ=255)
                                                agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACCCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa                                                                                                    <  2:1534038/148‑1 (MQ=255)
                                                agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGACAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa                                                                                                    <  2:3172919/148‑1 (MQ=255)
                                                 tgaagaaGGTTTTCCAGCAACGTTTGAAGATCTGGGGTAACAGATCGGTCAGGGTGAGTGTTTGGTGACCGTCACGCCGGAATATAACTCCTGGGTACCGGGTGGGCTGAAAAATCCCATCGACTGGCTTTCCCGCCTGCCGGATCaa                                                                                                   <  1:1063319/147‑1 (MQ=255)
                                                 ggAAGAAGGTTTTCCAGCAACGTTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa                                                                                                   <  1:1063555/148‑1 (MQ=255)
                                                 ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa                                                                                                   >  1:2434941/1‑148 (MQ=255)
                                                 ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa                                                                                                   <  1:2152670/148‑1 (MQ=255)
                                                  gaacaaGGTTTGCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc                                                                                                  <  2:2959389/148‑1 (MQ=255)
                                                    agatgGTTTTCCACCAACGTTTTAAGCTCTGGCGGAACATATCCGTCAGGCTGAGGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCg                                                                                                <  1:3156158/148‑1 (MQ=255)
                                                    agaagGTTTTCCAGCAAGGGTTGAAGCTCTGGCGGAACAGATCCGTCTGGCTGAGGGTGTGGTGATCGTCCCGCCGTAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCg                                                                                                <  1:3156484/148‑1 (MQ=255)
                                                    agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCg                                                                                                >  1:1772206/1‑148 (MQ=255)
                                                     gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGc                                                                                               >  1:77830/1‑148 (MQ=255)
                                                     gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGc                                                                                               >  2:1531455/1‑148 (MQ=255)
                                                      aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCt                                                                                              >  1:2819452/1‑148 (MQ=255)
                                                        ggTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTgg                                                                                            >  1:252301/1‑148 (MQ=255)
                                                         gTTTTCCAGCAACGGTTGAAGCTCGGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGc                                                                                           <  2:625407/148‑1 (MQ=255)
                                                          ttttCCAGCACGGTTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGAGCGTCACGCGGTAATATAACTACTCGGTACCGGGTGGGCTGAAAAAGGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGcc                                                                                          <  1:1147466/148‑1 (MQ=255)
                                                               cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCTGGCTGACGGTGTGGTTATCGTCACGCCGTAATATAACTACTCGGTACTGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaa                                                                                     <  1:351308/148‑1 (MQ=255)
                                                               cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaa                                                                                     <  1:2856335/148‑1 (MQ=255)
                                                                aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTaaa                                                                                    >  1:2256750/1‑148 (MQ=255)
                                                                     cGGTTGAAGCTCTGGCGGAACAGAGCCGTCAGGCTGACGGTGGGGTGATCGTCACGCCGGAATATAACTCCTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGt                                                                               <  2:2855898/148‑1 (MQ=255)
                                                                      ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTa                                                                              <  2:3020911/148‑1 (MQ=255)
                                                                       gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTAt                                                                             >  2:1452095/1‑148 (MQ=255)
                                                                       gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTACACCGGTAt                                                                             >  2:1686373/1‑148 (MQ=255)
                                                                         tGAAGCTCGGGCGCAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTg                                                                           <  1:892392/148‑1 (MQ=255)
                                                                          gAAGCTCTGGCGGAACAGATCCGTCAGGCGGAGGGGGTGGTGATCGTCACGCCGGAATATAACTACGCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGa                                                                          <  1:1790018/148‑1 (MQ=255)
                                                                           aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAt                                                                         <  2:2747100/148‑1 (MQ=255)
                                                                           aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAt                                                                         <  2:1881124/148‑1 (MQ=255)
                                                                            aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAtt                                                                        <  2:1903904/148‑1 (MQ=255)
                                                                            aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAtt                                                                        <  2:234397/148‑1 (MQ=255)
                                                                             gCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTc                                                                       >  1:888907/1‑148 (MQ=255)
                                                                             gCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTc                                                                       >  1:2849134/1‑148 (MQ=255)
                                                                              ctctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCa                                                                      >  2:1176129/1‑148 (MQ=255)
                                                                               tctGGCGTAACAGATCCGTCTGGCTTAGGGTGTGGTTATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTTAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAg                                                                     <  2:3234831/148‑1 (MQ=255)
                                                                                ctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGa                                                                    <  2:3216042/148‑1 (MQ=255)
                                                                                ctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGa                                                                    >  1:1280533/1‑148 (MQ=255)
                                                                                ctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGa                                                                    >  1:1280984/1‑148 (MQ=255)
                                                                                 tGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGAc                                                                   >  1:1422899/1‑148 (MQ=255)
                                                                                   gCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCa                                                                 >  1:779510/1‑148 (MQ=255)
                                                                                      gAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCt                                                              >  1:3221046/1‑148 (MQ=255)
                                                                                         cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                           <  2:1266965/148‑1 (MQ=255)
                                                                                         cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                           <  1:3275801/148‑1 (MQ=255)
                                                                                         cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                           <  2:3249396/148‑1 (MQ=255)
                                                                                           gATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATg                                                         <  2:358148/148‑1 (MQ=255)
                                                                                               cGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCg                                                     >  1:1228188/1‑148 (MQ=255)
                                                                                                  cAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGa                                                  <  1:739078/148‑1 (MQ=255)
                                                                                                             tgtgGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcg                                       >  2:1681227/1‑148 (MQ=255)
                                                                                                               tggtgATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcgcg                                     <  2:2434941/148‑1 (MQ=255)
                                                                                                                     tCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCa                               >  2:1582099/1‑148 (MQ=255)
                                                                                                                        tCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTa                            <  1:2657317/148‑1 (MQ=255)
                                                                                                                         cACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTAt                           <  2:207500/148‑1 (MQ=255)
                                                                                                                           cGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCa                         >  2:1297628/1‑148 (MQ=255)
                                                                                                                           cGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCa                         >  2:2312662/1‑148 (MQ=255)
                                                                                                                              cGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCt                      >  1:3000793/1‑148 (MQ=255)
                                                                                                                                gAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                    >  1:1982435/1‑148 (MQ=255)
                                                                                                                                gAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                    >  2:1829159/1‑148 (MQ=255)
                                                                                                                                gAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                    >  2:1829158/1‑148 (MQ=255)
                                                                                                                                gAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                    >  2:1020308/1‑148 (MQ=255)
                                                                                                                                  atatAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgcgc                  >  2:814210/1‑148 (MQ=255)
                                                                                                                                    ataACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCACTGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCa                >  1:2290929/1‑148 (MQ=255)
                                                                                                                                    ataACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCa                >  1:2291073/1‑148 (MQ=255)
                                                                                                                                    ataACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCa                >  2:620103/1‑148 (MQ=255)
                                                                                                                                    ataACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCa                >  1:932158/1‑148 (MQ=255)
                                                                                                                                     taACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAg               <  2:2111581/148‑1 (MQ=255)
                                                                                                                                      aaCTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGa              >  1:1256383/1‑148 (MQ=255)
                                                                                                                                          actCGGTACCGGGGGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCt          <  2:460126/148‑1 (MQ=255)
                                                                                                                                           ctCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTg         <  1:2199957/148‑1 (MQ=255)
                                                                                                                                             cGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGt       <  2:77830/148‑1 (MQ=255)
                                                                                                                                                tACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGtttt    >  2:1060435/1‑148 (MQ=255)
                                                                                                                                                  ccGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTcc  <  2:2562957/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCC  >  CP009273/3888127‑3888420

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