Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,280 G→T 10.6% A23D (GCC→GAC)  PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,2800GT10.6% 166.9 / 7.4 56A23D (GCC→GAC) PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base T (3/3);  total (30/26)
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

TGCTGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTC  >  NC_002947/3057134‑3057425
                                                                                                                                                   |                                                                                                                                                 
tGCTGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGc                                                                                                                                                  >  1:11761/1‑148 (MQ=255)
 gCTGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGc                                                                                                                                                  >  6:168447/1‑147 (MQ=255)
   tGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAgc                                                                                                                                              >  4:211342/1‑149 (MQ=255)
       gtgtGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGcc                                                                                                                                           <  2:11761/148‑1 (MQ=255)
              agcagcAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCGTTCTCGTTGGAGACCCCGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGcc                                                                                                                                   >  6:208737/1‑149 (MQ=255)
               gcagcaAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGccc                                                                                                                                  <  3:113472/149‑1 (MQ=255)
                   caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGa                                                                                                                                >  1:5537/1‑147 (MQ=255)
                   caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGAca                                                                                                                              <  5:168447/149‑1 (MQ=255)
                           gcgGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTg                                                                                                                      <  8:169049/149‑1 (MQ=255)
                                  cTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTa                                                                                                                 >  8:130603/1‑147 (MQ=255)
                                     gCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGCACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGc                                                                                                             <  1:76128/149‑1 (MQ=255)
                                       ccACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAgg                                                                                                          >  1:207648/1‑149 (MQ=255)
                                       ccACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAgg                                                                                                          <  3:202702/149‑1 (MQ=255)
                                         aCCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGc                                                                                                        >  5:60169/1‑149 (MQ=255)
                                          ccGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGcc                                                                                                       >  4:141414/1‑149 (MQ=255)
                                            ggCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCg                                                                                                     >  5:44213/1‑149 (MQ=255)
                                               aGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTg                                                                                                  >  3:266268/1‑149 (MQ=255)
                                                        gTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAgggg                                                                                         >  4:103998/1‑149 (MQ=255)
                                                        gTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAgggg                                                                                         <  2:207648/149‑1 (MQ=255)
                                                         tGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAggggg                                                                                        >  2:19055/1‑149 (MQ=255)
                                                           aCTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcg                                                                                      <  6:60169/149‑1 (MQ=255)
                                                             ttCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                    >  2:92656/1‑149 (MQ=255)
                                                             ttCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                    >  5:194758/1‑149 (MQ=255)
                                                                               gATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                   <  3:211342/148‑1 (MQ=255)
                                                                                   aggctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTAttgtt                                                              >  7:59846/1‑149 (MQ=255)
                                                                                      ctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTtgtg                                                           <  1:19055/149‑1 (MQ=255)
                                                                                      ctgaggctgTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTtgtg                                                           <  5:94627/149‑1 (MQ=255)
                                                                                               ttgtCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCgg                                                   <  7:130603/148‑1 (MQ=255)
                                                                                                  tCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGt                                                 >  3:35187/1‑147 (MQ=255)
                                                                                                     ttCTCGTTGGAGACACAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAggg                                            >  7:150997/1‑149 (MQ=255)
                                                                                                       ctcGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAggg                                            <  3:141414/147‑1 (MQ=255)
                                                                                                          gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTa                                         <  1:92656/147‑1 (MQ=255)
                                                                                                          gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGc                                       <  3:103998/149‑1 (MQ=255)
                                                                                                          gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTCTTGTGGTTGTCGGGTAAGGGGTAGc                                       >  6:269943/1‑149 (MQ=255)
                                                                                                               agaCCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTacac                                  <  4:35187/149‑1 (MQ=255)
                                                                                                                   ccAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCa                              >  1:85173/1‑149 (MQ=255)
                                                                                                                    cAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCAc                             <  3:118529/149‑1 (MQ=255)
                                                                                                                       gCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTc                          >  3:323457/1‑149 (MQ=255)
                                                                                                                           tGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGc                      <  5:208737/149‑1 (MQ=255)
                                                                                                                             gCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                    >  1:125485/1‑149 (MQ=255)
                                                                                                                              cAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTg                   <  6:44213/149‑1 (MQ=255)
                                                                                                                               aGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTgg                                                          <  2:261234/109‑1 (MQ=255)
                                                                                                                                 ccATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGt                                <  8:186039/133‑1 (MQ=255)
                                                                                                                                 ccATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGt                                >  7:186039/1‑133 (MQ=255)
                                                                                                                                   aTGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCa                >  3:77636/1‑147 (MQ=255)
                                                                                                                                      gtgtGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc           >  4:308571/1‑149 (MQ=255)
                                                                                                                                       tgtgACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc           >  6:134257/1‑148 (MQ=255)
                                                                                                                                        ccagcctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGcc           >  4:37918/6‑148 (MQ=38)
                                                                                                                                           gcctgctgTCGGGCGCCCTGGCGCAGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                    <  3:290165/134‑1 (MQ=255)
                                                                                                                                           gcctgctgTCGGGCGCCCTGGCGCAGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCtt                    >  4:290165/3‑136 (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)
                                                                                                                                                   |                                                                                                                                                 
TGCTGTCGTGTGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTC  >  NC_002947/3057134‑3057425

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: