Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,282 C→G 10.7% L22L (CTG→CTC) ‡ PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,2820CG10.7% 162.2 / 8.6 56L22L (CTG→CTC) ‡PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base G (3/3);  total (29/27)
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.83e-01

TGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTT  >  NC_002947/3057137‑3057430
                                                                                                                                                  |                                                                                                                                                    
tGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAgc                                                                                                                                                   >  4:211342/1‑149 (MQ=255)
    gtgtGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGcc                                                                                                                                                <  2:11761/148‑1 (MQ=255)
           agcagcAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCGTTCTCGTTGGAGACCCCGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGcc                                                                                                                                        >  6:208737/1‑149 (MQ=255)
            gcagcaAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGccc                                                                                                                                       <  3:113472/149‑1 (MQ=255)
                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                                                                <  1:19055/149‑1 (MQ=255)
                                                                                   ctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTtgtg                                                                <  5:94627/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)
                                                                                                                                                 ccAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAtt   <  4:323457/149‑1 (MQ=255)
                                                                                                                                                  cAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTGGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAttt  >  7:108727/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTT  >  NC_002947/3057137‑3057430

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