Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,391,522 G→A 10.7% 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,5220GA10.7% 173.2 / 9.5 65W244* (TGG→TGA) ‡PP_2994FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base A (5/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391376‑3391670
                                                                                                                                                    |                                                                                                                                                    
gAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGc                                                                                                                                                     <  4:369052/148‑1 (MQ=255)
      cGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGa                                                                                                                                              <  1:69418/149‑1 (MQ=255)
      cGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGa                                                                                                                                              <  5:274911/149‑1 (MQ=255)
           gcgcCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTgg                                                                                                                                         <  5:364634/149‑1 (MQ=255)
                 gTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCg                                                                                                                                      >  8:240514/1‑146 (MQ=255)
                    gCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACa                                                                                                                                >  5:47649/1‑149 (MQ=255)
                      tCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                              >  1:118447/1‑149 (MQ=255)
                        gAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggt                                                                                                                            >  2:107155/1‑149 (MQ=255)
                         aaaTCGTCAAGGCGGGCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggtg                                                                                                                           <  6:267582/149‑1 (MQ=255)
                          aaTCGTCAAGGCGGTCCTGGCCGTGGTTGCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCTCGGAGGGTTCACCACCGAAGACGCGAAAAAGGGCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggtgg                                                                                                                          <  8:51707/149‑1 (MQ=255)
                           aTCGTCAAGGCGGTCCGGTCCGTGGTTTCTGTGGATTTTGCAGTGGCGCTCAAACTCAACTCCGCCGATTTCCAGCGCGGTGGGTTCACCACCGAAGACGCGAAAAAGGGCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         <  6:293836/149‑1 (MQ=255)
                           aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         <  1:107155/149‑1 (MQ=255)
                           aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         <  3:296578/149‑1 (MQ=255)
                            tCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAg                                                                                                                        >  6:311998/1‑149 (MQ=255)
                               tCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCTCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGAc                                                                                                                                             >  1:317935/1‑125 (MQ=255)
                               tCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGAc                                                                                                                                             <  2:317935/125‑1 (MQ=255)
                                     cGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                             >  1:54830/1‑135 (MQ=255)
                                     cGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                             <  2:54830/135‑1 (MQ=255)
                                      ggTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAg                                                                                                              <  7:33812/149‑1 (MQ=255)
                                                 tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                              >  1:352957/1‑122 (MQ=255)
                                                 tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                              <  2:352957/122‑1 (MQ=255)
                                                 tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                                   <  1:10225/149‑1 (MQ=255)
                                                          tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                                   >  8:105268/1‑140 (MQ=255)
                                                          tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                                   <  7:105268/140‑1 (MQ=255)
                                                          tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                           >  2:246343/1‑148 (MQ=255)
                                                          tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGCAGGCAGCTACGAAGTCCCCGcc                                                                                          <  7:240514/149‑1 (MQ=255)
                                                             aTTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATg                                                                                       >  5:187018/1‑149 (MQ=255)
                                                              ttttGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGc                                                                                      <  5:311998/149‑1 (MQ=255)
                                                                ttGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGACT‑‑GGCGCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCa                                                                                     <  4:2056/148‑1 (MQ=255)
                                                                 tgaaGTGCCGGTAAAACTCAACTCCGCCTATTTCCAGCCCGGTGGGTTAACCACCGATTACCCTAAAAAGGTCGTTTTTTCCTGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGc                                                                                      <  1:113071/145‑1 (MQ=255)
                                                                 tGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAg                                                                                    <  2:130605/148‑1 (MQ=255)
                                                                   cAGTGGCTGTCAAACTCAACTCCGGCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAATAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGt                                                                                 <  7:342972/149‑1 (MQ=255)
                                                                   cAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGt                                                                                 <  2:118447/149‑1 (MQ=255)
                                                                              aaaCTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGa                                                                      <  1:246343/149‑1 (MQ=255)
                                                                                   caactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATgg                                                                   >  2:57140/1‑147 (MQ=255)
                                                                                   caactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATgg                                                                   <  1:57140/147‑1 (MQ=255)
                                                                                        ccgccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGTTGACGCGCGCGATGGTCGCAcc                                                            >  8:58589/1‑149 (MQ=255)
                                                                                           ccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTg                                                         <  6:47649/149‑1 (MQ=255)
                                                                                              aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         >  5:221765/1‑82 (MQ=255)
                                                                                              aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                         <  6:221765/82‑1 (MQ=255)
                                                                                              aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCt                                                      <  3:213719/149‑1 (MQ=255)
                                                                                                ttCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTcg                                                    >  6:297484/1‑149 (MQ=255)
                                                                                                         ggAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTa                                           >  4:92254/1‑149 (MQ=255)
                                                                                                             ggTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtc                                       >  4:214791/1‑149 (MQ=255)
                                                                                                             ggTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtc                                       >  6:214620/1‑149 (MQ=255)
                                                                                                               ttCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCAACGCCATGCAGTGTGACGCGCGCGATCGAAGCACACTGGCTCGCGAAGCTTATTacgt                                     >  6:248494/1‑149 (MQ=255)
                                                                                                                 caccacCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                                    <  3:214791/148‑1 (MQ=255)
                                                                                                                  accaccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                                >  2:307204/1‑135 (MQ=255)
                                                                                                                  accaccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                                <  1:307204/135‑1 (MQ=255)
                                                                                                                    caccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGt                                <  5:214620/149‑1 (MQ=255)
                                                                                                                      ccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTc                              >  6:35773/1‑149 (MQ=255)
                                                                                                                       cGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCg                             >  1:178673/1‑149 (MQ=255)
                                                                                                                       cGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCg                             <  8:173901/149‑1 (MQ=255)
                                                                                                                           gACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCg                         >  5:379990/1‑149 (MQ=255)
                                                                                                                               cgAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAt                     >  5:140767/1‑149 (MQ=255)
                                                                                                                                gAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtc                    <  3:92254/149‑1 (MQ=255)
                                                                                                                                  aaaaGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtctc                  >  5:283680/1‑149 (MQ=255)
                                                                                                                                   ccaGGTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtctc                  >  2:77394/3‑148 (MQ=255)
                                                                                                                                     aGGTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTACCAGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCa                >  5:272698/1‑148 (MQ=255)
                                                                                                                                       gTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc             >  2:116430/1‑149 (MQ=255)
                                                                                                                                       gTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc             >  7:187372/1‑149 (MQ=255)
                                                                                                                                        tCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAcc            >  8:46931/1‑149 (MQ=255)
                                                                                                                                           ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc         <  2:178673/149‑1 (MQ=255)
                                                                                                                                           ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc         >  2:362380/1‑149 (MQ=255)
                                                                                                                                           ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc         <  1:362380/149‑1 (MQ=255)
                                                                                                                                                    gCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCaaaa  >  4:87278/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391376‑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: ‑