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,2850GT97.8% 573.2 / 3.7 181D871Y (GAC→TAC) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (79/98);  total (81/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 183 total aligned reads displayed.

CGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATC  >  CP009273/3017142‑3017428
                                                                                                                                               |                                                                                                                                               
cgcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGc                                                                                                                                             >  2:32411/1‑148 (MQ=255)
    ttctggTGATATCCTCCTGGAGAAATTCTACGCCTTGCTGTACTTCGTTACACATGTCAATCCAGAGTTAGCGTTGGGTCAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGTAGAGATCATTTATTACGCCGaa                                                                                                                                         <  2:3038398/144‑1 (MQ=255)
     aCGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAg                                                                                                                                        >  1:1209954/1‑148 (MQ=255)
       gggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGt                                                                                                                                      <  1:2008117/148‑1 (MQ=255)
        ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTc                                                                                                                                     >  2:2256705/1‑148 (MQ=255)
        ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTc                                                                                                                                     <  2:2334629/148‑1 (MQ=255)
        ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTc                                                                                                                                     >  2:464059/1‑148 (MQ=255)
          tGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAc                                                                                                                                   >  1:2466268/1‑148 (MQ=255)
          tGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAc                                                                                                                                   <  2:1897458/148‑1 (MQ=255)
           gAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAcc                                                                                                                                  >  1:3553350/1‑148 (MQ=255)
           gAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAcc                                                                                                                                  <  2:1092357/148‑1 (MQ=255)
            aaaTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAccc                                                                                                                                 >  2:2801866/1‑148 (MQ=255)
              aTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGt                                                                                                                               <  1:3477559/148‑1 (MQ=255)
                ccGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTa                                                                                                                             <  1:3280456/148‑1 (MQ=255)
                ccGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTa                                                                                                                             <  1:2136347/148‑1 (MQ=255)
                ccGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTa                                                                                                                             >  2:2748675/1‑148 (MQ=255)
                 cGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTaa                                                                                                                            <  1:1533650/148‑1 (MQ=255)
                       gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTg                                                                                                                      >  1:2765984/1‑148 (MQ=255)
                        aTAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGa                                                                                                                     <  1:192461/148‑1 (MQ=255)
                        aTAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGa                                                                                                                     <  1:1237212/148‑1 (MQ=255)
                        aTAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGa                                                                                                                     <  2:2368531/148‑1 (MQ=255)
                         tAAATTCTTCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCAGAGTATGAGCGAAGAGATCATTTATTTCGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                    <  2:3276942/148‑1 (MQ=255)
                          aaaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAtt                                                                                                                   <  2:1978948/148‑1 (MQ=255)
                          aaaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAtt                                                                                                                   >  1:3210137/1‑148 (MQ=255)
                          aaaTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAtt                                                                                                                   >  1:3210274/1‑148 (MQ=255)
                             ttCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTAc                                                                                                                >  2:1004990/1‑148 (MQ=255)
                             ttCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTAc                                                                                                                >  2:1005137/1‑148 (MQ=255)
                              tCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACg                                                                                                               >  1:80567/1‑148 (MQ=255)
                              tCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACg                                                                                                               >  1:80322/1‑148 (MQ=255)
                                tACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCa                                                                                                             <  2:684269/148‑1 (MQ=255)
                                tACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCa                                                                                                             <  2:684268/148‑1 (MQ=255)
                                 aCGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAg                                                                                                            >  1:2746856/1‑148 (MQ=255)
                                   gCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGtt                                                                                                          <  1:965501/148‑1 (MQ=255)
                                    ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTa                                                                                                         <  1:133303/148‑1 (MQ=255)
                                     cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                        >  1:416998/1‑148 (MQ=255)
                                     cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                        <  1:1439525/148‑1 (MQ=255)
                                         cTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgcg                                                                                                    >  1:3333624/1‑148 (MQ=255)
                                            gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAc                                                                                                 >  1:1045775/1‑148 (MQ=255)
                                            gACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAc                                                                                                 >  1:283986/1‑148 (MQ=255)
                                              cTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCg                                                                                               >  1:1679018/1‑148 (MQ=255)
                                               tGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGa                                                                                              >  2:748098/1‑148 (MQ=255)
                                                  ggTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGaaaa                                                                                           <  1:3214546/148‑1 (MQ=255)
                                                   gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAAt                                                                                          <  1:742285/148‑1 (MQ=255)
                                                     acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTg                                                                                        <  1:3178567/148‑1 (MQ=255)
                                                           gTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACa                                                                                  <  2:3523675/148‑1 (MQ=255)
                                                            tCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAt                                                                                 >  2:1524417/1‑148 (MQ=255)
                                                             cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                >  2:1924881/1‑148 (MQ=255)
                                                              aaTCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTc                                                                               >  2:2708396/1‑148 (MQ=255)
                                                                tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                             <  2:977482/148‑1 (MQ=255)
                                                                tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                             <  2:51830/148‑1 (MQ=255)
                                                                tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                             >  2:242876/1‑148 (MQ=255)
                                                                   agagTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGt                                                                          >  2:1025521/1‑148 (MQ=255)
                                                                   agagTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGt                                                                          >  1:3438152/1‑148 (MQ=255)
                                                                     agTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGa                                                                        >  1:1568439/1‑148 (MQ=255)
                                                                     agTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGa                                                                        <  2:1285/148‑1 (MQ=255)
                                                                       ttAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGAtt                                                                      <  2:2409193/148‑1 (MQ=255)
                                                                        tAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGAttt                                                                     <  2:3553350/148‑1 (MQ=255)
                                                                             ttGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcgcg                                                                >  1:3271963/1‑148 (MQ=255)
                                                                                 gACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtgt                                                            >  2:2538208/1‑148 (MQ=255)
                                                                                 gACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtgt                                                            >  2:1943329/1‑148 (MQ=255)
                                                                                 gACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtgt                                                            >  2:1600112/1‑148 (MQ=255)
                                                                                 gACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTAACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCtgt                                                            >  2:1600192/1‑148 (MQ=255)
                                                                                  aCAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACTTTCCGCGTGATTTCCGCGCtgtg                                                           <  2:1072506/148‑1 (MQ=255)
                                                                                      aTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTgg                                                       >  1:1412397/1‑148 (MQ=255)
                                                                                      aTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTgg                                                       >  1:3186183/1‑148 (MQ=255)
                                                                                        cTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg                                                     >  1:2257575/1‑148 (MQ=255)
                                                                                         tACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGc                                                    <  1:1760507/148‑1 (MQ=255)
                                                                                                ccACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcga                                             >  2:111844/1‑148 (MQ=255)
                                                                                                   cccTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACCGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTCGTGCCGAGagacga                                          >  1:2743631/1‑148 (MQ=255)
                                                                                                       gcatGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                      <  1:729043/145‑1 (MQ=255)
                                                                                                       gcatGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                      <  1:3088135/145‑1 (MQ=255)
                                                                                                       gCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaaa                                      <  2:3210176/148‑1 (MQ=255)
                                                                                                        cGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAg                                     >  2:450282/1‑148 (MQ=255)
                                                                                                          tGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTc                                   >  2:1866180/1‑148 (MQ=255)
                                                                                                          tGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTc                                   >  2:1747429/1‑148 (MQ=255)
                                                                                                           gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  <  1:1461250/148‑1 (MQ=255)
                                                                                                           gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  >  1:43641/1‑148 (MQ=255)
                                                                                                           gCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                  <  2:784316/148‑1 (MQ=255)
                                                                                                             tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGc                                <  1:3028701/148‑1 (MQ=255)
                                                                                                             tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGc                                <  1:2658817/148‑1 (MQ=255)
                                                                                                                cGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCg                             <  2:1679018/148‑1 (MQ=255)
                                                                                                                 ggCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGt                            <  1:1585019/148‑1 (MQ=255)
                                                                                                                  gCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGtt                           <  2:416998/148‑1 (MQ=255)
                                                                                                                    cacaGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                         <  1:1469775/148‑1 (MQ=255)
                                                                                                                    cacaGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCg                         <  1:748098/148‑1 (MQ=255)
                                                                                                                     acaGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCgg                        >  2:3094958/1‑148 (MQ=255)
                                                                                                                       aGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCgggg                      >  2:3147855/1‑148 (MQ=255)
                                                                                                                          aTGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGa                   <  2:3118304/148‑1 (MQ=255)
                                                                                                                            gAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                 <  1:2792620/148‑1 (MQ=255)
                                                                                                                               cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg              >  2:1665007/1‑148 (MQ=255)
                                                                                                                                 aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt            >  2:1790518/1‑148 (MQ=255)
                                                                                                                                  agagATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc           <  2:3394318/148‑1 (MQ=255)
                                                                                                                                   gagaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct          <  1:1005032/148‑1 (MQ=255)
                                                                                                                                   gagaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct          <  1:1004990/148‑1 (MQ=255)
                                                                                                                                    agaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc         >  1:3101572/1‑148 (MQ=255)
                                                                                                                                    agaTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc         >  1:1326912/1‑148 (MQ=255)
                                                                                                                                        cATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGaa     <  2:2265286/148‑1 (MQ=255)
                                                                                                                                        cATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGaa     <  2:80567/148‑1 (MQ=255)
                                                                                                                                          tttATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt   >  1:170979/1‑148 (MQ=255)
                                                                                                                                           ttattaCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATc  >  2:2898311/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATC  >  CP009273/3017142‑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: