Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,704 C→A 10.3% P20Q (CCG→CAG)  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,7040CA10.3% 165.2 / 7.4 56P20Q (CCG→CAG) PP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base A (3/3);  total (27/29)
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.77e-01

CCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCA  >  NC_002947/1245558‑1245851
                                                                                                                                                  |                                                                                                                                                    
ccAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccgg                                                                                                                                                    >  8:280077/1‑149 (MQ=255)
  aaTCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccggcc                                                                                                                                                  >  7:123254/1‑149 (MQ=255)
         ccGGAAAACGCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAc                                                                                                                                           <  2:251023/149‑1 (MQ=255)
          cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACt                                                                                                                                          >  2:58238/1‑149 (MQ=255)
                 ccccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                                   <  7:280077/149‑1 (MQ=255)
                  cccGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTcc                                                                                                                                  <  6:177966/149‑1 (MQ=255)
                    cGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                                   >  8:196041/1‑146 (MQ=255)
                    cGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTc                                                                                                                                   >  3:284157/1‑146 (MQ=255)
                        gCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAcgc                                                                                                                            <  8:123254/149‑1 (MQ=255)
                          gAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAcgcg                                                                                                                           >  1:340482/1‑148 (MQ=255)
                                ggCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGgag                                                                                                                    <  4:284157/149‑1 (MQ=255)
                                      cAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCCGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACAtt                                                                                                              >  5:4079/1‑149 (MQ=255)
                                       aaGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTg                                                                                                             <  2:340482/149‑1 (MQ=255)
                                        aGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGc                                                                                                            >  8:154283/1‑149 (MQ=255)
                                             tCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTgg                                                                                                        <  5:71622/148‑1 (MQ=255)
                                              ccGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTgg                                                                                                        >  4:188074/1‑147 (MQ=255)
                                               cGGACAGCGTCCGATAGCCTGCTGGTAAGCCGGCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGaa                                                                                                      <  7:323079/148‑1 (MQ=255)
                                                ggACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAgg                                                                                                    >  2:45560/1‑149 (MQ=255)
                                                    aGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTc                                                                                                  >  4:52689/1‑147 (MQ=255)
                                                      cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                             <  4:158814/149‑1 (MQ=255)
                                                      cGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                             <  7:18739/149‑1 (MQ=255)
                                                           gATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAc                                                                                        <  6:4079/149‑1 (MQ=255)
                                                            aTAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGAc                                                                                        >  1:232598/1‑148 (MQ=255)
                                                             tAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCa                                                                                      <  3:52689/149‑1 (MQ=255)
                                                                    ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                             >  3:52111/1‑135 (MQ=255)
                                                                    ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAgg                                                                                             <  4:52111/135‑1 (MQ=255)
                                                                    ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAg                                                                               <  7:154283/149‑1 (MQ=255)
                                                                       gTAAGCCGTCAGATTTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgtgt                                                                            >  1:304893/1‑149 (MQ=255)
                                                                       gTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgtg                                                                             <  2:232598/148‑1 (MQ=255)
                                                                            ccGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCt                                                                       >  8:73830/1‑149 (MQ=255)
                                                                            ccGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCt                                                                       <  3:188074/149‑1 (MQ=255)
                                                                              gTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGa                                                                     >  8:287360/1‑149 (MQ=255)
                                                                               tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGGAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGag                                                                    >  8:202899/1‑149 (MQ=255)
                                                                                   aaTTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagccag                                                                >  4:30544/1‑149 (MQ=255)
                                                                                          ggAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATc                                                          >  1:7910/1‑148 (MQ=255)
                                                                                           gAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGa                                                        <  2:247285/149‑1 (MQ=255)
                                                                                            aGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       <  1:58238/149‑1 (MQ=255)
                                                                                                gccCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGtt                                                   >  4:310034/1‑149 (MQ=255)
                                                                                                 cccGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGttt                                                  <  7:196041/149‑1 (MQ=255)
                                                                                                  ccGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGGAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTc                                                 <  7:202899/149‑1 (MQ=255)
                                                                                                     tGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCgg                                               >  2:67656/1‑148 (MQ=255)
                                                                                                          cAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTTCCTGGAAGGTCGCAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                          <  1:90208/149‑1 (MQ=255)
                                                                                                           aGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtc                                                                           >  2:95909/1‑114 (MQ=255)
                                                                                                                  ccTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTc                                  <  7:287360/148‑1 (MQ=255)
                                                                                                                  ccTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCg                                 >  6:70256/1‑149 (MQ=255)
                                                                                                                   cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                >  3:327385/1‑149 (MQ=255)
                                                                                                                        ccTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTc                           <  3:113603/149‑1 (MQ=255)
                                                                                                                            ttGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCgg                         <  4:327385/147‑1 (MQ=255)
                                                                                                                              gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAg                     >  1:339079/1‑149 (MQ=255)
                                                                                                                                 ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc                   >  5:95913/1‑148 (MQ=255)
                                                                                                                                   agtgcGCCGGCCGGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGag                >  1:200769/5‑149 (MQ=255)
                                                                                                                                       cGCCGGCCGGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                          <  6:259511/119‑1 (MQ=255)
                                                                                                                                       cGCCGGCCGGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                          >  5:259511/1‑119 (MQ=255)
                                                                                                                                        gccggccgGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       >  3:36471/1‑105 (MQ=255)
                                                                                                                                        gccggccgGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       <  4:36471/105‑1 (MQ=255)
                                                                                                                                          cAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCTGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCgg         <  1:45560/149‑1 (MQ=255)
                                                                                                                                             ccTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGATCACATTGCCGGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGgtgt      <  7:216069/149‑1 (MQ=255)
                                                                                                                                              cTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGtt     <  7:73830/149‑1 (MQ=255)
                                                                                                                                                 ccggccggCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCa  <  1:67656/149‑1 (MQ=255)
                                                                                                                                                 ccggccggCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCa  <  5:70256/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCA  >  NC_002947/1245558‑1245851

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: