Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,255 Δ1 bp 6.4% coding (382/585 nt) nfuA → Fe/S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,2550A. 6.4% 241.7 / 10.4 78coding (382/585 nt)nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base . (2/3);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

GCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGC  >  NC_002947/2719108‑2719402
                                                                                                                                                   |                                                                                                                                                     
gcgcTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACCAAGACAGCCCCATCAACGAACGAATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgccac                                                                                                                                                      >  7:14512/1‑148 (MQ=255)
   cTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgccagcca                                                                                                                                                   >  8:333429/1‑149 (MQ=255)
    tGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCAt                                                                                                                                                  >  1:349348/1‑149 (MQ=255)
       tCGACTACGCCACCGACCGCATGGGTGGCCAGCTGAACATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATAAACGAACGTATCAACTACTACCTACAGACAGAGATCAACCCCGGACTGGCCAGCCAAggc                                                                                                                                               >  8:23503/1‑149 (MQ=255)
         gACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATggctgg                                                                                                                                            >  8:206569/1‑146 (MQ=255)
            taatccACAGACCGCATGGGTGGCCAGCTGACAAACAAGGCGCAAAACGCCAAGGTGCCGATTGTCAAAGAAGACAGCAAGATCAACGAAAGTATAAACTAGAACCTACAGACAGAGATCAACCCCGGCCTTGCCACACATggcggc                                                                                                                                           >  6:373359/5‑147 (MQ=255)
              cGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAgg                                                                                                                                       <  1:13274/149‑1 (MQ=255)
              cGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAgg                                                                                                                                       <  3:323075/149‑1 (MQ=255)
               gCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGCCCGAGGTCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGt                                                                                                                                      >  5:394603/1‑149 (MQ=255)
                    cgaccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGcc                                                                                                                                 <  3:77638/149‑1 (MQ=255)
                    cgaccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGcc                                                                                                                                 >  4:77638/1‑149 (MQ=255)
                    cgaccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggccgggg                                                                                                                                 <  7:206569/149‑7 (MQ=255)
                      accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtg                                                                                                                               <  5:140498/149‑1 (MQ=255)
                      accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtg                                                                                                                               <  3:269923/149‑1 (MQ=255)
                            tGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                        <  7:333429/149‑1 (MQ=255)
                                tGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTgg                                                                                                                    <  2:219002/149‑1 (MQ=255)
                                tGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTgg                                                                                                                    >  1:389828/1‑149 (MQ=255)
                                  gCCATCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGt                                                                                                                   <  6:282166/148‑1 (MQ=255)
                                  gCCAGCTGACCATCTAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGATGACAGCCCGATGAACGAACGTATCAACTACTACCTTCAGACCTAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  <  8:14512/149‑1 (MQ=255)
                                  gCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGACGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  >  2:389341/1‑149 (MQ=255)
                                   ccAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGTCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  <  7:23503/148‑1 (MQ=255)
                                   ccAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  >  7:15856/1‑148 (MQ=255)
                                       cTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGa                                                                                                             <  2:349348/149‑1 (MQ=255)
                                                aaGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTg                                                                                                     <  8:15856/148‑1 (MQ=255)
                                                 aGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGcc                                                                                                   >  2:384282/1‑149 (MQ=255)
                                                 aGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGcc                                                                                                   <  4:382556/149‑1 (MQ=255)
                                                      ccAAACGCCAAGGTGCCGATGGTCAACGAAGACACCCCGATCAACGAACGTATCAACTACTACCTACATACCGAGATCAACCCCGGCCTGGCCATCCATGGGGGCCATGTGTTCCTGGTGGGACGTGGTCGATGACGGCATTGCCGTAc                                                                                               <  1:229365/149‑1 (MQ=255)
                                                       cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTg                                                                                             <  4:18696/149‑1 (MQ=255)
                                                       cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTg                                                                                             >  1:345848/1‑149 (MQ=255)
                                                         aaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                        >  3:233943/1‑120 (MQ=255)
                                                         aaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                        <  4:233943/120‑1 (MQ=255)
                                                           cGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGt                                                                                         >  4:127212/1‑149 (MQ=255)
                                                            gCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGtt                                                                                        >  8:366317/1‑149 (MQ=255)
                                                              cAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCg                                                                                      >  6:175080/1‑149 (MQ=255)
                                                                  gTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGCCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtgg                                                                                  >  6:239237/1‑149 (MQ=255)
                                                                  gTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtgg                                                                                  >  2:251169/1‑149 (MQ=255)
                                                                  gTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtgg                                                                                  <  5:251919/149‑1 (MQ=255)
                                                                   tGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                 <  8:205048/149‑1 (MQ=255)
                                                                     ccGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGgtggtg                                                                                <  7:366317/148‑1 (MQ=255)
                                                                       gATGGTCAACGAAGACAGCCCGATCAACGAACGTATAAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGACCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCt                                                                             >  3:89716/1‑149 (MQ=255)
                                                                            tCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAg                                                                        >  7:276864/1‑149 (MQ=255)
                                                                              aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGtt                                                                                        >  5:149707/1‑131 (MQ=255)
                                                                              aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGtt                                                                                        <  6:149707/131‑1 (MQ=255)
                                                                              aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAggg                                                                      >  5:357899/1‑149 (MQ=255)
                                                                               aCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGt                                                                     >  5:350557/1‑149 (MQ=255)
                                                                                cGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTAc                                                                                               >  4:360515/1‑122 (MQ=255)
                                                                                cGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTAc                                                                                               <  3:360515/122‑1 (MQ=255)
                                                                                 gAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTg                                                                   >  3:4797/1‑149 (MQ=255)
                                                                                   aGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                                 <  3:127212/149‑1 (MQ=255)
                                                                                        gCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAg                                                            >  6:286456/1‑149 (MQ=255)
                                                                                         cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  >  2:217287/1‑94 (MQ=255)
                                                                                         cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTc                                                                                                                  <  1:217287/94‑1 (MQ=255)
                                                                                         cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAgg                                                           >  4:100319/1‑149 (MQ=255)
                                                                                              tCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                      <  2:345848/149‑1 (MQ=255)
                                                                                                   gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGAc                                                                                                            <  5:303340/90‑1 (MQ=255)
                                                                                                   gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGAc                                                                                                            >  6:303340/1‑90 (MQ=255)
                                                                                                    aacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAcc                                                <  8:276864/149‑1 (MQ=255)
                                                                                                         aTCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTga                                            <  1:389341/148‑1 (MQ=255)
                                                                                                         aTCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTg                                             >  2:347096/1‑147 (MQ=255)
                                                                                                          tCAACTACTACCTACAGACCGAGAGCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                           <  6:350557/148‑1 (MQ=255)
                                                                                                             actactACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaa                                       <  2:389828/149‑1 (MQ=255)
                                                                                                              ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                 >  1:145795/1‑106 (MQ=255)
                                                                                                              ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                 <  2:145795/106‑1 (MQ=255)
                                                                                                                   acctacAGACCGAGATCAACCCAGGAATGCCCAGCCATGCCGGCAAGGTGAGCCTGGTGGACGTGGTCGATGACGCAATTGCCGTAATGAAGTAGGGTGGTGCATGCCAGGGTTGCGGCCAGCCCGAGGTCACCCGGAAGGAAGACAt                                  >  7:248027/1‑148 (MQ=255)
                                                                                                                     ctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATc                                 >  8:48904/1‑147 (MQ=255)
                                                                                                                        cAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                             <  1:347096/148‑1 (MQ=255)
                                                                                                                           aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAc                         >  5:234277/1‑149 (MQ=255)
                                                                                                                            ccGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAcc                        <  5:286456/149‑1 (MQ=255)
                                                                                                                              gagaTCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCt                      >  7:99070/1‑149 (MQ=255)
                                                                                                                                    aaCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGa                <  1:384282/149‑1 (MQ=255)
                                                                                                                                     aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                             >  1:283018/1‑135 (MQ=255)
                                                                                                                                     aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                             <  2:283018/135‑1 (MQ=255)
                                                                                                                                     aCCCCGGCCTGGCCAGCCATGGCGACCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTTGTGGCTGCCAGGGTTGCGGCCCGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAg               <  4:89716/149‑1 (MQ=255)
                                                                                                                                      cgccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                      <  1:215716/107‑1 (MQ=255)
                                                                                                                                      cgccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                      >  2:215716/3‑109 (MQ=255)
                                                                                                                                       cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAg                              <  8:217719/132‑1 (MQ=255)
                                                                                                                                       cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAg                              >  7:217719/1‑132 (MQ=255)
                                                                                                                                           gCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTAtt         <  7:48904/149‑1 (MQ=255)
                                                                                                                                             cTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCa                          >  7:98717/1‑130 (MQ=255)
                                                                                                                                             cTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCa                          <  8:98717/130‑1 (MQ=255)
                                                                                                                                                gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCTTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa    <  5:175080/149‑1 (MQ=255)
                                                                                                                                                gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCg      >  6:185510/1‑147 (MQ=255)
                                                                                                                                                gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa    <  3:100319/149‑1 (MQ=255)
                                                                                                                                                  cagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGc  >  6:233471/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGC  >  NC_002947/2719108‑2719402

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: