Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,708 C→T 7.2% A21A (GCC→GCT)  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,7080CT 7.2% 200.7 / 4.6 69A21A (GCC→GCT) PP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/3);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

ATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAG  >  NC_002947/1245561‑1245837
                                                                                                                                                   |                                                                                                                                 
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                                                                                                                    >  4:69211/1‑148 (MQ=255)
               cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                    >  5:263258/1‑148 (MQ=255)
                         gCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCAt                                                                                                         <  2:46282/149‑1 (MQ=255)
                             ggCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAc                                                                                                               >  7:282742/1‑139 (MQ=255)
                             ggCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAc                                                                                                               <  8:282742/139‑1 (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                                                                               <  2:166621/149‑1 (MQ=255)
                                                   cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               >  7:142339/1‑149 (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                                                                               <  8:219241/124‑1 (MQ=255)
                                                                            tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                               >  7:219241/1‑124 (MQ=255)
                                                                            tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAcc                                                                         <  8:22023/130‑1 (MQ=255)
                                                                            tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAcc                                                                         >  7:22023/1‑130 (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                                                                              >  2:183077/1‑93 (MQ=255)
                                                                                                            ttACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAggg                                                                              <  1:183077/93‑1 (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                            <  2:6240/125‑1 (MQ=255)
                                                                                                                              ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                            >  1:6240/1‑125 (MQ=255)
                                                                                                                               cTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc     >  4:249651/1‑147 (MQ=255)
                                                                                                                               cTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   <  6:293085/149‑1 (MQ=255)
                                                                                                                                tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGag  >  2:202408/1‑149 (MQ=255)
                                                                                                                                  ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   >  2:32067/1‑146 (MQ=255)
                                                                                                                                  ttCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa   <  1:32067/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                 
ATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAG  >  NC_002947/1245561‑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: ‑