Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,290 G→T 13.3% L20M (CTG→ATG) ‡ PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,2900GT13.3% 155.9 / 15.6 60L20M (CTG→ATG) ‡PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base T (4/4);  total (30/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 = 4.78e-01

AGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTT  >  NC_002947/3057148‑3057434
                                                                                                                                               |                                                                                                                                                
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                                                                                                                   >  1:207648/1‑149 (MQ=255)
                         ccACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAgg                                                                                                                   <  3:202702/149‑1 (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)
                                                                                                                                               |                                                                                                                                                
AGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTG‑ACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTT  >  NC_002947/3057148‑3057434

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