Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,806,160 A→T 9.8% S593T (TCG→ACG)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,806,1600AT 9.8% 257.8 / 16.7 92S593T (TCG→ACG) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (49/34);  new base T (4/5);  total (53/39)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTC  >  NC_002947/4806013‑4806308
                                                                                                                                                   |                                                                                                                                                      
ccGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAAc                                                                                                                                                      <  1:1134204/150‑1 (MQ=255)
        ggCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAgg                                                                                                                                              >  2:794739/1‑150 (MQ=255)
        ggCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAgg                                                                                                                                              <  1:1194427/150‑1 (MQ=255)
         gCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacga                                                                                                                                             <  2:2289415/150‑1 (MQ=255)
         gCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacga                                                                                                                                             <  2:2268426/150‑1 (MQ=255)
              gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCg                                                                                                                                        >  1:772089/1‑150 (MQ=255)
              gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCg                                                                                                                                        >  1:772284/1‑150 (MQ=255)
               ccAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGc                                                                                                                                       <  2:2078724/150‑1 (MQ=255)
                  gCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGa                                                                                                                                    <  2:1664847/150‑1 (MQ=255)
                  gCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGa                                                                                                                                    <  2:3342/150‑1 (MQ=255)
                   cTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAAGTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCTTCGGCCGAACGGCCGAGGAAATCGCTGac                                                                                                                                   <  1:1066176/150‑1 (MQ=255)
                        cGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGAccacca                                                                                                                              >  1:1029216/1‑150 (MQ=255)
                          aCCCATCGAGGGCTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCaaa                                                                                                                            <  2:1318806/150‑1 (MQ=255)
                            ccAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGc                                                                                                                          >  2:2306229/1‑150 (MQ=255)
                            ccAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGc                                                                                                                          >  2:1833458/1‑150 (MQ=255)
                               gcgcGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTg                                                                                                                       >  1:1712820/1‑150 (MQ=255)
                                   gggTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAg                                                                                                                   >  1:2097208/1‑150 (MQ=255)
                                     gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                 <  2:5468/150‑1 (MQ=255)
                                     gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                 >  1:1939522/1‑150 (MQ=255)
                                      ttGCGCCGCGCGCCGTCCCGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAgatt                                                                                                                <  1:1264364/150‑4 (MQ=255)
                                       tGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTgg                                                                                                               <  2:575527/150‑1 (MQ=255)
                                          gcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACa                                                                                                            >  1:675672/1‑150 (MQ=255)
                                             gcgcGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGaccacc                                                                                                                               >  1:1049406/1‑128 (MQ=255)
                                             gcgcGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGaccacc                                                                                                                               <  2:1049406/128‑1 (MQ=255)
                                                      ccTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCa                                                                                                <  2:1029216/150‑1 (MQ=255)
                                                       ttGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCGGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                 >  2:1151420/2‑132 (MQ=255)
                                                       ttGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCGGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                 <  1:1151420/131‑1 (MQ=255)
                                                       cTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAg                                                                                               >  2:481275/1‑150 (MQ=255)
                                                                   gctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGat                                                                                   >  1:1503294/1‑150 (MQ=255)
                                                                     tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                                 >  1:1235263/1‑150 (MQ=255)
                                                                     tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                                 >  1:611007/1‑150 (MQ=255)
                                                                      taTCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTg                                                                                        <  1:1920755/141‑1 (MQ=255)
                                                                      gATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTg                                                                                        >  2:1920755/1‑142 (MQ=255)
                                                                       aTCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATgg                                                                               <  2:1486771/150‑1 (MQ=255)
                                                                           ggTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCa                                                                         <  1:2306229/150‑1 (MQ=255)
                                                                             tACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                      >  1:2022673/1‑105 (MQ=255)
                                                                             tACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                      <  2:2022673/105‑1 (MQ=255)
                                                                                gggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTg                                                                    >  1:2419835/1‑150 (MQ=255)
                                                                                 ggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGa                                                                   >  2:2212333/1‑150 (MQ=255)
                                                                                   tttCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTg                                                                                                                <  2:1846648/105‑1 (MQ=255)
                                                                                   tttCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTg                                                                                                                >  1:1846648/1‑105 (MQ=255)
                                                                                       cGCTGGTGCTCAATTGCAGCTCTGGCCCTTCGGTGTCAGAACTCCGCACCTCTTCAGCCGAACGGCCGAGGCAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGc                                                             <  1:533423/150‑1 (MQ=255)
                                                                                       cGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGAtgctgc                                                             >  2:1039004/1‑150 (MQ=255)
                                                                                        gCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGCCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGcc                                                            >  2:83305/1‑150 (MQ=255)
                                                                                             tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCgg                                                       >  1:295136/1‑150 (MQ=255)
                                                                                             tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCgg                                                       >  1:995507/1‑150 (MQ=255)
                                                                                             tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCgg                                                       >  1:1286891/1‑150 (MQ=255)
                                                                                             tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCgc                                                       >  2:2088577/1‑149 (MQ=255)
                                                                                              gCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGgc                                                      >  1:516477/1‑150 (MQ=255)
                                                                                                tCAATTGCAGCTCTGGCCATTCGGTGGCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAgg                                                                       >  1:1638235/1‑131 (MQ=255)
                                                                                                tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAgg                                                                       <  2:1638235/131‑1 (MQ=255)
                                                                                                tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGgcgc                                                    >  1:804481/1‑150 (MQ=255)
                                                                                                tCAATTGCAGCTCTGGCCATTCGGTGACTGAAGTCTTCACACCGCCGGCCCAACGGCCAAGGCAATCGCTAACCAGAAAAAAATGGGGGGCCGACATGCCGGGGCCCAGCTCGGTGCCGATATGG‑‑CGCAAGCCGATGCGCCCAGcgccgc                                                    >  2:511456/1‑150 (MQ=255)
                                                                                                    ttGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGaa                                                <  2:239728/150‑1 (MQ=255)
                                                                                                       cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCACAGTGATGGAGGTGGACATGCCGGTGCACAGCTCGGTTCCGTTAGGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                             >  1:534611/1‑150 (MQ=255)
                                                                                                        aGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCg                                            >  1:1322747/1‑150 (MQ=255)
                                                                                                         gcgctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGa                                           <  2:2248678/147‑1 (MQ=255)
                                                                                                          ctctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGac                                          >  1:547594/1‑150 (MQ=255)
                                                                                                          ctctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCCGGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGACCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGac                                          >  1:1023565/1‑150 (MQ=255)
                                                                                                               gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                     <  1:1644242/150‑1 (MQ=255)
                                                                                                                  aTTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGc                                                                      <  2:442765/114‑1 (MQ=255)
                                                                                                                  aTTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGc                                                                      >  1:442765/1‑114 (MQ=255)
                                                                                                                       gTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                      <  2:115571/63‑1 (MQ=255)
                                                                                                                       gTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                      >  1:115571/1‑63 (MQ=255)
                                                                                                                        tGACTGCAGTCTTCACCTCTTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgcg                            <  2:1023565/150‑1 (MQ=255)
                                                                                                                          aCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGcc                          >  2:747652/1‑150 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCTGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                        >  2:2024737/1‑150 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                        >  1:1791352/1‑150 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                        <  2:675672/150‑1 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                        >  2:934287/1‑150 (MQ=255)
                                                                                                                               aGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAACGCTGTTAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGa                     >  1:1879302/1‑150 (MQ=255)
                                                                                                                                 tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc                   >  1:1436462/1‑150 (MQ=255)
                                                                                                                                   ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                 <  2:995507/150‑1 (MQ=255)
                                                                                                                                   ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCCAGCCTGGTACGGTTCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                 >  1:1002217/1‑150 (MQ=255)
                                                                                                                                    tCACCTCGTCGGCTGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGAGGGCATGGCGGTGACCAGCTCGCTGCCGATAAGG‑‑CGCAGGCTGATTCGGCCCTCGGCGCCGGACTCGACAATGGCCGTCGACGGCGCgggaacgg                >  2:870708/1‑144 (MQ=255)
                                                                                                                                     cACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGt               >  1:2074319/1‑150 (MQ=255)
                                                                                                                                     cACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGt               >  2:642512/1‑150 (MQ=255)
                                                                                                                                      aCCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTa              <  1:709397/150‑1 (MQ=255)
                                                                                                                                        ttcctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCg                                                               <  1:100373/95‑1 (MQ=255)
                                                                                                                                        ttcctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCg                                                               >  2:100373/5‑99 (MQ=255)
                                                                                                                                        cTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCg                                                                                     <  1:1154864/79‑1 (MQ=255)
                                                                                                                                        cTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCACAGCTCGGTGCCg                                                                                     >  2:1154864/1‑79 (MQ=255)
                                                                                                                                           gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa         <  2:2377341/150‑1 (MQ=255)
                                                                                                                                           gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa         <  1:794739/150‑1 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAgg                                                                       <  2:1601440/87‑1 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAgg                                                                       >  1:1601440/2‑88 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc                   <  1:1098737/139‑1 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc                   >  2:1098737/2‑140 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa         >  1:1154286/2‑150 (MQ=255)
                                                                                                                                               gCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCTATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTcc     <  2:1286891/150‑1 (MQ=255)
                                                                                                                                               gCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTcc     <  2:2097208/150‑1 (MQ=255)
                                                                                                                                                  gAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTc  >  1:2130178/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGG‑‑CGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTC  >  NC_002947/4806013‑4806308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: