Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,167 C→A 5.5% G168C (GGC→TGC)  PP_1829 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,054,1670CA 5.5% 351.5 / 7.1 109G168C (GGC→TGC) PP_1829alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (54/49);  new base A (3/3);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 117 total aligned reads displayed.

GAAACACCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGC  >  NC_002947/2054023‑2054314
                                                                                                                                                |                                                                                                                                                   
gAAACACCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGc                                                                                                                                                >  1:617387/1‑150 (MQ=255)
 aaaCACCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGcc                                                                                                                                               >  1:153702/1‑150 (MQ=255)
   acacCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCga                                                                                                                                             <  2:15335/150‑2 (MQ=255)
    cacCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCAGCGGCCGGCCGCTggccgag                                                                                                                                           <  1:390732/150‑7 (MQ=255)
       cTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCAGCggccggccgc                                                                                                                                         >  2:960595/1‑148 (MQ=255)
       cTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCAGCggccggccgc                                                                                                                                         <  1:960595/150‑3 (MQ=255)
        tGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAgg                                                                                                                                        >  1:814985/1‑150 (MQ=255)
        tGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAgg                                                                                                                                        <  2:1036524/150‑1 (MQ=255)
            gcAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGtga                                                                                                                                    >  1:2286949/1‑150 (MQ=255)
             cAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGtgat                                                                                                                                   <  1:895091/150‑1 (MQ=255)
                cggcgCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGAc                                                                                                                                >  1:961398/1‑150 (MQ=255)
                     cTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTgcgc                                                                                                                           <  2:1035631/150‑1 (MQ=255)
                          acaTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAggg                                                                                                                                       >  1:1760134/1‑133 (MQ=255)
                          acaTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAggg                                                                                                                                       <  2:1760134/133‑1 (MQ=255)
                          acaTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGgcg                                                                                                                      >  1:2186298/1‑150 (MQ=255)
                           caTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTggccgag                                                                                                                                           <  1:1444283/128‑7 (MQ=255)
                            aTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAgg                                                                                                                                        >  2:841957/1‑130 (MQ=255)
                            aTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAgg                                                                                                                                        <  1:841957/130‑1 (MQ=255)
                            aTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGcc                                                                                                                    >  2:1524703/1‑150 (MQ=255)
                            aTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGcc                                                                                                                    >  1:409901/1‑150 (MQ=255)
                                  gCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTg                                                                                                              >  2:536152/1‑150 (MQ=255)
                                   cTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGc                                                                                                             <  2:1348626/150‑1 (MQ=255)
                                     ggTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTc                                                                                                           >  1:918675/1‑150 (MQ=255)
                                        cTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGGCCCAAGGGAAGCGTTTGTTCGTGCAGCACCGCGGCAAGCGGCAGGACGCTGCCGGGCTGTCTAATGACTTCGCTGGCGCCAGCCTGCTcggc                                                                                                        <  1:2351173/150‑1 (MQ=255)
                                         tGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGTGGCCGGCAGCggccggccgct                                                                                                                                        <  2:2096215/117‑4 (MQ=255)
                                         tGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGTGGCCGGCAGCggccggccgct                                                                                                                                        >  1:2096215/1‑114 (MQ=255)
                                           cccGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCCGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCggcggc                                                                                                     >  1:1234535/1‑150 (MQ=255)
                                           cccGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCggcggc                                                                                                     >  1:1561175/1‑150 (MQ=255)
                                                 cTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCt                                                                                               >  2:2374055/1‑150 (MQ=255)
                                                  tGCATGCGGTCGAGCAGGGCACTGCCAACGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTg                                                                                              <  2:509766/150‑1 (MQ=255)
                                                     aGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCa                                                                                           >  1:2233479/1‑150 (MQ=255)
                                                      ggCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGcc                                                                                                    <  2:59854/140‑1 (MQ=255)
                                                      ggCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGcc                                                                                                    >  1:59854/1‑140 (MQ=255)
                                                         ggTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTcgg                                                                                                         >  1:187549/1‑132 (MQ=255)
                                                         ggTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTcgg                                                                                                         <  2:187549/132‑1 (MQ=255)
                                                         ctacGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCa                                                                                       >  2:787338/4‑150 (MQ=255)
                                                          gTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAg                                                                                      >  1:595239/1‑150 (MQ=255)
                                                          gTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAg                                                                                      >  1:644832/1‑150 (MQ=255)
                                                          gTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAg                                                                                      <  2:2255772/150‑1 (MQ=255)
                                                          gTCGAGCAGGGCACGGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTgg                                                                                                                        <  1:1115992/116‑1 (MQ=255)
                                                           tCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAAt                                                                                         <  2:2211636/146‑1 (MQ=255)
                                                             gAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTgc                                                                                   >  2:2310493/1‑150 (MQ=255)
                                                               gCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTgcgc                                                                                 >  2:2327336/1‑150 (MQ=255)
                                                                cAGGGCACTGCCACCGGCCAGTTGCGTGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGcc                                                                                <  1:1699904/150‑1 (MQ=255)
                                                                cAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTg                                                                                                                         >  2:1448693/1‑109 (MQ=255)
                                                                cAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTg                                                                                                                         <  1:1448693/109‑1 (MQ=255)
                                                                    gCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCg                                                                            <  2:1037510/150‑1 (MQ=255)
                                                                    gCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCg                                                                            <  2:595239/150‑1 (MQ=255)
                                                                        tGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCAGCggccggccg                                                                                                                                          <  1:1106295/84‑2 (MQ=255)
                                                                        tGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCAGCggccggccg                                                                                                                                          >  2:1106295/1‑83 (MQ=255)
                                                                          ccACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTgg                                                                      >  1:2315881/1‑150 (MQ=255)
                                                                          ccACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTgg                                                                      >  2:217610/1‑150 (MQ=255)
                                                                          ccACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTgg                                                                      <  2:409901/150‑1 (MQ=255)
                                                                          ccACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCCTCGCTGTgg                                                                      <  2:2233479/150‑1 (MQ=255)
                                                                                gCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATc                                                                <  1:521475/150‑1 (MQ=255)
                                                                                  cAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGGCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAg                                                              >  1:2254018/1‑150 (MQ=255)
                                                                                   aGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTCGCCAATCAgg                                                             <  2:1886334/150‑1 (MQ=255)
                                                                                    gTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTgcgc                                                                                                                           >  1:520856/1‑87 (MQ=255)
                                                                                    gTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTgcgc                                                                                                                           <  2:520856/87‑1 (MQ=255)
                                                                                     ttGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGAt                                                           <  1:217610/150‑1 (MQ=255)
                                                                                          cgggcggCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGAGCCAGCCTGCTCGCCCGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGTTCATCg                                                      <  2:340869/150‑1 (MQ=255)
                                                                                            ggcggcAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCggg                                                    <  2:1374584/150‑1 (MQ=255)
                                                                                               ggcAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCg                                                 >  2:140988/1‑150 (MQ=255)
                                                                                                 cAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGaa                                               >  2:2259618/1‑150 (MQ=255)
                                                                                                     cgcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAgcgc                                           <  2:948853/150‑1 (MQ=255)
                                                                                                      gcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGaa                                               >  2:2209160/1‑145 (MQ=255)
                                                                                                      gcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGaa                                               <  1:2209160/145‑1 (MQ=255)
                                                                                                      gcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGcgcg                                          >  1:1559416/1‑150 (MQ=255)
                                                                                                      gcTGGGCCAGGTGTTCGTGCAGCACCTCGGCCAGCGGGCGGCCGCTGCCGGGCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGACAGGCTGGCAATCAGTGCGCCTGCGCTGTGGCCAATCAGGATCAGCGGGCCGAAgcgct                                          <  2:935934/150‑2 (MQ=255)
                                                                                                          ggCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTc                                      >  2:1688912/1‑150 (MQ=255)
                                                                                                          ggCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTc                                      >  2:6254/1‑150 (MQ=255)
                                                                                                            ccAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCgg                                    >  2:1805761/1‑150 (MQ=255)
                                                                                                             cAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGt                                                                                     <  1:1735522/100‑1 (MQ=255)
                                                                                                             cAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGt                                                                                     >  2:1735522/1‑100 (MQ=255)
                                                                                                               gttgttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCtt                                 >  2:1311088/1‑150 (MQ=255)
                                                                                                                 tgttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTg                               <  1:138388/150‑1 (MQ=255)
                                                                                                                  gttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGgcgc                                                                                                                     >  2:2235035/1‑63 (MQ=255)
                                                                                                                  gttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGgcgc                                                                                                                     <  1:2235035/63‑1 (MQ=255)
                                                                                                                   ttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCGTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAg                             >  1:117837/1‑150 (MQ=255)
                                                                                                                   ttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAg                             <  2:293361/150‑1 (MQ=255)
                                                                                                                      gggcagcgCCTGGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCGTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCtt                          <  2:117837/148‑1 (MQ=255)
                                                                                                                       tGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAgg                                                                                                 <  1:992510/78‑1 (MQ=255)
                                                                                                                       tGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTg                         <  2:2363224/150‑1 (MQ=255)
                                                                                                                        gcagcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAACGCGCGTCGCCTTTGAGCTTGc                        <  2:644832/150‑1 (MQ=255)
                                                                                                                         cagcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCg                       <  2:1221496/150‑1 (MQ=255)
                                                                                                                          agcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgtg                                                                       >  1:1046276/1‑101 (MQ=255)
                                                                                                                          agcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgtg                                                                       <  2:1046276/101‑1 (MQ=255)
                                                                                                                              ccTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTcc                  <  1:1805761/150‑1 (MQ=255)
                                                                                                                               cTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCa                 >  2:2016570/1‑150 (MQ=255)
                                                                                                                               cTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCa                 >  1:335452/1‑150 (MQ=255)
                                                                                                                                    ccAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGa            <  1:2149354/150‑1 (MQ=255)
                                                                                                                                    ccAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGa            <  1:787338/150‑1 (MQ=255)
                                                                                                                                    accctggccggcaGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCttt                                <  2:824632/125‑1 (MQ=255)
                                                                                                                                     cAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTACGAAACAGAAGCGCCCAATCAGGATGAGCGGGCCGAAGACCGGGCCGGCTTTGAGCTTGCGGCTCCAGGCGAc           >  1:1251343/1‑150 (MQ=255)
                                                                                                                                      aGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCATGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGTTCGGCTTTGAGCTTGCGGCTCCAGGCGAcc          >  1:1947936/1‑150 (MQ=255)
                                                                                                                                            cggccgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTTAGCTTGCGGCTCCAGGCGACCACATCg    <  1:1768512/150‑1 (MQ=255)
                                                                                                                                            cggccgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCg    >  2:2030247/1‑150 (MQ=255)
                                                                                                                                             ggccgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGTATCAGCGGGCCGAAGCACGCGTCGGCTTTGAGCTTGCAGCTCCAGGCCACCACATctg   >  2:219816/1‑148 (MQ=255)
                                                                                                                                                cgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCtt                          <  2:418495/124‑1 (MQ=255)
                                                                                                                                                cgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCtt                          >  1:418495/1‑124 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GAAACACCTGGCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGC  >  NC_002947/2054023‑2054314

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