Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,417,746 C→A 10.8% E99* (GAA→TAA) ‡ PP_3030 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,417,7460CA10.8% 177.8 / 7.8 64E99* (GAA→TAA) ‡PP_3030hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (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

GTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGCGA  >  NC_002947/3417603‑3417894
                                                                                                                                               |                                                                                                                                                      
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                                                                             <  7:306279/149‑1 (MQ=255)
                                                                    gTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACtt                                                                             >  1:67585/1‑149 (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                                         <  6:60187/119‑1 (MQ=255)
                                                                                                                                   tctctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGATTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                         >  5:60187/4‑122 (MQ=255)
                                                                                                                                     tctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                         >  6:15912/2‑136 (MQ=255)
                                                                                                                                     tctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                         <  5:15912/135‑1 (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)
                                                                                                                                               cGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGCGa  >  1:167289/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
GTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGC‑‑CTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGCGA  >  NC_002947/3417603‑3417894

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