Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,701 T→G 6.9% L19R (CTG→CGG)  PP_1087 → OmpA family outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,245,7010TG 6.9% 228.8 / 6.5 72L19R (CTG→CGG) PP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base G (2/3);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AAGCCCCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAG  >  NC_002947/1245553‑1245837
                                                                                                                                                    |                                                                                                                                        
aaGCCCCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCt                                                                                                                                          >  1:166621/1‑149 (MQ=255)
aaGCCCCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCt                                                                                                                                          >  3:95961/1‑149 (MQ=255)
    cccAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGATCGCCAGCCTgccg                                                                                                                                      <  3:2166/149‑1 (MQ=255)
     ccAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccgg                                                                                                                                     >  1:23284/1‑149 (MQ=255)
        aTCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccggccg                                                                                                                                  <  7:265991/149‑1 (MQ=255)
           cTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGgcg                                                                                                                               >  4:122446/1‑149 (MQ=255)
            tACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGgcg                                                                                                                               >  7:260165/1‑148 (MQ=255)
            tACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGgcg                                                                                                                               >  6:127831/1‑148 (MQ=255)
            tACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGgcg                                                                                                                               >  4:142149/1‑148 (MQ=255)
            tACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTGCGCCGGCCGGCAGGCTGgcg                                                                                                                               >  5:139042/1‑148 (MQ=255)
             aCCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCa                                                                                                                             >  5:108134/1‑149 (MQ=255)
               cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAca                                                                                                                           <  8:260165/149‑2 (MQ=255)
               cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAca                                                                                                                           <  3:66028/149‑2 (MQ=255)
                 gAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGa                                                                                                                         <  6:185872/149‑1 (MQ=255)
                 gAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTGCGCCGGCCGGCAGGCTGGCGAAcagg                                                                                                                         <  6:139042/149‑4 (MQ=255)
                  aaaaCCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGAc                                                                                                                        >  3:82965/1‑149 (MQ=255)
                  aaaaCCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCAGTTCGCCAGCCTGCCGGCCGGCGCACTGAc                                                                                                                        <  6:232288/149‑1 (MQ=255)
                    aaCCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACtt                                                                                                                      <  2:181197/149‑1 (MQ=255)
                     aCCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACttt                                                                                                                     <  7:161116/149‑1 (MQ=255)
                       cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTcc                                                                                                                   <  2:23284/149‑1 (MQ=255)
                       cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                    >  5:263258/1‑148 (MQ=255)
                       cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                    >  4:69211/1‑148 (MQ=255)
                                 gCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCAt                                                                                                         <  2:46282/149‑1 (MQ=255)
                                     ggCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAc                                                                                                               <  8:282742/139‑1 (MQ=255)
                                     ggCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAc                                                                                                               >  7:282742/1‑139 (MQ=255)
                                          tCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACAt                                                                                                >  7:5568/1‑149 (MQ=255)
                                            aaGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTg                                                                                              >  8:3718/1‑149 (MQ=255)
                                             aGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGc                                                                                             <  4:95961/149‑1 (MQ=255)
                                               ttttCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCt                                                                                           <  5:127831/149‑1 (MQ=255)
                                               ttttCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCt                                                                                           <  7:226363/149‑1 (MQ=255)
                                                           cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               >  2:195674/1‑149 (MQ=255)
                                                           cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               >  7:142339/1‑149 (MQ=255)
                                                           cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               <  2:166621/149‑1 (MQ=255)
                                                                gATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAc                                                                          <  3:69211/149‑1 (MQ=255)
                                                                    gCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAg                                                                       <  8:142339/148‑1 (MQ=255)
                                                                     ccTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGGGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTTCCTGGAAGGTCGAGGGTGACCAGtt                                                                     <  2:113922/149‑1 (MQ=255)
                                                                         ctgGTTAGCCGGCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGCAAGGTCGAGGGTGACCAGTTCg                                                                   <  7:198837/147‑1 (MQ=255)
                                                                         ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGa                                                                  <  3:142149/148‑1 (MQ=255)
                                                                            gTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgt                                                                <  4:82965/147‑1 (MQ=255)
                                                                                   gTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGa                                                       <  6:108134/149‑1 (MQ=255)
                                                                                    tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               >  7:219241/1‑124 (MQ=255)
                                                                                    tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               <  8:219241/124‑1 (MQ=255)
                                                                                    tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAcc                                                                         >  7:22023/1‑130 (MQ=255)
                                                                                    tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAcc                                                                         <  8:22023/130‑1 (MQ=255)
                                                                                        aaTTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagccag                                                  >  4:185968/1‑149 (MQ=255)
                                                                                         aTTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagccagc                                                 <  7:3718/149‑1 (MQ=255)
                                                                                                   ccgccCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAgg                                       >  6:138709/1‑149 (MQ=255)
                                                                                                     gccCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGtt                                     >  5:293085/1‑149 (MQ=255)
                                                                                                       ccGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTc                                   <  1:195674/149‑1 (MQ=255)
                                                                                                         gTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCgg                                 >  2:120837/1‑149 (MQ=255)
                                                                                                         gTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCgg                                 <  1:148724/149‑1 (MQ=255)
                                                                                                              ccAGCGTTGCCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            >  3:121918/1‑149 (MQ=255)
                                                                                                              ccAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            >  4:177849/1‑149 (MQ=255)
                                                                                                                  cGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTAc                        <  3:177849/149‑1 (MQ=255)
                                                                                                                    ttACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAggg                                                                              <  1:183077/93‑1 (MQ=255)
                                                                                                                    ttACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAggg                                                                              >  2:183077/1‑93 (MQ=255)
                                                                                                                    ttACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGt                      <  6:263258/149‑1 (MQ=255)
                                                                                                                      aCCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGtt                     <  8:5568/148‑1 (MQ=255)
                                                                                                                        cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                  <  3:185968/149‑1 (MQ=255)
                                                                                                                        cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                  <  3:122446/149‑1 (MQ=255)
                                                                                                                            cccTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGt              >  4:55616/1‑149 (MQ=255)
                                                                                                                                 ttGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCggg          >  7:134717/1‑148 (MQ=255)
                                                                                                                                   gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            >  7:115263/1‑128 (MQ=255)
                                                                                                                                   gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            <  8:115263/128‑1 (MQ=255)
                                                                                                                                    ccccTGTTCGCCAGCCTGCCGGCCGCCGCACTGACTTTCCAGACGCGCAGGGAGAACATTGCCTGGACGGGCGAGGGTGACCAGTTCGAGTGACGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAgg      <  6:263515/148‑1 (MQ=255)
                                                                                                                                     cccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc     <  5:138709/149‑1 (MQ=255)
                                                                                                                                      ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATGGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc     >  7:216862/1‑148 (MQ=255)
                                                                                                                                      ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            >  1:6240/1‑125 (MQ=255)
                                                                                                                                      ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            <  2:6240/125‑1 (MQ=255)
                                                                                                                                       cTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc     >  4:249651/1‑147 (MQ=255)
                                                                                                                                       cTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   <  6:293085/149‑1 (MQ=255)
                                                                                                                                        tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGag  >  2:202408/1‑149 (MQ=255)
                                                                                                                                          ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   <  1:32067/146‑1 (MQ=255)
                                                                                                                                          ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   >  2:32067/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
AAGCCCCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAG  >  NC_002947/1245553‑1245837

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