Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,299 A→C 13.3% S17A (TCG→GCG) ‡ PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,2990AC13.3% 154.9 / 14.3 59S17A (TCG→GCG) ‡PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base C (4/4);  total (29/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 = 1.00e+00

CAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAGGGATGCAG  >  NC_002947/3057153‑3057446
                                                                                                                                                   |                                                                                                                                                   
caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGa                                                                                                                                                     >  1:5537/1‑147 (MQ=255)
caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGAca                                                                                                                                                   <  5:168447/149‑1 (MQ=255)
        gcgGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTg                                                                                                                                           <  8:169049/149‑1 (MQ=255)
               cTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTa                                                                                                                                      >  8:130603/1‑147 (MQ=255)
                  gCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGCACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGc                                                                                                                                  <  1:76128/149‑1 (MQ=255)
                    ccACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAgg                                                                                                                               <  3:202702/149‑1 (MQ=255)
                    ccACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAgg                                                                                                                               >  1:207648/1‑149 (MQ=255)
                      aCCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGc                                                                                                                             >  5:60169/1‑149 (MQ=255)
                       ccGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGcc                                                                                                                            >  4:141414/1‑149 (MQ=255)
                         ggCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCg                                                                                                                          >  5:44213/1‑149 (MQ=255)
                            aGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTg                                                                                                                       >  3:266268/1‑149 (MQ=255)
                                     gTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAgggg                                                                                                              >  4:103998/1‑149 (MQ=255)
                                     gTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAgggg                                                                                                              <  2:207648/149‑1 (MQ=255)
                                      tGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAggggg                                                                                                             >  2:19055/1‑149 (MQ=255)
                                        aCTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcg                                                                                                           <  6:60169/149‑1 (MQ=255)
                                          ttCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                                         >  2:92656/1‑149 (MQ=255)
                                          ttCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                                         >  5:194758/1‑149 (MQ=255)
                                                            gATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                                        <  3:211342/148‑1 (MQ=255)
                                                                aggctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTAttgtt                                                                                   >  7:59846/1‑149 (MQ=255)
                                                                   ctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTtgtg                                                                                <  5:94627/149‑1 (MQ=255)
                                                                   ctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTtgtg                                                                                <  1:19055/149‑1 (MQ=255)
                                                                            ttgtCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCgg                                                                        <  7:130603/148‑1 (MQ=255)
                                                                               tCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGt                                                                      >  3:35187/1‑147 (MQ=255)
                                                                                  ttCTCGTTGGAGACACAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAggg                                                                 >  7:150997/1‑149 (MQ=255)
                                                                                    ctcGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAggg                                                                 <  3:141414/147‑1 (MQ=255)
                                                                                       gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTa                                                              <  1:92656/147‑1 (MQ=255)
                                                                                       gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGc                                                            <  3:103998/149‑1 (MQ=255)
                                                                                       gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTCTTGTGGTTGTCGGGTAAGGGGTAGc                                                            >  6:269943/1‑149 (MQ=255)
                                                                                            agaCCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTacac                                                       <  4:35187/149‑1 (MQ=255)
                                                                                                ccAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCa                                                   >  1:85173/1‑149 (MQ=255)
                                                                                                 cAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCAc                                                  <  3:118529/149‑1 (MQ=255)
                                                                                                    gCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTc                                               >  3:323457/1‑149 (MQ=255)
                                                                                                        tGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGc                                           <  5:208737/149‑1 (MQ=255)
                                                                                                          gCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                                         >  1:125485/1‑149 (MQ=255)
                                                                                                           cAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTg                                        <  6:44213/149‑1 (MQ=255)
                                                                                                            aGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTgg                                                                               <  2:261234/109‑1 (MQ=255)
                                                                                                              ccATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGt                                                     >  7:186039/1‑133 (MQ=255)
                                                                                                              ccATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGt                                                     <  8:186039/133‑1 (MQ=255)
                                                                                                                aTGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCa                                     >  3:77636/1‑147 (MQ=255)
                                                                                                                   gtgtG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc                                >  4:308571/1‑149 (MQ=255)
                                                                                                                    tgtg‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc                                >  6:134257/1‑148 (MQ=255)
                                                                                                                     ccagcctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc                                >  4:37918/6‑148 (MQ=38)
                                                                                                                        gcctgctgTCGGGCGCCCTGGCGCAGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                                         >  4:290165/3‑136 (MQ=255)
                                                                                                                        gcctgctgTCGGGCGCCCTGGCGCAGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                                         <  3:290165/134‑1 (MQ=255)
                                                                                                                         aCTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAgg                            <  8:59846/148‑1 (MQ=255)
                                                                                                                          ctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTacac                                                       >  2:187571/1‑120 (MQ=38)
                                                                                                                          ctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTacac                                                       <  1:187571/120‑1 (MQ=38)
                                                                                                                          ctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCt                          <  3:37918/149‑1 (MQ=255)
                                                                                                                             ctgGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTc                       <  6:194758/149‑1 (MQ=255)
                                                                                                                              tgGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGATTGTCATAGCCCAGGCTGTc                       >  7:27473/1‑148 (MQ=255)
                                                                                                                                gCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTac                                                         <  2:158273/112‑1 (MQ=255)
                                                                                                                                gCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTac                                                         >  1:158273/1‑112 (MQ=255)
                                                                                                                                 cgggcgcCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCtt                                                                                         <  4:101563/76‑1 (MQ=37)
                                                                                                                                 cgggcgcCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCtt                                                                                         >  3:101563/4‑79 (MQ=37)
                                                                                                                                 ccAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAtt                   <  4:323457/149‑1 (MQ=255)
                                                                                                                                  cAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTGGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAttt                  >  7:108727/1‑149 (MQ=255)
                                                                                                                                      gcCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGtt              <  5:134257/149‑1 (MQ=255)
                                                                                                                                      gcCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGtt              <  4:77636/149‑1 (MQ=255)
                                                                                                                                                 cGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGGTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAGGGATGCa   <  8:27473/149‑1 (MQ=255)
                                                                                                                                                  gACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAGGGATGCAg  >  3:229592/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAGGGATGCAG  >  NC_002947/3057153‑3057446

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: ‑