Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,499,711 G→T 5.6% D538Y (GAC→TAC)  PP_2193 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,499,7110GT 5.6% 226.9 / 3.5 71D538Y (GAC→TAC) PP_2193ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

CTTCAACGCCGAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACGGCCACCAG  >  NC_002947/2499563‑2499857
                                                                                                                                                    |                                                                                                                                                  
cTTCAACGCCGAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCg                                                                                                                                                    <  2:19356/149‑1 (MQ=255)
 ttCAACGCCGAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCAGGAGCTGGAAAAAGACAGCAAGGCCGCCCGCTACCCCAGCTACAGGGTGGATGGTGTGTCGGTGTATGCCAAGCTGGACTGGCATGCCACCGa                                                                                                                                                   >  6:76453/1‑149 (MQ=255)
    aaCGCCGAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACaa                                                                                                                                                >  8:237488/1‑149 (MQ=255)
      cGCCGAGCAGACCACCTTCAACGCACAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTAACGGCTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAAcc                                                                                                                                              >  5:438577/1‑149 (MQ=255)
          gAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGAc                                                                                                                                          <  6:33399/149‑1 (MQ=255)
                accaccTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTc                                                                                                                                    <  7:237488/149‑1 (MQ=255)
                  caccTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCg                                                                                                                                  >  8:53359/1‑149 (MQ=255)
                    ccTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCggt                                                                                                                                >  8:378132/1‑149 (MQ=255)
                     cTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCggtt                                                                                                                               <  8:382232/149‑2 (MQ=255)
                       tCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCggtggt                                                                                                                             >  5:236088/1‑149 (MQ=255)
                            gCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCcg                                                                                                                        >  5:269241/1‑149 (MQ=255)
                             cccAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCcgg                                                                                                                       <  7:378132/149‑1 (MQ=255)
                             cccAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCcgg                                                                                                                       <  3:382525/149‑1 (MQ=255)
                               cAGCCCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCcggcg                                                                                                                     <  2:133113/149‑1 (MQ=255)
                                  aCCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGgca                                                                                                                  <  1:311524/149‑1 (MQ=255)
                                   ccTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGgcag                                                                                                                 <  6:411365/149‑1 (MQ=255)
                                   ccTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGgcag                                                                                                                 <  6:92376/149‑1 (MQ=255)
                                       gTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGa                                                                                                             <  1:378756/149‑1 (MQ=255)
                                           gAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGAGGCTGCAGATGGa                                                                                                         >  7:287801/1‑149 (MQ=255)
                                              tCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGgcagc                                                                                                                <  4:341616/139‑1 (MQ=255)
                                              tCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGgcagc                                                                                                                >  3:341616/1‑139 (MQ=255)
                                              tCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGt                                                                                                      <  6:363293/149‑1 (MQ=255)
                                                 ggTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGa                                                                                                   >  4:356568/1‑149 (MQ=255)
                                                  gTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGAt                                                                                                  >  3:279001/1‑149 (MQ=255)
                                                  gTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGAt                                                                                                  <  6:236088/149‑1 (MQ=255)
                                                                 aaaaaGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGc                                                                                            <  1:148487/140‑1 (MQ=255)
                                                                 aaaaaGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGc                                                                                            >  2:148487/1‑140 (MQ=255)
                                                                  aaaaGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAgg                                                                                  >  5:170593/1‑149 (MQ=255)
                                                                         aGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGc                                                                            >  4:49104/1‑148 (MQ=255)
                                                                              ggCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTg                                                                      <  7:105265/149‑1 (MQ=255)
                                                                              ggCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTg                                                                      >  4:255392/1‑149 (MQ=255)
                                                                                  gACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGgt                                                                                >  2:146229/1‑135 (MQ=255)
                                                                                  gACCGCTACCCCAGCTACCGGGTGGATGGTGAGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGgt                                                                                <  1:146229/135‑1 (MQ=255)
                                                                                     cgatACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGATCGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGAGGGACGTTGATTTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTaa                                                               <  1:277351/146‑1 (MQ=255)
                                                                                     cGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTaa                                                               <  4:279001/149‑1 (MQ=255)
                                                                                       cTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTgg                                                                     <  7:30293/141‑1 (MQ=255)
                                                                                       cTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTgg                                                                     >  8:30293/1‑141 (MQ=255)
                                                                                          ccccAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGa                                                            <  6:438577/147‑1 (MQ=255)
                                                                                          ccccAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGa                                                            <  3:49104/147‑1 (MQ=255)
                                                                                          ccccAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGa                                                            <  7:53359/147‑1 (MQ=255)
                                                                                          ccccAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGa                                                            <  3:356568/147‑1 (MQ=255)
                                                                                            ccAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTAt                                                        <  1:160719/149‑1 (MQ=255)
                                                                                             cAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTAt                                                        >  3:330237/1‑148 (MQ=255)
                                                                                               gCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAAACTGACATTGTCCGGTGGTGCCCGGCGGCATCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAg                                                     <  2:127187/149‑1 (MQ=255)
                                                                                                   ccGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGt                                                   <  3:255392/147‑1 (MQ=255)
                                                                                                   ccGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTc                                                  >  2:149066/1‑148 (MQ=255)
                                                                                                          gATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGt                                           <  6:170593/148‑1 (MQ=255)
                                                                                                                 tgtCGGTGTATGCCCAGCTGGACTGGCATGCCACCTACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAAc                                               <  2:95364/137‑1 (MQ=255)
                                                                                                                 tgtCGGTGTATGCCCAGCTGGACTGGCATGCCACCTACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAAc                                               <  2:101433/137‑1 (MQ=255)
                                                                                                                 tgtCGGTGTATGCCCAGCTGGACTGGCATGCCACCTACAACCTGACAATGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAAc                                               >  1:101433/1‑137 (MQ=255)
                                                                                                                 tgtCGGTGTATGCCCAGCTGGACTGGCATGCCACCTACAACCTGACAATGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAAc                                               >  1:95364/1‑137 (MQ=255)
                                                                                                                 tgtCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTc                                   >  4:120033/1‑149 (MQ=255)
                                                                                                                     ggTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGt                                <  3:235123/148‑1 (MQ=255)
                                                                                                                        gtATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCa                            >  2:151492/1‑149 (MQ=255)
                                                                                                                         tATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCAt                           <  6:152320/149‑1 (MQ=255)
                                                                                                                          aTGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATc                          >  1:122933/1‑149 (MQ=255)
                                                                                                                            gCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTa                        >  7:155980/1‑149 (MQ=255)
                                                                                                                            gCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTa                        >  3:242294/1‑149 (MQ=255)
                                                                                                                              ccAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGc                              <  7:286754/141‑1 (MQ=255)
                                                                                                                              ccAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGc                              >  8:286754/1‑141 (MQ=255)
                                                                                                                                aGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGt                                >  2:177269/1‑137 (MQ=255)
                                                                                                                                aGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGt                                <  1:177269/137‑1 (MQ=255)
                                                                                                                                   tGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACtt                  >  2:141539/1‑148 (MQ=255)
                                                                                                                                     gACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTc                 <  1:61414/147‑1 (MQ=255)
                                                                                                                                      aCTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTc                 >  1:357402/1‑146 (MQ=255)
                                                                                                                                          gCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAAc           <  4:82075/148‑1 (MQ=255)
                                                                                                                                          gCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACg          >  7:207958/1‑149 (MQ=255)
                                                                                                                                          gCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACg          >  2:283939/1‑149 (MQ=255)
                                                                                                                                           cATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACgg         <  6:269241/149‑1 (MQ=255)
                                                                                                                                               ccACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCa                                                                   >  4:139742/1‑87 (MQ=255)
                                                                                                                                               ccACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCa                                                                   <  3:139742/87‑1 (MQ=255)
                                                                                                                                                  ccGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACGGCCACCAg  <  4:242294/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTCAACGCCGAGCAGACCACCTTCAACGCCCAGACCTCGTCCGAGTCGGGTGGCCTGGAGCTGGAAAAAGACAGCAAGGCCGACCGCTACCCCAGCTACCGGGTGGATGGTGTGTCGGTGTATGCCCAGCTGGACTGGCATGCCACCGACAACCTGACATTGTCCGGTGGTGCCCGGCGGCAGCAGATGGACGTTGATGTCAGCGACTTCAAGGGTGTGCCAGGTGGCAGTAACGACTATCAGGTCAACCTGTTCAACGTCGGTGCCATCTACGACTTCAAGAACGGCCACCAG  >  NC_002947/2499563‑2499857

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