Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,154,350 T→A 6.7% E141V (GAA→GTA)  PP_1909 ← septum formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,154,3500TA 6.7% 170.5 / 2.9 60E141V (GAA→GTA) PP_1909septum formation protein
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ATCAGCGGCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCT  >  NC_002947/2154203‑2154494
                                                                                                                                                   |                                                                                                                                                 
aTCAGCGGCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGtt                                                                                                                                                   >  3:322299/1‑148 (MQ=255)
 tCAGCGGCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCa                                                                                                                                                 >  6:360598/1‑149 (MQ=255)
     cGGCAGCCAGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCa                                                                                                                                             <  6:290197/149‑1 (MQ=255)
       gCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAgt                                                                                                                                           <  1:244254/149‑1 (MQ=255)
       gCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAgt                                                                                                                                           <  5:262214/149‑1 (MQ=255)
           cccGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCa                                                                                                                                       <  2:171962/149‑1 (MQ=255)
            ccGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCAc                                                                                                                                      >  2:192299/1‑149 (MQ=255)
                  aGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTaa                                                                                                                                 <  5:360598/148‑1 (MQ=255)
                   ggCTGGTCGCGGCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTaaaa                                                                                                                               <  7:14421/149‑1 (MQ=255)
                        gtcgcgtcgcACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTATGCTGCCACGTAACGATCGACATGAACCCTGCTCAGTTCACGCAGTGTCAc                                                                                                                                      >  4:297397/1‑137 (MQ=255)
                        gtagcgtcgcACCCGTGCGTGCTCTGGAACAAGCTCACCCCCATCCCCTCCGCATTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCAc                                                                                                                                      <  3:297397/134‑1 (MQ=255)
                         tcgcgtcgcACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑Cg                                                                                                                        <  4:29291/149‑1 (MQ=255)
                           gcgtcgcACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCAc                                                                                                                                                >  8:201865/1‑124 (MQ=255)
                           gcgtcgcACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCAc                                                                                                                                                <  7:201865/124‑1 (MQ=255)
                                            gcTCTGGAACAAGCTCACCCCCACCCCCTCCGCCTTGAAACAACCAGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCAAGACCAGCTCCCTGCTCAGTTAACGCAGTGTCAAGGTAAAATGCA‑CGCAATACAACAGGCATTTcc                                                                                                     >  3:345892/1‑149 (MQ=255)
                                                  gAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCATCGCAATCCACCTGGCATTGCCCGGTg                                                                                                <  1:187056/149‑1 (MQ=255)
                                                  gAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTg                                                                                                             <  7:185667/135‑1 (MQ=255)
                                                  gAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTg                                                                                                <  4:322299/148‑1 (MQ=255)
                                                  gAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGATAACGCAGTGTCACGGTAAAAGGCA‑CGCAATACACCTg                                                                                                             >  8:185667/1‑135 (MQ=255)
                                                        gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTTCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGt                                                                                         <  4:345892/149‑1 (MQ=255)
                                                        gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTgg                                                                                               <  6:89724/143‑1 (MQ=255)
                                                        gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGt                                                                                         >  2:26644/1‑149 (MQ=255)
                                                        gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCAAGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTgg                                                                                               >  5:89724/1‑143 (MQ=255)
                                                              ccccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                     >  7:319184/1‑147 (MQ=255)
                                                              ccccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                     >  4:290720/1‑147 (MQ=255)
                                                               cccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagcag                                                                                  >  5:53567/1‑149 (MQ=255)
                                                                ccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                 >  3:7844/1‑149 (MQ=255)
                                                                ccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                 >  7:76699/1‑149 (MQ=255)
                                                                ccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTACACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                 >  4:76814/1‑149 (MQ=255)
                                                                ccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTACACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                 >  4:71785/1‑149 (MQ=255)
                                                                  aGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                               <  8:319184/149‑1 (MQ=255)
                                                                      ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTCCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGc                                                                           >  3:99564/1‑149 (MQ=255)
                                                                       cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTgg                                                                                                            >  7:171368/1‑115 (MQ=255)
                                                                       cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTgg                                                                                                            <  8:171368/115‑1 (MQ=255)
                                                                                 aaaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGa                                                                  >  1:46803/1‑147 (MQ=255)
                                                                                 aaaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTACACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaa                                                                 <  3:71785/148‑1 (MQ=255)
                                                                                 aaaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTACACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaa                                                                 <  3:76814/148‑1 (MQ=255)
                                                                                  aaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAg                                                               <  2:46803/149‑1 (MQ=255)
                                                                                        cTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCAcc                                                         >  3:42306/1‑149 (MQ=255)
                                                                                              aaTCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCtgctgc                                                   <  1:26644/149‑1 (MQ=255)
                                                                                                tCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGccc                                                 >  3:242750/1‑149 (MQ=255)
                                                                                                  cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACTGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCg                                                <  2:305696/148‑1 (MQ=255)
                                                                                                  cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGc                                               <  1:192299/149‑1 (MQ=255)
                                                                                                   aGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                              >  2:69245/1‑149 (MQ=255)
                                                                                                       gCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCAc                                          <  4:42306/149‑1 (MQ=255)
                                                                                                         tGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                        >  3:286818/1‑149 (MQ=255)
                                                                                                             ctgctgCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCag                                     <  3:290720/148‑1 (MQ=255)
                                                                                                               gctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCagtag                                  >  4:150140/1‑149 (MQ=255)
                                                                                                                   ccACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGc                              >  3:59837/1‑149 (MQ=255)
                                                                                                                     aCGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                            <  4:59837/149‑1 (MQ=255)
                                                                                                                          aCGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                        <  8:76699/148‑1 (MQ=255)
                                                                                                                          aCGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                        <  4:99564/148‑1 (MQ=255)
                                                                                                                            gCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAgcg                     >  3:168178/1‑149 (MQ=255)
                                                                                                                               cGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCt                             >  6:199575/1‑138 (MQ=255)
                                                                                                                               cGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCt                             <  5:199575/138‑1 (MQ=255)
                                                                                                                                         cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTAt        >  7:64586/1‑149 (MQ=255)
                                                                                                                                             ctcAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCg     <  6:53567/148‑1 (MQ=255)
                                                                                                                                              tcAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                             >  2:339478/1‑59 (MQ=255)
                                                                                                                                              tcAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                             <  1:339478/59‑1 (MQ=255)
                                                                                                                                              tcAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGc   >  1:304307/1‑149 (MQ=255)
                                                                                                                                                aGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCt  <  4:168178/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCAGCGGCAGCCCGATCAGGCTGGTCGCGTCGCACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCA‑CGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCT  >  NC_002947/2154203‑2154494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑