Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 925,355 A→T 13.3% L355Q (CTG→CAG)  cydD ← fused glutathione, cysteine exporter subunits of ABC superfamily: membrane component/ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273925,3550AT13.3% 432.9 / 52.6 158L355Q (CTG→CAG) cydDfused glutathione, cysteine exporter subunits of ABC superfamily: membrane component/ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base A (50/86);  new base T (13/8);  total (63/95)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 160 total aligned reads displayed.

CCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGAT  >  CP009273/925208‑925502
                                                                                                                                                   |                                                                                                                                                   
cccTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACa                                                                                                                                                     <  2:3225347/148‑1 (MQ=255)
 ccTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAg                                                                                                                                                    >  1:2046082/1‑148 (MQ=255)
 ccTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAg                                                                                                                                                    >  1:1823972/1‑148 (MQ=255)
    gATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGctc                                                                                                                                                 >  1:32005/1‑148 (MQ=255)
      tatGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGctcct                                                                                                                                               <  1:1166836/148‑1 (MQ=255)
      tatGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGctcct                                                                                                                                               >  2:3034632/1‑148 (MQ=255)
         gagagaAAACCAGAAAGCGCGTTCAGCCGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAATGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCgg                                                                                                                                            >  2:1312962/1‑148 (MQ=255)
          agagaAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGc                                                                                                                                           <  2:437177/148‑1 (MQ=255)
            agaAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCt                                                                                                                                         <  2:2471105/148‑1 (MQ=255)
                aaCCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAAt                                                                                                                                     >  1:1192921/1‑148 (MQ=255)
                 aCCAGAAAGCGCGTTCAGCAGTGAGCTTTTACATGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATg                                                                                                                                    >  1:781206/1‑148 (MQ=255)
                   cAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGt                                                                                                                                  >  2:1284721/1‑148 (MQ=255)
                    aGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTc                                                                                                                                 <  2:41766/148‑1 (MQ=255)
                      aaaGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCAc                                                                                                                               >  1:2043165/1‑148 (MQ=255)
                      aaaGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCAc                                                                                                                               >  2:1064297/1‑148 (MQ=255)
                          cgcgTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGa                                                                                                                           <  2:1675987/148‑1 (MQ=255)
                           gcgTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGAt                                                                                                                          <  1:2878094/148‑1 (MQ=255)
                               tCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAACGGCCCGGGCCGCGGGCCGGCCCGCGTTTTACCTTCGGGGGTGGTGTTTTTGTGGTCGTGGGCCTCAATGGTCACCGGATCGGt                                                                                                                      >  1:1516243/1‑148 (MQ=255)
                                  gcagTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGa                                                                                                                   <  2:2476632/148‑1 (MQ=255)
                                   cagTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGAt                                                                                                                  >  1:3169377/1‑148 (MQ=255)
                                    agTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATg                                                                                                                 <  1:1034798/148‑1 (MQ=255)
                                       gAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTa                                                                                                              <  1:1262524/148‑1 (MQ=255)
                                         gCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAAt                                                                                                            <  2:3189117/148‑1 (MQ=255)
                                           ttttACCTGAACCGCTGCGACCAACCAACACCGGACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTc                                                                                                          <  2:1877919/148‑1 (MQ=255)
                                               aCCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGcc                                                                                                      <  1:1161845/148‑1 (MQ=255)
                                                ccTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCt                                                                                                     <  2:2324735/148‑1 (MQ=255)
                                                 cTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTc                                                                                                    >  2:2067581/1‑148 (MQ=255)
                                                    aaCCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTcacc                                                                                                 >  1:3164155/1‑148 (MQ=255)
                                                      ccGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTcaccac                                                                                               >  2:727272/1‑148 (MQ=255)
                                                         cTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGtt                                                                                            >  2:2653465/1‑148 (MQ=255)
                                                         cTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGtt                                                                                            >  2:2653625/1‑148 (MQ=255)
                                                          tGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTg                                                                                           <  2:32005/148‑1 (MQ=255)
                                                           gCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGc                                                                                          >  2:2520497/1‑148 (MQ=255)
                                                                caaccaacACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATg                                                                                     <  2:2277451/148‑1 (MQ=255)
                                                                caaccaacACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATg                                                                                     <  2:77641/148‑1 (MQ=255)
                                                                    caacACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCg                                                                                 <  2:921936/148‑1 (MQ=255)
                                                                      acacCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAg                                                                               <  2:1516243/148‑1 (MQ=255)
                                                                      acacCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAg                                                                               >  1:392962/1‑148 (MQ=255)
                                                                          cGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCggg                                                                           <  1:2983517/148‑1 (MQ=255)
                                                                          cGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCggg                                                                           <  1:1338836/148‑1 (MQ=255)
                                                                            cACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATATACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGt                                                                         <  1:122200/148‑1 (MQ=255)
                                                                            cACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGt                                                                         <  1:122476/148‑1 (MQ=255)
                                                                             aCGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGtt                                                                        <  1:1647104/148‑1 (MQ=255)
                                                                              cGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGttt                                                                       <  2:1062302/148‑1 (MQ=255)
                                                                                 tGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCa                                                                    >  2:1740964/1‑148 (MQ=255)
                                                                                         tgGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGtt                                                            <  1:1088911/148‑1 (MQ=255)
                                                                                          ggCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGttt                                                           <  2:457650/148‑1 (MQ=255)
                                                                                            cAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTc                                                         >  2:2542947/1‑148 (MQ=255)
                                                                                            cAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTc                                                         >  2:2542871/1‑148 (MQ=255)
                                                                                              aagtaaagtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAg                                                       <  2:2135403/148‑1 (MQ=255)
                                                                                               agtaaagtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGa                                                      <  2:1617579/148‑1 (MQ=255)
                                                                                               agtaaagtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGa                                                      <  2:1617577/148‑1 (MQ=255)
                                                                                               agtaaagtTCAGCGGTCCGGCCAGCGTTTTACCGTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTGTTCAGa                                                      <  2:3255739/148‑1 (MQ=255)
                                                                                                 taaagtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACt                                                    >  1:770411/1‑148 (MQ=255)
                                                                                                  aaagtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTg                                                   <  1:190792/148‑1 (MQ=255)
                                                                                                    agtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                 <  1:2198296/148‑1 (MQ=255)
                                                                                                    agtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                 >  2:1118242/1‑148 (MQ=255)
                                                                                                    agtTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                 >  1:1843856/1‑148 (MQ=255)
                                                                                                                  gCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAg                                   <  1:510191/148‑1 (MQ=255)
                                                                                                                  gCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAg                                   >  1:1022256/1‑148 (MQ=255)
                                                                                                                   ccAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGc                                  <  1:2122359/148‑1 (MQ=255)
                                                                                                                   ccAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGc                                  <  2:2773689/148‑1 (MQ=255)
                                                                                                                   ccAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGc                                  <  1:3186762/148‑1 (MQ=255)
                                                                                                                       cGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTgg                              >  1:1875628/1‑148 (MQ=255)
                                                                                                                       cGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTgg                              >  1:2638768/1‑148 (MQ=255)
                                                                                                                            tACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCttt                         <  2:1176072/148‑1 (MQ=255)
                                                                                                                            tACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCttt                         <  2:2331958/148‑1 (MQ=255)
                                                                                                                            tACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCttt                         <  2:2332148/148‑1 (MQ=255)
                                                                                                                             aCCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTa                        >  2:3009163/1‑148 (MQ=255)
                                                                                                                             aCCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTa                        >  2:3009575/1‑148 (MQ=255)
                                                                                                                             aCCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTa                        <  2:2517023/148‑1 (MQ=255)
                                                                                                                             aCCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTa                        <  2:2517166/148‑1 (MQ=255)
                                                                                                                              ccTTCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAg                       <  1:2125235/148‑1 (MQ=255)
                                                                                                                              ccTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAg                       >  2:2871919/1‑148 (MQ=255)
                                                                                                                              ccTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAg                       >  2:1841876/1‑148 (MQ=255)
                                                                                                                              ccTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAg                       >  1:1723724/1‑148 (MQ=255)
                                                                                                                               cTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGc                      <  2:1192921/148‑1 (MQ=255)
                                                                                                                                 tCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCAt                    >  2:3101568/1‑148 (MQ=255)
                                                                                                                                 tCCGGCGACGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCAt                    >  2:3101702/1‑148 (MQ=255)
                                                                                                                                 tCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCAt                    <  1:1319648/148‑1 (MQ=255)
                                                                                                                                 tCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCAt                    >  2:1362372/1‑148 (MQ=255)
                                                                                                                                 tCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCAt                    >  1:2931235/1‑148 (MQ=255)
                                                                                                                                   cGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGa                  <  1:1835131/148‑1 (MQ=255)
                                                                                                                                    ggCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGAt                 <  1:2809383/148‑1 (MQ=255)
                                                                                                                                        aCGTGATAAACAGCTCCTCGGCCTCACTGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATaaaa             <  1:2566928/148‑1 (MQ=255)
                                                                                                                                        aCGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATaaaa             >  2:766866/1‑148 (MQ=255)
                                                                                                                                        aCGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATaaaa             <  1:156288/148‑1 (MQ=255)
                                                                                                                                         cGTGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAAc            >  1:917452/1‑148 (MQ=255)
                                                                                                                                         cGTGATAAACAGCTCCTCGGCCTGAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGTTAAAAc            <  2:943307/148‑1 (MQ=255)
                                                                                                                                          gTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACg           <  2:2935493/148‑1 (MQ=255)
                                                                                                                                          gTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACg           <  2:446823/148‑1 (MQ=255)
                                                                                                                                           tGATAAACTGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGt          <  1:2653625/148‑1 (MQ=255)
                                                                                                                                            gATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTa         >  2:2125977/1‑148 (MQ=255)
                                                                                                                                            gATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTa         >  1:2427625/1‑148 (MQ=255)
                                                                                                                                            gATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTa         >  1:2426303/1‑148 (MQ=255)
                                                                                                                                             aTAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTAc        <  2:1794440/148‑1 (MQ=255)
                                                                                                                                               aaaCAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCg      >  2:140148/1‑148 (MQ=255)
                                                                                                                                               aaaCAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCg      >  2:249495/1‑148 (MQ=255)
                                                                                                                                                 aCAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCgag    <  1:2525896/148‑1 (MQ=255)
                                                                                                                                                   aGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGAt  >  1:503765/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTTGGCCTGCTGGCAAAGTAAAGTTCAGCGGTCCGGCCAGCGTTTTACCTTCCGGCGACGTGATAAACAGCTCCTCGGCCTCAATGGTCACCGGATCGGTCGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGAT  >  CP009273/925208‑925502

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: