Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,751 A→G pseudogene (1102/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,7510AG99.8% 1713.0 / NA 432pseudogene (1102/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (203/228);  minor base T (0/1);  total (203/229)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 432 total aligned reads displayed.

CGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACC  >  NC_008463/3566603‑3566894
                                                                                                                                                    |                                                                                                                                                
cGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                                                                                                               >  2:8400715/1‑151 (MQ=255)
   cATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                                                                                                            >  1:4682353/1‑151 (MQ=255)
    aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                                                                                                           >  1:3604256/1‑151 (MQ=255)
       ccGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCa                                                                                                                                        <  1:97813/151‑1 (MQ=255)
        cGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAt                                                                                                                                       >  2:4654061/1‑151 (MQ=255)
         gCTGCGGGTGTTCTACTTCATGCTGGACCTGCCCAGCCACGACTTCCTCCACCGCAAGCATGGGTTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGTCAAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                                                                                                                      <  2:4216624/151‑1 (MQ=255)
          cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                                                                                                                     >  1:1268277/1‑151 (MQ=255)
          cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                                                                                                                     >  2:262949/1‑151 (MQ=255)
           tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgtg                                                                                                                                    >  1:7755038/1‑151 (MQ=255)
            gCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTgg                                                                                                                                   >  2:6143683/1‑151 (MQ=255)
             cGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTggg                                                                                                                                  >  1:6473699/1‑151 (MQ=255)
              gggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                                                                                                                 <  2:3109692/151‑1 (MQ=255)
              gggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                                                                                                                 >  2:253284/1‑151 (MQ=255)
                  gtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAGGTGGGagacc                                                                                                                             >  2:9877108/1‑147 (MQ=255)
                    tctcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACtt                                                                                                                           <  1:5559666/151‑1 (MQ=255)
                       ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCCCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCt                                                                                                                        <  2:763090/151‑1 (MQ=255)
                          tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGag                                                                                                                     >  2:932248/1‑150 (MQ=255)
                          tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                                                                                                                     >  1:8977342/1‑151 (MQ=255)
                          tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                                                                                                                     <  1:4012052/151‑1 (MQ=255)
                          tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                                                                                                                     <  2:5601215/151‑1 (MQ=255)
                          tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                                                                                                                     <  1:6068771/151‑1 (MQ=255)
                           cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                                                                                                                    <  1:6280628/151‑1 (MQ=255)
                            aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                                                                                                                   >  1:7997937/1‑151 (MQ=255)
                             tGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                                                                                                                  <  2:3388720/151‑1 (MQ=255)
                             tGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                                                                                                                  <  1:10380373/151‑1 (MQ=255)
                               cAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                                                                                                                <  1:2946776/151‑1 (MQ=255)
                                  gACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                                                                                                             >  2:1204014/1‑150 (MQ=255)
                                   aCCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcga                                                                                                            >  2:9469989/1‑149 (MQ=255)
                                    ccTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGac                                                                                                           <  1:7665525/151‑1 (MQ=255)
                                      tGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAGCCTGCGCGacca                                                                                                         >  2:2968634/1‑151 (MQ=255)
                                       gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                                                                                                        >  1:3951790/1‑151 (MQ=255)
                                       gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                                                                                                        <  1:4464176/151‑1 (MQ=255)
                                        cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                                                                                                       >  1:2673275/1‑151 (MQ=255)
                                        cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                                                                                                       >  2:4114618/1‑151 (MQ=255)
                                         ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct                                                                                                      <  1:8639994/151‑1 (MQ=255)
                                          cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctc                                                                                                     >  1:3119917/1‑151 (MQ=255)
                                          cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctc                                                                                                     <  2:7264601/151‑1 (MQ=255)
                                          cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctc                                                                                                     <  1:3297982/151‑1 (MQ=255)
                                          cAGCCACGACTTCCACCCCCGCCAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAAGTTCCTGATCCTGCGCGACCACctc                                                                                                     >  2:6533329/1‑151 (MQ=255)
                                           aGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctct                                                                                                    <  2:6286839/151‑1 (MQ=255)
                                             ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc                                                                                                  <  2:8034833/151‑1 (MQ=255)
                                              cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc                                                                                                 >  1:10021363/1‑151 (MQ=255)
                                               aCGACTTCCACCACCGCAAGCCGGGGGTGAACTACTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcg                                                                                                <  1:5478467/151‑1 (MQ=255)
                                                 gACTTCCACCACCGCAAGCCGGGGGTGACCGTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGCTGCCCGAGACCTGGGGCATGTGGGAGAACTTCCTGAGCCTGCGCGACCACCTCTGCcgcg                                                                                              >  2:1616755/1‑151 (MQ=255)
                                                 gACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg                                                                                              >  2:3290486/1‑151 (MQ=255)
                                                   cTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGc                                                                                            <  1:5823710/151‑1 (MQ=255)
                                                    ttCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc                                                                                           >  1:8441002/1‑151 (MQ=255)
                                                     tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa                                                                                          <  1:134151/151‑1 (MQ=255)
                                                      ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAg                                                                                         <  2:7404858/151‑1 (MQ=255)
                                                       caccacCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGt                                                                                        <  1:9340090/151‑1 (MQ=255)
                                                         accaccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTc                                                                                       <  1:2383266/147‑1 (MQ=255)
                                                         accaccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTc                                                                                       >  1:5730923/5‑151 (MQ=255)
                                                           caccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCCTGGGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCgg                                                                                     >  2:1344161/3‑151 (MQ=255)
                                                             ccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAc                                                                                   <  1:1763799/151‑1 (MQ=255)
                                                                cAAGCCGGGGGTGACCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCAGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCg                                                                                >  2:2420049/1‑151 (MQ=255)
                                                                cAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCg                                                                                <  2:5362313/151‑1 (MQ=255)
                                                                cAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCg                                                                                <  2:10402859/151‑1 (MQ=255)
                                                                  aGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGtt                                                                              >  1:9613331/1‑151 (MQ=255)
                                                                  aGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGtt                                                                              <  2:7849691/151‑1 (MQ=255)
                                                                     cGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAGGTGGGAGACCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTGCgg                                                                           >  2:6003325/1‑151 (MQ=255)
                                                                      gggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGa                                                                         >  1:4145456/1‑151 (MQ=255)
                                                                       ggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaa                                                                        <  1:5931964/151‑2 (MQ=255)
                                                                         ggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatc                                                                      <  1:7388963/151‑4 (MQ=255)
                                                                               cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGAc                                                                 >  1:7057560/1‑151 (MQ=255)
                                                                                ttCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACt                                                                >  1:10709967/1‑151 (MQ=255)
                                                                                      cgcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTAcc                                                          <  1:915871/151‑1 (MQ=255)
                                                                                        cATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgc                                                        <  1:2977480/151‑1 (MQ=255)
                                                                                         aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTAGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgcg                                                       <  2:3987022/151‑1 (MQ=255)
                                                                                           cgcgCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgag                                                     >  1:8238200/1‑151 (MQ=255)
                                                                                             cgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGaa                                                   >  2:967296/1‑151 (MQ=255)
                                                                                              gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                  <  2:10652238/151‑1 (MQ=255)
                                                                                              gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                  <  1:2263885/151‑1 (MQ=255)
                                                                                              gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                  >  2:6790969/1‑151 (MQ=255)
                                                                                               ccGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAcc                                                 >  2:1075912/1‑151 (MQ=255)
                                                                                                cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                <  2:57873/151‑1 (MQ=255)
                                                                                                cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                >  2:2183999/1‑151 (MQ=255)
                                                                                                  gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                              <  2:10017799/151‑1 (MQ=255)
                                                                                                          ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                      >  2:3552200/1‑151 (MQ=255)
                                                                                                          ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                      >  1:10562337/1‑151 (MQ=255)
                                                                                                             caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                   >  2:4379004/1‑151 (MQ=255)
                                                                                                              aacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgctc                                  >  1:10101970/1‑151 (MQ=255)
                                                                                                              aacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgctc                                  >  2:8752462/1‑151 (MQ=255)
                                                                                                                 gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg                               <  1:192647/150‑1 (MQ=255)
                                                                                                                   ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGCCTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag                             <  2:2893400/151‑1 (MQ=255)
                                                                                                                     cgcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAc                           >  2:4874215/1‑151 (MQ=255)
                                                                                                                     cgcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAc                           <  2:8987860/151‑1 (MQ=255)
                                                                                                                      gcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAcc                          <  1:6857409/151‑1 (MQ=255)
                                                                                                                      gcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAcc                          <  1:1498445/151‑1 (MQ=255)
                                                                                                                               gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                 <  1:9453142/150‑1 (MQ=255)
                                                                                                                               gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                 <  1:5292610/150‑1 (MQ=255)
                                                                                                                               gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                 >  2:868710/2‑151 (MQ=255)
                                                                                                                                aaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCa                >  2:9634466/1‑151 (MQ=255)
                                                                                                                                   ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc             <  1:3393931/151‑1 (MQ=255)
                                                                                                                                      ggCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGc          >  1:4361648/1‑151 (MQ=255)
                                                                                                                                        cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGccc        >  2:4952623/1‑151 (MQ=255)
                                                                                                                                         ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcccc       >  2:7115616/1‑151 (MQ=255)
                                                                                                                                          cGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCca      <  2:10153886/151‑1 (MQ=255)
                                                                                                                                           gATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcac     >  2:6925249/1‑151 (MQ=255)
                                                                                                                                            aTGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCCCCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCCTCGACTCGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCACCCCcaca    <  2:2115107/151‑1 (MQ=255)
                                                                                                                                             tGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCacac   >  1:4999840/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACC  >  NC_008463/3566603‑3566894

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: