Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,823 +A pseudogene (1174/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,8221.A99.0% 1614.8 / 5.8 401pseudogene (1173/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base A (189/208);  total (192/209)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 427 total aligned reads displayed.

GGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAC  >  NC_008463/3566673‑3566966
                                                                                                                                                      |                                                                                                                                                
   ggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aat                                                                                                                                               <  1:5318295/150‑3 (MQ=255)
    gTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aatc                                                                                                                                              <  2:3285974/150‑4 (MQ=255)
        aCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtcg                                                                                                                                           <  1:5039055/150‑1 (MQ=255)
        aCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtcg                                                                                                                                           >  1:474454/1‑150 (MQ=255)
         cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGa                                                                                                                                          >  2:3690180/1‑150 (MQ=255)
         cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGa                                                                                                                                          >  1:279967/1‑150 (MQ=255)
            cTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                       >  1:2927434/1‑150 (MQ=255)
            cTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                       <  1:1978628/150‑1 (MQ=255)
             tGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTgg                                                                                                                                      <  2:4570294/150‑1 (MQ=255)
             tGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTgg                                                                                                                                      >  2:8499414/1‑150 (MQ=255)
              ggCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGt                                                                                                                                     <  1:10413944/150‑1 (MQ=255)
                cgcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTAc                                                                                                                                   >  1:837358/1‑150 (MQ=255)
                   agcgcgCCGGAGCGCGCGGCCAACGTGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAGGTGGGAGAACTTCCTGATCCTGCGCGACCCCCTCTGCCGCGGCCAGTCGGACCCGTGCGGAATCGTCGAGTGGTAgcgc                                                                                                                                >  2:6906884/3‑150 (MQ=255)
                   aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgc                                                                                                                                >  2:1433930/1‑150 (MQ=255)
                          cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAcc                                                                                                                         <  1:8845621/150‑1 (MQ=255)
                          cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAcc                                                                                                                         >  1:6138875/1‑150 (MQ=255)
                           ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                                                                                        >  1:4291431/1‑150 (MQ=255)
                           ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                                                                                        >  2:528877/1‑150 (MQ=255)
                           ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                                                                                        >  2:3372765/1‑150 (MQ=255)
                            gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAt                                                                                                                       >  1:8030383/1‑150 (MQ=255)
                             aGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                                                                                                      >  2:1130651/1‑150 (MQ=255)
                             aGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                                                                                                      <  2:63092/150‑1 (MQ=255)
                              gcgcgcGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGc                                                                                                                     >  1:6327064/1‑150 (MQ=255)
                                gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGccc                                                                                                                   <  1:2552766/150‑1 (MQ=255)
                                 cgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCg                                                                                                                  <  2:2154186/150‑1 (MQ=255)
                                  gcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                                                                                                                 >  1:60455/1‑150 (MQ=255)
                                    ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgct                                                                                                               >  1:9491170/1‑150 (MQ=255)
                                    ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgct                                                                                                               <  1:646812/150‑1 (MQ=255)
                                     gCCATCGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                              <  2:9056501/150‑1 (MQ=255)
                                     gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                              >  2:6964554/1‑150 (MQ=255)
                                     gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                              >  1:8262478/1‑150 (MQ=255)
                                     gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                              >  1:8755737/1‑150 (MQ=255)
                                     gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                              >  1:3521046/1‑150 (MQ=255)
                                      ccAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAAGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcg                                                                                                             <  2:3501353/150‑1 (MQ=255)
                                       caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgc                                                                                                            >  2:4632929/1‑150 (MQ=255)
                                         acgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgctc                                                                                                          <  1:2237701/147‑1 (MQ=255)
                                          cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTcc                                                                                                         <  1:3714550/148‑1 (MQ=255)
                                          cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTcc                                                                                                         <  1:1097759/148‑1 (MQ=255)
                                           gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACGTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg                                                                                                        <  2:3315165/149‑1 (MQ=255)
                                            aGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg                                                                                                       >  1:9917717/1‑150 (MQ=255)
                                                gcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAc                                                                                                   >  1:3936480/1‑150 (MQ=255)
                                                gcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAGGCCCGGCTCGCTCCGGAGAc                                                                                                   <  2:2779719/150‑1 (MQ=255)
                                                   ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGa                                                                                                >  1:7248425/1‑150 (MQ=255)
                                                    cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt                                                                                               >  1:7663980/1‑150 (MQ=255)
                                                      cTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgc                                                                                             >  2:4008621/1‑150 (MQ=255)
                                                         gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgc                                                                                          >  2:1069010/2‑150 (MQ=255)
                                                          caGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                         <  1:10534103/149‑1 (MQ=255)
                                                            gTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCaa                                                                                       <  2:2263886/150‑1 (MQ=255)
                                                              tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                     <  1:3330359/150‑1 (MQ=255)
                                                              tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAGGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                     <  1:9389343/150‑1 (MQ=255)
                                                                ggCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                                                   >  2:4791587/1‑150 (MQ=255)
                                                                  cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcc                                                                                 >  1:7105413/1‑150 (MQ=255)
                                                                  cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcc                                                                                 <  1:9739981/150‑1 (MQ=255)
                                                                  cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcc                                                                                 >  1:3838855/1‑150 (MQ=255)
                                                                  cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcc                                                                                 <  1:5562511/150‑1 (MQ=255)
                                                                         accgagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACacc                                                                          <  2:8959468/144‑1 (MQ=255)
                                                                           cgagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaa                                                                        <  1:8593086/146‑1 (MQ=255)
                                                                           cgagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaa                                                                        <  2:9536512/146‑1 (MQ=255)
                                                                                ccTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACcaaccaac                                                                   <  1:10023068/150‑1 (MQ=255)
                                                                                 cTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                                                  <  2:6804764/150‑1 (MQ=255)
                                                                                 cTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                                                  >  1:183862/1‑150 (MQ=255)
                                                                                    gggCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGc                                                               >  2:5527398/1‑150 (MQ=255)
                                                                                    gggCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGc                                                               >  2:7697383/1‑150 (MQ=255)
                                                                                     ggCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCa                                                              <  2:10518428/150‑1 (MQ=255)
                                                                                      gCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCAt                                                             >  1:9984534/1‑150 (MQ=255)
                                                                                      gCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCAt                                                             >  1:1636468/1‑150 (MQ=255)
                                                                                        aTGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATcc                                                           <  2:10466043/150‑1 (MQ=255)
                                                                                         tgtgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATccc                                                          >  1:3179041/1‑150 (MQ=255)
                                                                                           tgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAg                                                        >  2:3783428/1‑150 (MQ=255)
                                                                                            gggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAgg                                                       >  2:118784/1‑150 (MQ=255)
                                                                                            gggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAgg                                                       >  2:132844/1‑150 (MQ=255)
                                                                                             ggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTCCCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGc                                                      >  2:822656/1‑150 (MQ=255)
                                                                                                 aaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTAcc                                                  <  2:708757/150‑1 (MQ=255)
                                                                                                   cTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGc                                                <  1:5864843/150‑1 (MQ=255)
                                                                                                      ccTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTg                                             <  2:31001/150‑1 (MQ=255)
                                                                                                       cTGATCCTGCGCGCCCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGATACCGCCTGa                                            <  2:5733378/150‑1 (MQ=255)
                                                                                                       cTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCCCACCAACCAACCCGCATCCCAGGCTACCGCCTGa                                            <  1:10219885/150‑1 (MQ=255)
                                                                                                         gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATc                                          >  2:7465859/1‑150 (MQ=255)
                                                                                                           tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcac                                        >  1:10897620/1‑150 (MQ=255)
                                                                                                             cTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCGGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcaccc                                      >  2:10289220/1‑148 (MQ=255)
                                                                                                               gcgcGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAg                                    <  1:6234188/150‑1 (MQ=255)
                                                                                                                  cgaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGaa                                 >  2:643348/3‑150 (MQ=255)
                                                                                                                        ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCa                           <  1:6898135/150‑1 (MQ=255)
                                                                                                                         ctctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACCCCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAt                          <  1:1071827/150‑1 (MQ=255)
                                                                                                                          tctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATg                         <  1:7391982/150‑1 (MQ=255)
                                                                                                                           ctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGa                        <  1:9110596/150‑1 (MQ=255)
                                                                                                                             gCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACCCCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAAt                      >  2:1941097/1‑150 (MQ=255)
                                                                                                                             gCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAAt                      >  2:8759062/1‑150 (MQ=255)
                                                                                                                             gCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAAt                      >  2:9225877/1‑150 (MQ=255)
                                                                                                                              ccGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATa                     <  1:3099538/150‑1 (MQ=255)
                                                                                                                                 cgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATTCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTc                  <  2:2031078/150‑1 (MQ=255)
                                                                                                                                 cgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTc                  <  2:7921731/150‑1 (MQ=255)
                                                                                                                                    ccAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGa               <  1:4713888/150‑1 (MQ=255)
                                                                                                                                    ccAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGa               >  2:2395341/1‑150 (MQ=255)
                                                                                                                                     cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAt              >  1:7878406/1‑150 (MQ=255)
                                                                                                                                     cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAt              >  2:6014431/1‑150 (MQ=255)
                                                                                                                                     cAGTCGGACCCGTTCGGAATCGTCGACTGGGACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAt              <  2:5087006/150‑1 (MQ=255)
                                                                                                                                       gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg            >  1:128524/1‑150 (MQ=255)
                                                                                                                                          ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGccc         >  2:163717/1‑150 (MQ=255)
                                                                                                                                               cGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATc    >  2:9625383/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAC  >  NC_008463/3566673‑3566966

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: ‑