Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,704 C→A 7.0% 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,7040CA 7.0% 181.7 / 4.1 57P20Q (CCG→CAG) PP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base A (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CCAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCA  >  NC_002947/1245558‑1245851
                                                                                                                                                  |                                                                                                                                                   
ccAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccgg                                                                                                                                                   >  1:299043/1‑149 (MQ=255)
 cAATCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTgccggc                                                                                                                                                  <  7:1172/149‑1 (MQ=255)
    tCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGccggccgg                                                                                                                                               <  7:169090/149‑1 (MQ=255)
    tCCTACCGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGccggccgg                                                                                                                                               >  3:231543/1‑149 (MQ=255)
          cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAca                                                                                                                                         <  8:354235/149‑2 (MQ=255)
          cGGAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACt                                                                                                                                         <  4:231543/149‑1 (MQ=255)
            gAAAACCCCGGGGCGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCGGCCGGCAGGCTGGCGAAcagg                                                                                                                                       <  3:102178/149‑4 (MQ=255)
                         cGAGCACGGCTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGAcgcg                                                                                                                          <  2:269544/149‑1 (MQ=255)
                              ctatcTATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATg                                                                                                                      <  7:175590/144‑1 (MQ=255)
                                   tATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGaa                                                                                                                 >  7:240405/1‑148 (MQ=255)
                                   tATCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGaa                                                                                                                 >  3:333424/1‑148 (MQ=255)
                                     tCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCCCCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCaag                                                                                                                      >  8:175590/1‑139 (MQ=255)
                                     tCAAGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACAt                                                                                                              >  8:201672/1‑149 (MQ=255)
                                       aaGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACAtt                                                                                                             <  4:254824/148‑1 (MQ=255)
                                       aaGTTTTCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTg                                                                                                            >  7:59769/1‑149 (MQ=255)
                                             tCCGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTgg                                                                                                       <  6:219945/148‑1 (MQ=255)
                                               cGGACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAg                                                                                                    >  4:305377/1‑149 (MQ=255)
                                                 gACAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGt                                                                                                  <  3:305377/149‑1 (MQ=255)
                                                  aCAGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTc                                                                                                 <  7:201672/149‑1 (MQ=255)
                                                    aGCGTCCGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGTGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGa                                                                                               <  2:172606/149‑1 (MQ=255)
                                                         ccGATAGCCTGCTGGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTg                                                                                          <  7:353063/149‑1 (MQ=255)
                                                           gATAGACTGCTGGTAAGCCGTCAGAATTTTTGGAGCCGCCCTTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCTGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGATTGTGAc                                                                                        <  6:311843/149‑1 (MQ=255)
                                                                  tgctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTc                                                                                  >  5:54785/1‑148 (MQ=255)
                                                                    ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTc                                                                                  <  6:54785/146‑1 (MQ=255)
                                                                    ctgGTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGa                                                                                <  8:240405/148‑1 (MQ=255)
                                                                       gTAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAgtgt                                                                            >  3:328072/1‑149 (MQ=255)
                                                                        tAAGCCGTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtc                                                                           <  8:59769/149‑1 (MQ=255)
                                                                              gTCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGa                                                                     >  7:101970/1‑149 (MQ=255)
                                                                               tCAGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGag                                                                    >  4:49634/1‑149 (MQ=255)
                                                                                 aGAATTTGTGGAGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGagcc                                                                  <  4:333424/149‑1 (MQ=255)
                                                                                            aGCCGCCCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       <  3:49634/149‑1 (MQ=255)
                                                                                              ccgccCGTGCGCCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTTACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAgg                                                     <  8:246239/149‑1 (MQ=255)
                                                                                                   cGTGCGCCAGCGTTACCTCGCCCGGGTGACCCTGTGCGCCAGCCGGCCGGCCGGCGCACTGACGTTCCAGACGCGCAGGGAGAACATTGCCGTGAAGGTCGAGGGTGACCAGTTCGAGTGGCGTCTGAGCCAGCCGATCGACGGTTTCg                                                <  6:86807/149‑1 (MQ=255)
                                                                                                        gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAg                                                      >  8:3094/1‑138 (MQ=255)
                                                                                                        gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAg                                                      <  7:3094/138‑1 (MQ=255)
                                                                                                        gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGc                                           >  5:370644/1‑149 (MQ=255)
                                                                                                        gcCAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGc                                           >  5:141997/1‑149 (MQ=255)
                                                                                                         ccAGCGTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCg                                          >  1:265370/1‑149 (MQ=255)
                                                                                                              gTTACCTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACg                                     >  7:171873/1‑149 (MQ=255)
                                                                                                                   cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                >  5:794/1‑149 (MQ=255)
                                                                                                                   cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                <  2:265370/149‑1 (MQ=255)
                                                                                                                   cTAGCCCTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                <  6:143801/149‑1 (MQ=255)
                                                                                                                        ccTGTTGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTc                           <  8:171873/149‑1 (MQ=255)
                                                                                                                            ttGACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCggg                        <  8:101970/148‑1 (MQ=255)
                                                                                                                              gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCg                                                <  5:345034/122‑1 (MQ=255)
                                                                                                                              gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCg                                                >  6:345034/1‑122 (MQ=255)
                                                                                                                              gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                <  4:86669/138‑1 (MQ=255)
                                                                                                                              gACCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGt                                >  3:86669/1‑138 (MQ=255)
                                                                                                                               aCCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAg                                                                                     <  1:305544/84‑1 (MQ=255)
                                                                                                                               aCCCTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAg                                                                                     >  2:305544/1‑84 (MQ=255)
                                                                                                                                 ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGc                   >  4:202201/1‑148 (MQ=255)
                                                                                                                                 ccTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCg                  >  7:38590/1‑149 (MQ=255)
                                                                                                                                  cTGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGa                 <  6:794/149‑1 (MQ=255)
                                                                                                                                   tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtcgt                                                                         <  3:180471/92‑1 (MQ=255)
                                                                                                                                   tGTTCGCCAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTgtcgt                                                                         >  4:180471/1‑92 (MQ=255)
                                                                                                                                   agtgcGCCGGCCGGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGag                >  3:165199/5‑149 (MQ=255)
                                                                                                                                        gccggccgGCAGGCTGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGcc           >  6:340851/1‑149 (MQ=255)
                                                                                                                                         ccAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       <  3:239403/104‑1 (MQ=255)
                                                                                                                                         ccAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGaa                                                       >  4:239403/1‑104 (MQ=255)
                                                                                                                                         ccAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCg          >  4:129989/1‑149 (MQ=255)
                                                                                                                                          cAGCCTGCCGGCCGGCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCgg         >  6:77083/1‑149 (MQ=255)
                                                                                                                                                 ccggccggCGCACTGACTTTCCAGACGCGCATGGAGAACATTGCCTGGAAGGTCGAGGGTGACCAGTTCGAGTGTCGTCTGAGCCAGCCGATCGAAGGTTTCGGCAGCGGTACGTTCGTGCGTCGGGCAGGCGAGCAGCCGGTGTTTCa  <  8:38590/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: ‑