Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,017,285 G→T 100% D871Y (GAC→TAC)  xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,017,2850GT99.4% 577.8 / ‑4.2 177D871Y (GAC→TAC) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (82/94);  total (83/94)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 179 total aligned reads displayed.

CACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATC  >  CP009273/3017140‑3017428
                                                                                                                                                 |                                                                                                                                               
cacGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTAc                                                                                                                                               >  2:536439/1‑148 (MQ=255)
    cgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCg                                                                                                                                           <  2:1641795/148‑1 (MQ=255)
     gTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                          <  2:1118783/148‑1 (MQ=255)
      tACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGaa                                                                                                                                         >  1:8121/1‑148 (MQ=255)
      tACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGaa                                                                                                                                         >  2:2146192/1‑148 (MQ=255)
      tACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGaa                                                                                                                                         <  1:3104631/148‑1 (MQ=255)
          ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTc                                                                                                                                     <  2:2088385/148‑1 (MQ=255)
                aTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGt                                                                                                                               <  2:2546569/148‑1 (MQ=255)
                  ccGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTa                                                                                                                             >  2:1754973/1‑148 (MQ=255)
                   cGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTaa                                                                                                                            <  2:1855293/148‑1 (MQ=255)
                      cTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAcgc                                                                                                                         >  2:2891528/1‑148 (MQ=255)
                       tGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAcgcg                                                                                                                        >  1:1665837/1‑148 (MQ=255)
                         gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTg                                                                                                                      >  1:173593/1‑148 (MQ=255)
                         gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTg                                                                                                                      >  2:2159462/1‑148 (MQ=255)
                         gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTg                                                                                                                      >  2:2344384/1‑148 (MQ=255)
                         gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTg                                                                                                                      >  2:516471/1‑148 (MQ=255)
                          aTAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGa                                                                                                                     <  2:3458171/148‑1 (MQ=255)
                             aaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                  <  1:307421/148‑1 (MQ=255)
                                tCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACg                                                                                                               <  2:1062901/148‑1 (MQ=255)
                                   aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAg                                                                                                            >  1:1358897/1‑148 (MQ=255)
                                   aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAg                                                                                                            <  2:1333459/148‑1 (MQ=255)
                                   aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAg                                                                                                            >  1:1358898/1‑148 (MQ=255)
                                    cGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGt                                                                                                           >  1:1940599/1‑148 (MQ=255)
                                      ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTa                                                                                                         <  2:1683130/148‑1 (MQ=255)
                                       cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                  >  1:3380064/1‑138 (MQ=255)
                                       cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                  <  2:3380064/138‑1 (MQ=255)
                                       cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                        >  1:1400738/1‑148 (MQ=255)
                                       cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                        >  1:902807/1‑148 (MQ=255)
                                         tGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACgg                                                                                                      >  2:1902117/1‑148 (MQ=255)
                                         tGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACgg                                                                                                      >  2:1163862/1‑148 (MQ=255)
                                              gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAc                                                                                                 >  2:3064598/1‑148 (MQ=255)
                                                 tGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGa                                                                                              >  2:1505631/1‑148 (MQ=255)
                                                    ggTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGaaaa                                                                                           <  2:1961004/148‑1 (MQ=255)
                                                     gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAt                                                                                          <  2:1254039/148‑1 (MQ=255)
                                                     gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAt                                                                                          <  1:1925482/148‑1 (MQ=255)
                                                      tACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                         >  2:694698/1‑148 (MQ=255)
                                                      tACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                         >  2:694621/1‑148 (MQ=255)
                                                      tACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                         >  1:1718867/1‑148 (MQ=255)
                                                      tACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                         <  1:291384/148‑1 (MQ=255)
                                                      tACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAtt                                                                                         <  1:2086846/148‑1 (MQ=255)
                                                       acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTg                                                                                        <  1:211063/148‑1 (MQ=255)
                                                        cacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTgg                                                                                       >  2:3524161/1‑148 (MQ=255)
                                                        cacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTgg                                                                                       >  2:3524169/1‑148 (MQ=255)
                                                           cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGa                                                                                    <  1:1834992/148‑1 (MQ=255)
                                                           cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGa                                                                                    <  1:2416780/148‑1 (MQ=255)
                                                           cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGa                                                                                    <  1:2417307/148‑1 (MQ=255)
                                                            ggTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGAc                                                                                   <  2:1593441/148‑1 (MQ=255)
                                                             gTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACa                                                                                  <  1:3408577/148‑1 (MQ=255)
                                                               cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                <  1:101877/148‑1 (MQ=255)
                                                               cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                >  2:3207581/1‑148 (MQ=255)
                                                               cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                <  1:101813/148‑1 (MQ=255)
                                                               cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                <  1:369995/148‑1 (MQ=255)
                                                                  tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                             >  1:83128/1‑148 (MQ=255)
                                                                  tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                             >  1:1548464/1‑148 (MQ=255)
                                                                   ccAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcgc                                                                            >  1:3080013/1‑148 (MQ=255)
                                                                      gagTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTg                                                                         >  2:926874/1‑148 (MQ=255)
                                                                        gTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGAt                                                                       <  2:3293500/148‑1 (MQ=255)
                                                                           aGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTc                                                                    <  2:2692666/148‑1 (MQ=255)
                                                                           aGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTc                                                                    <  2:2691867/148‑1 (MQ=255)
                                                                            gCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTcc                                                                   >  2:1043033/1‑148 (MQ=255)
                                                                            gCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTcc                                                                   >  1:1282851/1‑148 (MQ=255)
                                                                             cGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcg                                                                  >  1:596071/1‑148 (MQ=255)
                                                                               ttGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTTTGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcgcg                                                                <  2:2895645/148‑1 (MQ=255)
                                                                               ttGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcgcg                                                                >  1:2512486/1‑148 (MQ=255)
                                                                               ttGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcgcg                                                                <  1:3193769/148‑1 (MQ=255)
                                                                                     cAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGc                                                          <  1:2997792/148‑1 (MQ=255)
                                                                                           tACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGc                                                    >  1:2818672/1‑148 (MQ=255)
                                                                                            aCGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGcc                                                   >  2:1721291/1‑148 (MQ=255)
                                                                                             cGGTGCCACCCTTCTTTCTATCGGACACAGTATGATCGAAGAGATCATTTTTTACGCCGAATGTCAGCCGTTATCGCGTTATTTACGCAGTTACTGCTCACCGAAAATTGGTTACATTCCGCGTGATTTCCGCTCTGTGCTGGTGCCg                                                  <  2:2825922/148‑1 (MQ=255)
                                                                                              ggTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGa                                                 <  1:3461219/148‑1 (MQ=255)
                                                                                               gTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAg                                                <  2:2410950/148‑1 (MQ=255)
                                                                                               gTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAg                                                >  1:1341591/1‑148 (MQ=255)
                                                                                                  ccACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcga                                             <  2:609947/148‑1 (MQ=255)
                                                                                                         gCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                      >  2:726821/1‑148 (MQ=255)
                                                                                                          cGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAg                                     <  2:1308092/148‑1 (MQ=255)
                                                                                                          cGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAg                                     <  2:29659/148‑1 (MQ=255)
                                                                                                           gTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGt                                    <  2:327404/148‑1 (MQ=255)
                                                                                                             gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  >  2:2411552/1‑148 (MQ=255)
                                                                                                             gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  >  1:267060/1‑148 (MQ=255)
                                                                                                             gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  >  2:1379148/1‑148 (MQ=255)
                                                                                                               tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGc                                <  2:2183688/148‑1 (MQ=255)
                                                                                                               tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGc                                <  1:2594776/148‑1 (MQ=255)
                                                                                                                aTCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGcc                               >  2:1091685/1‑148 (MQ=255)
                                                                                                                aTCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGcc                               <  1:584674/148‑1 (MQ=255)
                                                                                                                aTCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGcc                               <  1:584440/148‑1 (MQ=255)
                                                                                                                  cGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCg                             <  2:2437796/148‑1 (MQ=255)
                                                                                                                             tGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaa                  <  1:3442581/148‑1 (MQ=255)
                                                                                                                              gAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                 <  1:792920/148‑1 (MQ=255)
                                                                                                                              gAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                 >  1:1855596/1‑148 (MQ=255)
                                                                                                                               aGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAAt                <  2:3080013/148‑1 (MQ=255)
                                                                                                                                 cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg              >  1:949453/1‑148 (MQ=255)
                                                                                                                                   aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt            <  1:803092/148‑1 (MQ=255)
                                                                                                                                     gagaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct          >  2:2421778/1‑148 (MQ=255)
                                                                                                                                     gagaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct          >  1:961604/1‑148 (MQ=255)
                                                                                                                                       gaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCg        >  2:936508/1‑148 (MQ=255)
                                                                                                                                        aTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg       >  1:224260/1‑148 (MQ=255)
                                                                                                                                        aTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg       >  1:1936539/1‑148 (MQ=255)
                                                                                                                                            tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt   <  1:3352378/148‑1 (MQ=255)
                                                                                                                                            tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt   >  1:677435/1‑148 (MQ=255)
                                                                                                                                             ttattaCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATc  >  2:3430771/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATC  >  CP009273/3017140‑3017428

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: