Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,587,817 A→C 6.3% L444R (CTG→CGG) ‡ parE ← DNA topoisomerase IV subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,587,8170AC 6.3% 198.9 / 4.6 63L444R (CTG→CGG) ‡parEDNA topoisomerase IV subunit B
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base C (2/2);  total (34/29)
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.98e-01

GCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGA  >  NC_002947/5587670‑5587964
                                                                                                                                                   |                                                                                                                                                   
gCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAg                                                                                                                                                    >  4:114923/1‑149 (MQ=255)
          aGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATg                                                                                                                                          >  3:181809/1‑149 (MQ=255)
          aGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATg                                                                                                                                          >  1:322087/1‑149 (MQ=255)
          aGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCCTGCCGCTGCGCGGCAAGAtc                                                                                                                                          <  5:249976/149‑2 (MQ=255)
                  ccATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGa                                                                                                                                   <  8:207090/148‑1 (MQ=255)
                   cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                                                                                 >  6:282896/1‑149 (MQ=255)
                   cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                                                                                 >  3:6700/1‑149 (MQ=255)
                   cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                                                                                 >  1:73550/1‑149 (MQ=255)
                   cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                                                                                 >  3:91880/1‑149 (MQ=255)
                    aTAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                                                                                >  4:305377/1‑149 (MQ=255)
                    aTAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                                                                                <  6:153877/149‑1 (MQ=255)
                     tAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTc                                                                                                                               <  3:305377/149‑1 (MQ=255)
                       gcgcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                                                                              <  3:298001/148‑1 (MQ=255)
                         gcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTg                                                                                                                           <  4:181809/149‑1 (MQ=255)
                          cAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGt                                                                                                                          >  7:180196/1‑149 (MQ=255)
                                       tCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCtt                                                                                                             <  2:216873/149‑1 (MQ=255)
                                         ggTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTgg                                                                                                           <  5:282896/149‑1 (MQ=255)
                                           ttGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGcc                                                                                                         >  4:18590/1‑149 (MQ=255)
                                                      tcgtcgCCTATGGCCACGGGAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGGATGATGGCCTGGAATTCCTTGTCCCGCTCCTGCTTGGCCGAACCACCgg                                                                                               <  8:111473/148‑1 (MQ=255)
                                                         aCGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGa                                                                                           <  3:114923/149‑1 (MQ=255)
                                                         aCGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCCTGCCGCTGCGCGGCAAGAtc                                                                                                                                          <  3:121150/102‑2 (MQ=255)
                                                         aCGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCCTGCCGCTGCGCGGCAAGAtc                                                                                                                                          >  4:121150/1‑101 (MQ=255)
                                                            ccAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAac                                                                                                     >  2:307155/1‑136 (MQ=255)
                                                            ccAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAac                                                                                                     <  1:307155/136‑1 (MQ=255)
                                                                    cACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAcc                                                                                     <  6:53419/144‑1 (MQ=255)
                                                                    cACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAcc                                                                                     >  5:53419/1‑144 (MQ=255)
                                                                    cACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcga                                                                                <  3:18590/149‑1 (MQ=255)
                                                                       ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgac                                                                               <  2:336341/147‑1 (MQ=255)
                                                                       ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgac                                                                               >  1:336341/1‑147 (MQ=255)
                                                                        gCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAg                                                                            >  1:57655/1‑149 (MQ=255)
                                                                               tgtgCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTc                                                                                                                         >  5:106768/1‑97 (MQ=255)
                                                                               tgtgCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTc                                                                                                                         <  6:106768/97‑1 (MQ=255)
                                                                                 tgCAGTTCCTGTGTGGCCAGGACTTCGCCCCCGTCCTCTTCACTGGTGTTGAGGATCTTGCCGCGCAGCTGCTGGATGGCCTGGATTTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCAGAGTGAGGCTCGACGAGGAACAAGTc                                                                   <  4:91880/149‑1 (MQ=255)
                                                                                 tgCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTc                                                                   >  4:326831/1‑149 (MQ=255)
                                                                                              tggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCggg                                                      >  2:212585/1‑149 (MQ=255)
                                                                                               ggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGt                                                     >  2:106502/1‑149 (MQ=255)
                                                                                                gcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCggg                                                      <  7:225846/147‑1 (MQ=255)
                                                                                                  cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACCACTCGGCACGCATCGGGTCCt                                                  >  6:346692/1‑149 (MQ=255)
                                                                                                  cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCAGCACGCATCGGGTCCt                                                  >  8:249961/1‑149 (MQ=255)
                                                                                                  cAGGACTTCGCCCCCGCCCACTTCCCAGGTGTTGAGGATCTTGCCGCCCAGCGCCAGGATGGCCTCGAATTCCTTGTCCCGCGCCTCCTTGGCCGAACCACCGGCCGAGTCACCCTCGACCAGGAACAACACGGCACCCATCgcgaccg                                                  >  3:23346/1‑143 (MQ=255)
                                                                                                   aGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGTGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg                                                                           >  6:73381/1‑123 (MQ=255)
                                                                                                   aGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGTGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg                                                                           <  5:73381/123‑1 (MQ=255)
                                                                                                   aGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                 >  6:59447/1‑149 (MQ=255)
                                                                                                   aGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                 >  8:24221/1‑149 (MQ=255)
                                                                                                   aGGACTTCGCCCCCGACCACTTCCGAGGGGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAGCCTCGACGAGGAACAACTCGGCACGCATCGGGGCCTg                                                 <  4:23346/149‑1 (MQ=255)
                                                                                                        ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTg                                            >  2:254463/1‑149 (MQ=255)
                                                                                                        ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTg                                            >  3:291145/1‑149 (MQ=255)
                                                                                                             ccccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACaa                                       >  5:241895/1‑149 (MQ=255)
                                                                                                               ccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATc                                     <  2:322087/149‑1 (MQ=255)
                                                                                                                 ggccACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCgg                                   <  8:160090/147‑1 (MQ=255)
                                                                                                                    cACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGc                                  >  2:121863/1‑147 (MQ=255)
                                                                                                                    cACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCa                                <  2:73550/149‑1 (MQ=255)
                                                                                                                              gtgtTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCggg                      <  1:212585/149‑1 (MQ=255)
                                                                                                                                   ggccaTCCTGCCGCTGCGCGGCAAGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATg                 >  5:55671/5‑149 (MQ=255)
                                                                                                                                   gAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGTGCCTGCtt                                                                                                             >  8:8914/1‑57 (MQ=255)
                                                                                                                                   gAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGTGCCTGCtt                                                                                                             <  7:8914/57‑1 (MQ=255)
                                                                                                                                      gATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGc                          <  5:110446/137‑1 (MQ=255)
                                                                                                                                      gATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGc                          >  6:110446/1‑137 (MQ=255)
                                                                                                                                              cgcgCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTgg      >  7:87032/1‑149 (MQ=255)
                                                                                                                                               gcgcAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTggg     >  3:245749/1‑149 (MQ=255)
                                                                                                                                                cgcAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGt    <  4:6700/149‑1 (MQ=255)
                                                                                                                                                  cAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGa  <  2:57655/149‑1 (MQ=255)
                                                                                                                                                  cAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCAGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGa  <  7:249961/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGA  >  NC_002947/5587670‑5587964

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