Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,391,522 G→A 8.3% W244* (TGG→TGA) ‡ PP_2994 → FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,391,5220GA 8.3% 176.6 / 4.1 60W244* (TGG→TGA) ‡PP_2994FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (33/22);  new base A (2/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391374‑3391670
                                                                                                                                                    |                                                                                                                                                    
tGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTgg                                                                                                                                                      >  1:92833/1‑149 (MQ=255)
         gggCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAc                                                                                                                                             <  4:185328/149‑1 (MQ=255)
             gcgcCTGATGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCAGCCGATTTCCAGCGCTGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTgg                                                                                                                                         <  5:240481/149‑1 (MQ=255)
                   gTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGAGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCg                                                                                                                                      >  1:236175/1‑146 (MQ=255)
                       cTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAt                                                                                                                               >  5:367006/1‑149 (MQ=255)
                       cTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCAAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAt                                                                                                                               >  5:422209/1‑149 (MQ=255)
                        tCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAcc                                                                                                                                            <  8:380434/135‑1 (MQ=255)
                        tCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAcc                                                                                                                                            >  7:380434/1‑135 (MQ=255)
                         cGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                             >  6:397231/1‑149 (MQ=255)
                           aaaTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAtggtg                                                                                                                           >  6:244527/1‑149 (MQ=255)
                                          tCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGc                                                                                                            <  6:317220/149‑1 (MQ=255)
                                              ggCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGACCtt                                                                                                                                          >  4:332309/1‑114 (MQ=255)
                                              ggCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGACCtt                                                                                                                                          <  3:332309/115‑2 (MQ=255)
                                                      tttCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         <  5:112354/124‑1 (MQ=255)
                                                      tttCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         >  6:112354/1‑124 (MQ=255)
                                                        tCTGTGGATTTTGCGGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGAc                                                                                                                                 >  6:137825/1‑114 (MQ=255)
                                                        tCTGTGGATTTTGCGGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGAc                                                                                                                                 <  5:137825/114‑1 (MQ=255)
                                                          tgtgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCg                                                                                            >  3:224726/1‑149 (MQ=255)
                                                           gtgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                           >  8:304602/1‑149 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAg                                                                                                                                   <  4:116010/108‑1 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAg                                                                                                                                   >  3:116010/1‑108 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAg                                                                                                  >  6:13750/1‑141 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAg                                                                                                  <  5:13750/141‑1 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                           >  8:89842/1‑148 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                           <  5:244527/148‑1 (MQ=255)
                                                            tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                           <  2:236175/148‑1 (MQ=255)
                                                                  ttGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCa                                                                                     <  4:224726/148‑1 (MQ=255)
                                                                    gCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAgg                                                                                   <  7:304602/148‑1 (MQ=255)
                                                                      aGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTg                                                                                <  6:367006/149‑1 (MQ=255)
                                                                          gCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATc                                                                                                                   <  5:222902/110‑1 (MQ=255)
                                                                          gCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATc                                                                                                                   >  6:222902/1‑110 (MQ=255)
                                                                               cAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCg                                                                                                                  <  3:100243/106‑1 (MQ=255)
                                                                               cAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCg                                                                                                                  >  4:100243/1‑106 (MQ=255)
                                                                                aaaCTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGAcgcgcgcg                                                                       <  7:89842/148‑1 (MQ=255)
                                                                                aaaCTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGa                                                                      >  4:338670/1‑149 (MQ=255)
                                                                                 aactcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGa                                                                                                               <  8:148573/107‑1 (MQ=255)
                                                                                  actcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAgg                                                                                   >  4:407998/1‑134 (MQ=255)
                                                                                  actcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAgg                                                                                   <  3:407998/134‑1 (MQ=255)
                                                                                    tcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGt                                                                  >  7:148757/1‑149 (MQ=255)
                                                                                    tcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGt                                                                  >  7:377319/1‑149 (MQ=255)
                                                                                        ctcCGCCGATTTCCATCGCGGAGAGTTCAACACAGAAGACGCCAAAAAGGTCGTTGAGTAGCTAAACGACCTGAGCGTAGACATGGTGGAGTTATAGGTAGGCAGATACGAAGTAAACGAAATACAGGATGACGCGCGCTATGGTcac                                                               >  8:257558/1‑146 (MQ=255)
                                                                                           cgccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCAccc                                                           >  2:396889/1‑149 (MQ=255)
                                                                                            gccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCACCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCt                                                          >  2:521/1‑149 (MQ=255)
                                                                                             ccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTg                                                         >  3:107514/1‑149 (MQ=255)
                                                                                                     cAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGa                                                 <  3:338670/149‑1 (MQ=255)
                                                                                                        cgcgGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCCCGAAGc                                              <  8:377319/149‑1 (MQ=255)
                                                                                                                 ttCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgt                                     >  8:233835/1‑149 (MQ=255)
                                                                                                                   caccacCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    <  6:287043/148‑1 (MQ=255)
                                                                                                                   caccacCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    <  5:397231/148‑1 (MQ=255)
                                                                                                                    accaccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGTTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    >  1:280206/1‑147 (MQ=255)
                                                                                                                    accaccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    >  4:257716/1‑147 (MQ=255)
                                                                                                                     ccaccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    <  4:10924/146‑1 (MQ=255)
                                                                                                                     ccaccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    >  3:10924/1‑146 (MQ=255)
                                                                                                                      caccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGt                                                                  <  7:66917/115‑1 (MQ=255)
                                                                                                                      caccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGt                                                                  >  8:66917/1‑115 (MQ=255)
                                                                                                                      caccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGt                                >  7:80770/1‑149 (MQ=255)
                                                                                                                            aGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGccc                          >  8:364236/1‑149 (MQ=255)
                                                                                                                                gcgAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGa                      >  5:184128/1‑149 (MQ=255)
                                                                                                                                   aaaaaGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtct                   >  5:107613/1‑149 (MQ=255)
                                                                                                                                   aaaaaGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtct                   >  3:299051/1‑149 (MQ=255)
                                                                                                                                      aaGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCACCGACATGCAGGGTGACTCGCGCGATGGTCGCACCCTGGCTAGAGAAGCTAATTTCGTCGAGTTCGCCCGTGAtctgc                 >  6:77461/1‑146 (MQ=255)
                                                                                                                                                    gCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACTCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCaaaa  >  8:125159/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391374‑3391670

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