Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,399,900 T→A 5.5% L126Q (CTG→CAG)  oprL → peptidoglycan‑associated lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,399,9000TA 5.5% 227.8 / 2.7 72L126Q (CTG→CAG) oprLpeptidoglycan‑associated lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base A (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATGCGCGCTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCC  >  NC_002947/1399758‑1400046
                                                                                                                                              |                                                                                                                                                  
aTGCGCGCTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTgcagaa                                                                                                                                              >  6:254932/1‑144 (MQ=255)
 tGCGCGCTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTg                                                                                                                                             <  3:119886/149‑1 (MQ=255)
 tGCGCGCTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTg                                                                                                                                             <  2:208465/149‑1 (MQ=255)
     cgcTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAgg                                                                                                                                         <  4:88700/149‑1 (MQ=255)
            gACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGttt                                                                                                                                   <  5:335305/148‑1 (MQ=255)
            gACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTc                                                                                                                                  >  5:120233/1‑149 (MQ=255)
                ttCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCt                                                                                                                              <  1:60608/149‑1 (MQ=255)
                ttCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCt                                                                                                                              >  2:189170/1‑149 (MQ=255)
                   aCGCCAAGGACCTGAAAGCCTACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACCAGATGGCGCTGGGTGAGCGGCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGc                                                                                                                            <  5:121901/148‑1 (MQ=255)
                     gCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCt                                                                                                                           >  8:172093/1‑147 (MQ=255)
                        aaGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCt                                                                                                                      <  1:189170/149‑1 (MQ=255)
                           gACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTgg                                                                                                                    >  4:44322/1‑148 (MQ=255)
                            aCCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGaa                                                                                                                  >  5:62569/1‑149 (MQ=255)
                               tGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCAgg                                                                                                                                                  <  6:20836/114‑1 (MQ=255)
                               tGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCAgg                                                                                                                                                  <  2:201645/114‑1 (MQ=255)
                               tGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCAgg                                                                                                                                                  >  1:201645/1‑114 (MQ=255)
                               tGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCAgg                                                                                                                                                  >  5:20836/1‑114 (MQ=255)
                                gAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTgg                                                                                                              >  3:62633/1‑149 (MQ=255)
                                    gCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTc                                                                                                          <  3:44322/149‑1 (MQ=255)
                                      cAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTcc                                                                                                         >  5:45726/1‑148 (MQ=255)
                                       aaCGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTa                                                                                                       >  5:162579/1‑149 (MQ=255)
                                          ggCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACgg                                                                                                    <  4:62633/149‑1 (MQ=255)
                                             aaCCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCg                                                                                                  >  5:186081/1‑148 (MQ=255)
                                              aCCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGaa                                                                                                <  4:159076/149‑1 (MQ=255)
                                                cGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTGAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAaga                                                                                              <  1:116690/149‑1 (MQ=255)
                                                  tgttgtTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGc                                                                                            <  7:172093/149‑1 (MQ=255)
                                                  tgttgtTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGc                                                                                            <  6:120233/149‑1 (MQ=255)
                                                          tGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGtt                                                                                    >  5:111795/1‑149 (MQ=255)
                                                          tGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTCTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGtt                                                                                    >  1:343206/1‑149 (MQ=255)
                                                            gAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGATGGAACTGGTTTCCTACGGCGAAGAGCGTCAAGTTg                                                                                   >  7:122538/1‑148 (MQ=255)
                                                                    cacCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGc                                                                           >  3:94523/1‑148 (MQ=255)
                                                                       cgacgaGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAacg                                                                       >  4:241436/1‑149 (MQ=255)
                                                                         acgaGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGc                                                                                  >  3:320014/1‑136 (MQ=255)
                                                                         acgaGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGc                                                                                  <  4:320014/136‑1 (MQ=255)
                                                                         acgaGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAacgac                                                                     <  6:62569/149‑1 (MQ=255)
                                                                           gaGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAAcgacga                                                                   <  6:111795/149‑1 (MQ=255)
                                                                            aGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTCTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAg                                                                  <  2:343206/149‑1 (MQ=255)
                                                                                   cACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCt                                                           <  4:94523/149‑1 (MQ=255)
                                                                                      ccGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTggg                                                        >  3:327692/1‑149 (MQ=255)
                                                                                        gcgAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAac                                                                        <  6:341761/131‑1 (MQ=255)
                                                                                          gAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGaa                                                                                                >  4:101901/1‑105 (MQ=255)
                                                                                          gAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGaa                                                                                                <  3:101901/105‑1 (MQ=255)
                                                                                          gAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAAcgacga                                                                   >  4:111628/1‑134 (MQ=255)
                                                                                             tACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTggg                                                        <  2:265330/142‑1 (MQ=255)
                                                                                             tACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTggg                                                        >  1:265330/1‑142 (MQ=255)
                                                                                                aacaTGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACGTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCg                                              >  1:245262/1‑149 (MQ=255)
                                                                                                aacaTGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAg                                                                  >  1:137564/1‑129 (MQ=255)
                                                                                                aacaTGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCg                                              >  2:229882/1‑149 (MQ=255)
                                                                                                  caTGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTc                                            >  1:328587/1‑149 (MQ=255)
                                                                                                       ctctGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTa                                       <  6:45726/149‑1 (MQ=255)
                                                                                                        tctGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAg                                                   <  8:61068/136‑1 (MQ=255)
                                                                                                        tctGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAg                                                   >  7:61068/1‑136 (MQ=255)
                                                                                                        tctGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTcgcg                                         >  6:53233/1‑146 (MQ=255)
                                                                                                              tGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGc                                <  6:162579/149‑1 (MQ=255)
                                                                                                                  cgtcgtGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACt                                  <  1:222959/143‑1 (MQ=255)
                                                                                                                  cgtcgtGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACt                                  >  2:222959/1‑143 (MQ=255)
                                                                                                                    tcgtGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCgtaagt                          <  5:53233/149‑1 (MQ=255)
                                                                                                                    tcgtGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCgtaagt                          <  4:327692/149‑1 (MQ=255)
                                                                                                                      gtGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCgtaagt                          >  1:15577/1‑147 (MQ=255)
                                                                                                                         cGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTc                                                     <  4:5542/117‑1 (MQ=255)
                                                                                                                         cGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTc                                                     >  3:5542/1‑117 (MQ=255)
                                                                                                                           aaaGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTcgc                                          <  3:243081/126‑1 (MQ=255)
                                                                                                                           aaaGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTcgc                                          >  4:243081/1‑126 (MQ=255)
                                                                                                                             aGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTa                 <  2:328587/149‑1 (MQ=255)
                                                                                                                              gCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTAt                <  3:241436/149‑1 (MQ=255)
                                                                                                                                 tctCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGc              <  7:325201/146‑1 (MQ=255)
                                                                                                                                  ttCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCgt                              <  2:291719/131‑1 (MQ=255)
                                                                                                                                  ttCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCgt                              >  1:291719/1‑131 (MQ=255)
                                                                                                                                  ttCAGCGCTACCTGGTTCTGCAGGGCGTTTCCACTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCg                                              >  7:99821/1‑115 (MQ=255)
                                                                                                                                     aGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGc              >  8:325201/1‑144 (MQ=255)
                                                                                                                                       cgcTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTAtgtg       >  4:160678/1‑149 (MQ=255)
                                                                                                                                            ccTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTTGTGTCCTACGGCGATGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTgccgcc  <  8:122538/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ATGCGCGCTCTGGACGTTCACGCCAAGGACCTGAAAGCCAACGGCAACCGTGTTGTTCTGGAAGGTAACACCGACGAGCGCGGCACCCGCGAGTACAACATGGCTCTGGGTGAGCGTCGTGCGAAAGCCGTTCAGCGCTACCTGGTTCTGCAGGGCGTTTCCCCTGCTCAGCTGGAACTGGTTTCCTACGGCGAAGAGCGTCCAGTTGCTACCGGCAACGACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCC  >  NC_002947/1399758‑1400046

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