Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,587,765 G→A 5.5% A461A (GCC→GCT)  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,7650GA 5.5% 282.3 / 2.8 91A461A (GCC→GCT) parEDNA topoisomerase IV subunit B
Reads supporting (aligned to +/- strand):  ref base G (42/44);  new base A (3/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GAACAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGT  >  NC_002947/5587618‑5587913
                                                                                                                                                      |                                                                                                                                                        
gaacAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCctggctgg                                                                                                                                                         <  3:371066/149‑1 (MQ=255)
  acAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCtggctggc                                                                                                                                                        >  2:278117/1‑148 (MQ=255)
   cAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCtggctggc                                                                                                                                                        >  3:263793/1‑147 (MQ=255)
   cAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCa                                                                                                                                                      <  7:60934/149‑1 (MQ=255)
    aGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAg                                                                                                                                                     >  3:367225/1‑149 (MQ=255)
    aGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAg                                                                                                                                                     >  3:361162/1‑149 (MQ=255)
      cgcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGa                                                                                                                                                   >  6:420541/1‑149 (MQ=255)
      cgcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGa                                                                                                                                                   <  6:181118/149‑1 (MQ=255)
       gcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCCAGCCAGGaa                                                                                                                                                  <  7:81485/149‑2 (MQ=255)
       gcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACGTCCTGGCCAGCCAGGaa                                                                                                                                                  <  3:133559/149‑2 (MQ=255)
         acacaGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACtt                                                                                                                                                >  1:155267/1‑149 (MQ=255)
           acaGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCg                                                                                                                                              <  7:225772/149‑1 (MQ=255)
               cAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGcccc                                                                                                                                          >  2:152688/1‑149 (MQ=255)
               cAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCAGCGAGGATGCAGATCTTACCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGcccc                                                                                                                                          >  4:145839/1‑149 (MQ=255)
                           gtgCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGTTCGGTTGCACCCG‑GA‑TCTACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCTTCCA‑‑CTTcc                                                                                                                             <  6:421584/148‑1 (MQ=255)
                            tgCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTcc‑‑c                                                                                                                          <  4:145047/148‑1 (MQ=255)
                              cAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAgg                                                                                                                       >  7:97332/1‑149 (MQ=255)
                              cAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAg                                                                                                                        <  1:152688/148‑1 (MQ=255)
                                  ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGATAGAGGACGCCAATGGCCAGCGGCAATGTTGTGCACTTACTGGCTGTCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGgt                                                                                                                      >  1:353818/1‑149 (MQ=255)
                                  ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGgtgt                                                                                                                    <  1:269587/148‑1 (MQ=255)
                                  ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                                   <  1:410984/149‑1 (MQ=255)
                                   cATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCa                                                                                                                                    <  5:31682/135‑1 (MQ=255)
                                   cATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCa                                                                                                                                    >  6:31682/1‑135 (MQ=255)
                                    aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                                   >  8:251526/1‑147 (MQ=255)
                                    aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                                   >  2:318600/1‑147 (MQ=255)
                                    aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGa                                                                                                                 <  5:420541/149‑1 (MQ=255)
                                    aTCGGAGTCAGCGTCGGCGAGGAGGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGa                                                                                                                 <  8:355321/149‑1 (MQ=255)
                                     tCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAg                                                                                                                >  6:27816/1‑149 (MQ=255)
                                         aGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATc                                                                                                            <  6:278123/149‑1 (MQ=255)
                                           tCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCtt                                                                                                          >  4:101753/1‑149 (MQ=255)
                                           tCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCtt                                                                                                          >  5:168783/1‑149 (MQ=255)
                                              gCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGcc                                                                                                       >  2:141532/1‑149 (MQ=255)
                                                  cGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCcgcg                                                                                                    >  6:182536/1‑148 (MQ=255)
                                                    gCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAg                                                                                                 >  5:217462/1‑149 (MQ=255)
                                                      gAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCg                                                                                               >  3:211701/1‑149 (MQ=255)
                                                          aTGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAg                                                                                           <  7:219374/149‑1 (MQ=255)
                                                           tGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAgg                                                                                          <  6:217462/149‑1 (MQ=255)
                                                            gCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGa                                                                                         >  5:77339/1‑149 (MQ=255)
                                                             cAGATCTTACCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGa                                                                                         <  3:145839/148‑1 (MQ=255)
                                                                aTCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGc                                                                                     >  6:430153/1‑149 (MQ=255)
                                                                  cTTGCCATAGCGCAGCTGCGCAAGGTCGGTCGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   >  7:35248/1‑149 (MQ=255)
                                                                     gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   <  7:251526/146‑1 (MQ=255)
                                                                     gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                                 >  7:348423/1‑148 (MQ=255)
                                                                      ccATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGa                                                                                <  6:168783/148‑1 (MQ=255)
                                                                       cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   >  1:369631/1‑144 (MQ=255)
                                                                       cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   <  2:369631/144‑1 (MQ=255)
                                                                       cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                              <  5:430153/149‑1 (MQ=255)
                                                                       cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                              >  6:327473/1‑149 (MQ=255)
                                                                        aTAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                             >  3:50339/1‑149 (MQ=255)
                                                                        aTAGCGCAGCTGCGCAAGGTCGGTCGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                             <  8:35248/149‑1 (MQ=255)
                                                                           gcgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                           <  4:361162/148‑1 (MQ=255)
                                                                           gcgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                           <  4:367225/148‑1 (MQ=255)
                                                                            cgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCtt                                                                         <  5:75031/149‑1 (MQ=255)
                                                                             gcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑Ctt                                                                                                                               >  7:158122/1‑96 (MQ=255)
                                                                             gcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑Ctt                                                                                                                               <  8:158122/96‑1 (MQ=255)
                                                                              cAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGCCGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCACCGGCAGGATGGCCTGGAATTCCTTGt                                                                       >  3:284801/1‑149 (MQ=255)
                                                                                       aaGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCt                                                              >  2:274058/1‑149 (MQ=255)
                                                                                       aaGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCt                                                              >  2:18604/1‑149 (MQ=255)
                                                                                         ggTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGc                                                            <  4:50339/149‑1 (MQ=255)
                                                                                          gTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCt                                                           <  4:211701/149‑1 (MQ=255)
                                                                                              gTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                                       >  7:255809/1‑149 (MQ=255)
                                                                                              gTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                                       <  2:353818/149‑1 (MQ=255)
                                                                                               ttGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGcc                                                      <  3:27280/149‑1 (MQ=255)
                                                                                                 gCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                                 <  2:233290/120‑1 (MQ=255)
                                                                                                 gCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                                 >  1:233290/1‑120 (MQ=255)
                                                                                                  cACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   >  6:146728/1‑117 (MQ=255)
                                                                                                  cACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                   <  5:146728/117‑1 (MQ=255)
                                                                                                    cccG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAac                                                  >  4:57213/1‑148 (MQ=255)
                                                                                                     ccG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAacca                                                >  4:111207/1‑149 (MQ=255)
                                                                                                      cG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccac                                               <  4:152497/149‑1 (MQ=255)
                                                                                                       g‑gA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccacc                                              <  1:17097/149‑1 (MQ=255)
                                                                                                          a‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCgg                                            <  3:55327/149‑1 (MQ=255)
                                                                                                            tCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGc                                           <  1:278117/149‑1 (MQ=255)
                                                                                                            tCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGc                                           <  8:255809/149‑1 (MQ=255)
                                                                                                              gACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                          <  3:111207/148‑1 (MQ=255)
                                                                                                              gACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                          <  4:263793/148‑1 (MQ=255)
                                                                                                               aCGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                          >  4:35186/1‑147 (MQ=255)
                                                                                                                 gCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCCGGCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                          >  1:348781/1‑147 (MQ=255)
                                                                                                                   cAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCa                                    >  6:422151/1‑149 (MQ=255)
                                                                                                                   cAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCa                                    >  8:252663/1‑149 (MQ=255)
                                                                                                                      tGGCCA‑GGTCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAccc                                 <  3:56889/149‑1 (MQ=255)
                                                                                                                         ccA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcg                              <  3:35186/149‑1 (MQ=255)
                                                                                                                          cA‑CGGCAATTTTTTGCAATTCCTGGCTGGCCAGGAATTCGCCCCCGTCCACTCTTCC‑‑CCGGACTTCTCCCATTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTCGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTc                               <  1:141532/149‑1 (MQ=255)
                                                                                                                              ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgacga                          <  2:155267/149‑1 (MQ=255)
                                                                                                                                cAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg                        <  5:422151/149‑1 (MQ=255)
                                                                                                                                 aaTGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg                        >  4:94219/1‑148 (MQ=255)
                                                                                                                                   tgttgtGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaac                     <  2:348781/149‑1 (MQ=255)
                                                                                                                                   tgttgtGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaa                      >  1:174688/1‑148 (MQ=255)
                                                                                                                                         gCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg               >  8:264801/1‑149 (MQ=255)
                                                                                                                                          cACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg               <  7:264801/148‑1 (MQ=255)
                                                                                                                                            agtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa            >  1:153772/3‑149 (MQ=255)
                                                                                                                                            agtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa            >  1:218697/3‑149 (MQ=255)
                                                                                                                                             ttCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAc           >  8:55211/1‑149 (MQ=255)
                                                                                                                                             gtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGTGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAc           >  6:203607/2‑149 (MQ=255)
                                                                                                                                                      ggcCAGGACTTCGCCCCCCGCCA‑‑CTTCC‑‑CCGGTGGTGAGGATCTGGCCGCGCAGCGGCAGGATGGCCGGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGt  <  4:284801/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
GAACAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGT  >  NC_002947/5587618‑5587913

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