Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,417,745 T→A 10.6% E99V (GAA→GTA) ‡ PP_3030 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,417,7450TA10.6% 164.5 / 7.5 64E99V (GAA→GTA) ‡PP_3030hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (5/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGGGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGA  >  NC_002947/3417597‑3417889
                                                                                                                                                    |                                                                                                                                                  
aGATGGGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCtt                                                                                                                                                    <  1:344965/149‑1 (MQ=255)
      gTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCgcag                                                                                                                                               >  5:123016/1‑148 (MQ=255)
       tGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCgcagca                                                                                                                                             <  3:246461/149‑1 (MQ=255)
       tGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGATGATGGGGGGAATAATGCCTGCTGCTTCgcag                                                                                                                                               >  6:127072/1‑147 (MQ=255)
             taCCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACt                                                                                                                                       <  8:17623/149‑1 (MQ=255)
               ccccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACttt                                                                                                                                     >  4:182296/1‑149 (MQ=255)
                cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtttt                                                                                                                                    >  8:11688/1‑149 (MQ=255)
                  cACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTg                                                                                                                                  >  2:106568/1‑149 (MQ=255)
                  cACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTg                                                                                                                                  >  7:319977/1‑149 (MQ=255)
                        aGGCTGCGCGTGCGGGCGACTGTTAGCGATTAGCGTTGCAGCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGtttttt                                                                                                                            <  3:273363/149‑1 (MQ=255)
                          gCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGATACTTTTTGTTTTTTgc                                                                                                                          >  7:167279/1‑149 (MQ=255)
                             gcgcGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGATTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGCGGGAATAATGCATGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTgcgat                                                                                                                       >  6:209558/1‑147 (MQ=255)
                                    cGGGCGACTTTTAGCGATTAGCGTAGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGCCTGCTGATGGGGGCAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGacaa                                                                                                                <  5:127072/149‑1 (MQ=255)
                                     gggCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGacaac                                                                                                               <  8:319977/149‑1 (MQ=255)
                                        cGACTTTTTGGGATTAGGGTTGCATCAGCGGGGGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACg                                                                                                            <  8:167279/149‑1 (MQ=255)
                                          aCTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGa                                                                                                          <  7:11688/149‑1 (MQ=255)
                                            ttttAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACt                                                                                                        >  4:57687/1‑149 (MQ=255)
                                               tAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTcg                                                                                                     >  4:331494/1‑149 (MQ=255)
                                               tAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTc                                                                                                      <  3:172785/148‑1 (MQ=255)
                                                  cGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTcgcg                                                                                                   >  1:369485/1‑148 (MQ=255)
                                                    aTTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAg                                                                                                >  8:306279/1‑149 (MQ=255)
                                                       aGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCg                                                                                              <  2:369485/148‑1 (MQ=255)
                                                        gCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCg                                                                                              >  7:120413/1‑147 (MQ=255)
                                                          gTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGc                                                                                          >  2:62270/1‑149 (MQ=255)
                                                             gCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCTCCGCGAAgtcc                                                                                        <  4:183236/148‑4 (MQ=255)
                                                              cATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCgg                                                                                        >  8:146212/1‑147 (MQ=255)
                                                               aTCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAg                                                                                     >  7:279163/1‑149 (MQ=255)
                                                                tCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGt                                                                                    >  8:91895/1‑149 (MQ=255)
                                                                  aGCGGGCGGTGCCTCTCGCAGCTTGTGGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGTGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTgg                                                                                  <  4:38960/149‑1 (MQ=255)
                                                                   gcgggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGc                                                                                 <  2:241674/149‑1 (MQ=255)
                                                                          gTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACtt                                                                        >  1:67585/1‑149 (MQ=255)
                                                                          gTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACtt                                                                        <  7:306279/149‑1 (MQ=255)
                                                                                      gCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCt                                                            <  8:279163/149‑1 (MQ=255)
                                                                                          gttgtGGGGTATCAGTGCCTCGGAGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCCGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTcg                                                        >  6:365763/1‑149 (MQ=255)
                                                                                                       aGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTc                                            <  3:57687/148‑1 (MQ=255)
                                                                                                       aGTGCCTAGGCGGGGCTGCTGATGGGGGGAATAATGCATGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCttac                                           >  3:205447/1‑147 (MQ=255)
                                                                                                       aGTCCCTGGGCGGGGCTGCTGATGGGGGGATTAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTc                                            <  5:209558/148‑1 (MQ=255)
                                                                                                         tGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCt                                         >  5:207519/1‑149 (MQ=255)
                                                                                                               ggcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                    <  3:182296/148‑1 (MQ=255)
                                                                                                                gcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGt                                  >  8:116416/1‑149 (MQ=255)
                                                                                                                 cggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTc                                 <  6:123016/149‑1 (MQ=255)
                                                                                                                   gggCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAg                               <  3:331494/149‑1 (MQ=255)
                                                                                                                      ctgctgATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGAc                            <  2:67585/149‑1 (MQ=255)
                                                                                                                       tgctgATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACt                           >  1:199847/1‑149 (MQ=255)
                                                                                                                         cggATGGCGGTAATACTGCCTGCTGCTTCGCATCAGAGACCTTTTGTTTATTGCGCTTCCACAACACGGACTTCTCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTTCCCGCTTCCCTCCGGGGTCAGGACtt                          <  8:313253/146‑1 (MQ=255)
                                                                                                                          tgATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTtgt                        >  3:93969/1‑149 (MQ=255)
                                                                                                                               gggggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                    <  8:120413/148‑1 (MQ=255)
                                                                                                                               gggggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                    <  5:336483/148‑1 (MQ=255)
                                                                                                                               gggggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGtggtgg                   >  2:103344/1‑149 (MQ=255)
                                                                                                                                ggggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGAGAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCTACTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                    <  1:106568/149‑1 (MQ=255)
                                                                                                                                     ataatGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTcc             <  7:146212/149‑1 (MQ=255)
                                                                                                                                     ataatGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGCTATTcc             >  4:212829/1‑149 (MQ=255)
                                                                                                                                      taatGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGA‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTcc             >  4:264184/1‑148 (MQ=255)
                                                                                                                                        gtctctgctgCGAAGCAGCAGGGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat          >  5:14842/5‑149 (MQ=255)
                                                                                                                                        gtctctgctgCGAAGCAGCAGGGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCag          >  5:20040/5‑148 (MQ=255)
                                                                                                                                         tctctgctgCGAAGCAGCAGAGACTTTTTTTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat          >  3:266074/4‑148 (MQ=255)
                                                                                                                                         tctctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGATTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                    >  5:60187/4‑122 (MQ=255)
                                                                                                                                         tctctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGATTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                    <  6:60187/119‑1 (MQ=255)
                                                                                                                                           tctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                    <  5:15912/135‑1 (MQ=255)
                                                                                                                                           tctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                    >  6:15912/2‑136 (MQ=255)
                                                                                                                                           ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGCCTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgat       <  1:62270/149‑1 (MQ=255)
                                                                                                                                           ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgat       <  7:91895/149‑1 (MQ=255)
                                                                                                                                             tgctgcTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATgg     >  2:50120/1‑149 (MQ=255)
                                                                                                                                               ctgctTCGCAGCAGAGACTTTTTGTTTTTTGCGCTACGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTAGCCTGACCGCTTCCCTACGGGGTCAGGACTAGTGGTGGTATTCCCATGATggg    >  3:111159/1‑148 (MQ=255)
                                                                                                                                                  ctTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGa  <  2:224700/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGATGGGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGA  >  NC_002947/3417597‑3417889

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: