Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,135,129 T→A 5.2% K407I (AAA→ATA)  hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,135,1290TA 5.2% 777.3 / 24.0 247K407I (AAA→ATA) hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base T (138/96);  new base A (10/3);  total (148/99)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 251 total aligned reads displayed.

TGTGATCGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCCGTCGT  >  CP009273/3134983‑3135276
                                                                                                                                                  |                                                                                                                                                   
tgtgATCGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGttt                                                                                                                                                    <  2:179699/148‑1 (MQ=255)
  tgATCGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCt                                                                                                                                                  >  1:1882839/1‑148 (MQ=255)
   gATCGCCTGTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTg                                                                                                                                                 <  2:3530844/148‑1 (MQ=255)
     tCGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGAt                                                                                                                                               >  2:1533450/1‑148 (MQ=255)
      cGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATa                                                                                                                                              >  2:272398/1‑148 (MQ=255)
       gCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATaa                                                                                                                                             <  2:2154288/148‑1 (MQ=255)
             gCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCa                                                                                                                                       <  1:1514385/148‑1 (MQ=255)
              ccGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCaa                                                                                                                                      <  2:758178/148‑1 (MQ=255)
                   atTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTATCTGATAAATCGCAACGAtt                                                                                                                                 >  2:2289780/1‑148 (MQ=255)
                    ttGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGAttt                                                                                                                                <  2:383610/148‑1 (MQ=255)
                      ggTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATttca                                                                                                                              <  2:2765349/148‑1 (MQ=255)
                       gTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATttcat                                                                                                                             >  1:3174704/1‑148 (MQ=255)
                         gATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATttcattc                                                                                                                           <  1:3141316/148‑1 (MQ=255)
                         gATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATttcattc                                                                                                                           <  1:1388652/148‑1 (MQ=255)
                          aTCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATttcattca                                                                                                                          >  1:3276159/1‑148 (MQ=255)
                           tCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAg                                                                                                                         <  1:1002170/148‑1 (MQ=255)
                             aGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGtt                                                                                                                       >  1:350054/1‑148 (MQ=255)
                               tGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTg                                                                                                                     <  2:1015122/148‑1 (MQ=255)
                               tGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTg                                                                                                                     <  2:1276087/148‑1 (MQ=255)
                                   cTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCt                                                                                                                 >  1:677711/1‑148 (MQ=255)
                                   cTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCt                                                                                                                 >  2:2897277/1‑148 (MQ=255)
                                     gTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTgg                                                                                                               >  1:1763099/1‑148 (MQ=255)
                                       agtGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGt                                                                                                             <  1:1481307/146‑1 (MQ=255)
                                        ttGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTa                                                                                                            >  1:3005270/1‑148 (MQ=255)
                                         tGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAg                                                                                                           >  2:317330/1‑148 (MQ=255)
                                         tGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAg                                                                                                           >  1:2798717/1‑148 (MQ=255)
                                         tGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAg                                                                                                           <  1:2016494/148‑1 (MQ=255)
                                          ggTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGa                                                                                                          <  1:1987645/148‑1 (MQ=255)
                                            ttttGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGAtt                                                                                                        >  2:1235851/1‑148 (MQ=255)
                                             tttGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTc                                                                                                       >  2:1870380/1‑148 (MQ=255)
                                              ttGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCTACGATTTCATTCAGTTTGTTCTGGGTAGATTcg                                                                                                      <  1:2444549/148‑1 (MQ=255)
                                              ttGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTATCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTcg                                                                                                      <  1:1703838/148‑1 (MQ=255)
                                                gCAGGAGATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTcgcg                                                                                                    <  2:2757235/148‑1 (MQ=255)
                                                    gATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCt                                                                                                >  2:1973353/1‑148 (MQ=255)
                                                     atatCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGGGACCTg                                                                                               >  2:227939/1‑148 (MQ=255)
                                                       atCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGAGTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTgcc                                                                                             >  2:3279465/1‑148 (MQ=255)
                                                        tCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTgccg                                                                                            >  1:1822510/1‑148 (MQ=255)
                                                           tGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCa                                                                                         >  1:2804701/1‑148 (MQ=255)
                                                             cAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGt                                                                                       >  2:1106890/1‑148 (MQ=255)
                                                             cAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGt                                                                                       <  1:573777/148‑1 (MQ=255)
                                                                ttCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTc                                                                                    >  2:2206817/1‑148 (MQ=255)
                                                                 tCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCa                                                                                   >  1:1617139/1‑148 (MQ=255)
                                                                 tCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCa                                                                                   >  1:1783425/1‑148 (MQ=255)
                                                                 tCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCa                                                                                   >  1:2429139/1‑148 (MQ=255)
                                                                      gcagTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGc                                                                              <  1:2727152/148‑1 (MQ=255)
                                                                       cagTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCa                                                                             <  2:1882839/148‑1 (MQ=255)
                                                                       cagTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCa                                                                             >  1:2843038/1‑148 (MQ=255)
                                                                         gTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCata                                                                           >  2:1337610/1‑148 (MQ=255)
                                                                          tGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCatat                                                                          >  2:3386997/1‑148 (MQ=255)
                                                                          tGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTGTCACCAGCatat                                                                          >  2:2013008/1‑148 (MQ=255)
                                                                          tGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTATCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCatat                                                                          >  1:987836/1‑148 (MQ=255)
                                                                            aaCGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTATCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTa                                                                        >  2:2448185/1‑148 (MQ=255)
                                                                               ggTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCTCTTCCAGCGTGTTTCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGCTTCGCGACCTGCCGCCAGTTTCACCAGCTTATTAGcc                                                                     <  2:2580138/148‑1 (MQ=255)
                                                                               ggTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGcc                                                                     <  2:3625767/148‑1 (MQ=255)
                                                                                gTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCa                                                                    <  1:2083081/148‑1 (MQ=255)
                                                                                     aCCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGc                                                               >  2:218687/1‑148 (MQ=255)
                                                                                     aCCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGc                                                               <  2:3319823/148‑1 (MQ=255)
                                                                                        aataatACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGcccc                                                            >  1:274116/1‑148 (MQ=255)
                                                                                         ataataCGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGAAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGACCGTGGCCCCa                                                           >  1:332387/1‑148 (MQ=255)
                                                                                            ataCGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACtt                                                        <  2:1421251/148‑1 (MQ=255)
                                                                                            ataCGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACtt                                                        <  1:1221145/148‑1 (MQ=255)
                                                                                              aCGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCt                                                      <  2:919690/148‑1 (MQ=255)
                                                                                              aCGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCt                                                      <  2:2524295/148‑1 (MQ=255)
                                                                                               cGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTa                                                     <  2:1283228/148‑1 (MQ=255)
                                                                                               cGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTa                                                     <  2:444920/148‑1 (MQ=255)
                                                                                                 gCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTAcc                                                   >  1:1679847/1‑148 (MQ=255)
                                                                                                  cccAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCg                                                  >  2:3506721/1‑148 (MQ=255)
                                                                                                    cAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGtt                                                <  1:2278607/148‑1 (MQ=255)
                                                                                                    cAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGtt                                                >  2:2327779/1‑148 (MQ=255)
                                                                                                     aGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGttt                                               >  1:3627579/1‑148 (MQ=255)
                                                                                                           gAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGt                                         >  1:695844/1‑148 (MQ=255)
                                                                                                                cAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGaaa                                    <  2:688786/148‑1 (MQ=255)
                                                                                                                cAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGaaa                                    <  1:317330/148‑1 (MQ=255)
                                                                                                                 aGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAg                                   >  2:3102799/1‑148 (MQ=255)
                                                                                                                   cTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTc                                 >  2:1189401/1‑148 (MQ=255)
                                                                                                                   cTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATACGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTc                                 >  1:786681/1‑148 (MQ=255)
                                                                                                                       gCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTg                             <  1:1123842/148‑1 (MQ=255)
                                                                                                                       gCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTg                             <  2:15475/148‑1 (MQ=255)
                                                                                                                         cACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGAt                           >  1:2390184/1‑148 (MQ=255)
                                                                                                                         cACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGAt                           >  1:3546036/1‑148 (MQ=255)
                                                                                                                          aCTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGAtt                          >  2:1302977/1‑148 (MQ=255)
                                                                                                                                aGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCAccc                    >  1:2686829/1‑148 (MQ=255)
                                                                                                                                 gCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTcacaca                   >  2:3397697/1‑148 (MQ=255)
                                                                                                                                 gCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCa                   >  2:3186779/1‑148 (MQ=255)
                                                                                                                                  cGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAg                  >  2:2440369/1‑148 (MQ=255)
                                                                                                                                    tgtTGCCAGTCAGTATCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGa                >  1:769381/1‑148 (MQ=255)
                                                                                                                                      ttGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAat              <  2:3255085/148‑1 (MQ=255)
                                                                                                                                       tGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAata             >  2:2655366/1‑148 (MQ=255)
                                                                                                                                        gCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAatat            <  1:2897302/148‑1 (MQ=255)
                                                                                                                                        gCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAatat            <  1:451070/148‑1 (MQ=255)
                                                                                                                                          cagtcagtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATAttt          <  2:274116/148‑1 (MQ=255)
                                                                                                                                           agtcagtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTc         <  1:8117/148‑1 (MQ=255)
                                                                                                                                            gtcagtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTcc        >  2:2945664/1‑148 (MQ=255)
                                                                                                                                             tcagtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTccc       >  1:557458/1‑148 (MQ=255)
                                                                                                                                              cagtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcg      >  1:1206063/1‑148 (MQ=255)
                                                                                                                                                gtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcgtc    >  1:3816955/1‑148 (MQ=255)
                                                                                                                                                gtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcgtc    >  1:860027/1‑148 (MQ=255)
                                                                                                                                                gtTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcgtc    >  1:2773979/1‑148 (MQ=255)
                                                                                                                                                 tttCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcgtcg   >  2:3458558/1‑148 (MQ=255)
                                                                                                                                                 tttCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCcgtcg   <  2:2529427/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTGATCGCCTTTGCCGATATTGGTGATCAGTGCACTGTATTGGTTTTGCAGGATATCCTGCAATTCACAGCAGTGAACGGTACGACCAATAATACGGCCCAGCGTGGAGTGCAGCTGTGCCACTTCCAGCGTGTTGCCAGTCAGTTTCTGATAAATCGCAACGATTTCATTCAGTTTGTTCTGGGTAGATTCGCGACCTGCCGCCAGTTTCACCAGCATATTAGCCAGTGGCCCCACTTCTACCGTTTTGCCGTAGAAAGTCGGTGATTTCACCCAGGAATATTTCCCGTCGT  >  CP009273/3134983‑3135276

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: