Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,359 T→C 9.0% intergenic (+22/+639) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,552,3590TC 9.0% 215.2 / 9.1 78intergenic (+22/+639)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base C (3/4);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGG  >  NC_002947/5552215‑5552502
                                                                                                                                                |                                                                                                                                               
gtgtGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGc                                                                                                                                              <  8:48530/148‑1 (MQ=255)
 tgtgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCt                                                                                                                                            <  5:263223/149‑1 (MQ=255)
   tgGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCTGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTgg                                                                                                                                          >  2:216278/1‑149 (MQ=255)
    gggCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTggg                                                                                                                                         <  6:142458/149‑1 (MQ=255)
    gggCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTggg                                                                                                                                         <  1:306084/149‑1 (MQ=255)
     gtgacacactGCCAACGTCTAGTTCGCGGCGTCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGTTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGc                                                                                                                                        <  4:228880/146‑1 (MQ=255)
     ggCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTggg                                                                                                                                         >  2:112730/1‑148 (MQ=255)
     ggCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTggg                                                                                                                                         <  6:262115/148‑1 (MQ=255)
        acacacaGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGcac                                                                                                                                      >  4:261368/1‑148 (MQ=255)
          acacaGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGc                                                                                                                                        <  6:297930/144‑1 (MQ=255)
          acacaGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGc                                                                                                                                        >  5:297930/1‑144 (MQ=255)
               gCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGa                                                                                                                              <  2:54169/149‑1 (MQ=255)
                ccAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTAATTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAg                                                                                                                             >  7:159292/1‑149 (MQ=255)
                 cAACGTCGAGTTCGCGGCGGCGGCGCCGTATGTGGCGTTGTTGATTGTGGGGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGt                                                                                                                            <  3:217986/149‑1 (MQ=255)
                 cAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGt                                                                                                                            >  5:140777/1‑149 (MQ=255)
                  aaCGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTg                                                                                                                           <  6:189032/149‑1 (MQ=255)
                  aaCGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGCATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTg                                                                                                                           >  4:286766/1‑149 (MQ=255)
                      tCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACATGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGcc                                                                                                                         >  6:186480/1‑147 (MQ=255)
                           ttCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACtt                                                                                                                  >  5:261004/1‑149 (MQ=255)
                             cgcgGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTtca                                                                                                                >  5:253533/1‑149 (MQ=255)
                             cgcgGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTtc                                                                                                                 <  1:112730/148‑1 (MQ=255)
                              gcggcggcggcgCCGTATGCGGCGGTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTtcat                                                                                                               <  5:24026/149‑1 (MQ=255)
                                      ggcgCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGcaa                                                                                                                                      >  4:125783/1‑117 (MQ=255)
                                      ggcgCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGcaa                                                                                                                                      <  3:125783/118‑2 (MQ=255)
                                        cgcCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCa                                                                                                     <  3:85509/149‑1 (MQ=255)
                                           cGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACtt                                                                                                  <  1:40343/149‑1 (MQ=255)
                                                 cggcgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTaaaa                                                                                             <  2:161902/148‑1 (MQ=255)
                                                 cggcgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGCATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTaaaa                                                                                             <  3:286766/148‑1 (MQ=255)
                                                    cgTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAActct                                                                                         <  6:261004/149‑1 (MQ=255)
                                                        gttgATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCAc                                                                                      >  4:71061/1‑148 (MQ=255)
                                                           gATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAg                                                                                   <  6:140777/148‑1 (MQ=255)
                                                           gATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGc                                                                                  >  8:469/1‑149 (MQ=255)
                                                                tggtgTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACTAGGCTTGTAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGtt                                                                             <  5:83981/149‑1 (MQ=255)
                                                                tggtgTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGtt                                                                             <  3:71061/149‑1 (MQ=255)
                                                                tggtgTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGtt                                                                             <  6:54970/149‑1 (MQ=255)
                                                                   tgtCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGa                                                                           >  3:285441/1‑148 (MQ=255)
                                                                   tgtCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGa                                                                           >  5:1326/1‑148 (MQ=255)
                                                                    gtCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCa                                                                                       <  5:125773/135‑1 (MQ=255)
                                                                    gtCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCa                                                                                       >  6:125773/1‑135 (MQ=255)
                                                                     tCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCa                                                                                                     <  7:155037/120‑1 (MQ=255)
                                                                     tCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCa                                                                                                     >  8:155037/1‑120 (MQ=255)
                                                                       ggggTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGc                                                                       >  3:172376/1‑148 (MQ=255)
                                                                              ccGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAAc                                                                <  4:172376/148‑1 (MQ=255)
                                                                               cGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAACCTCTCACTAGCCAGTTGAAAGCAAAGAACaa                                                              <  8:159292/149‑1 (MQ=255)
                                                                               cGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAAc                                                                >  6:169006/1‑147 (MQ=255)
                                                                                 gtgtATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAg                                                             >  5:88203/1‑148 (MQ=255)
                                                                                    tATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTc                                                          <  6:1326/148‑1 (MQ=255)
                                                                                      tttGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTa                                                       <  4:285441/149‑1 (MQ=255)
                                                                                           tGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGAtac                                                  >  2:293208/1‑149 (MQ=255)
                                                                                                         tATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCcccaccc                                    <  3:261368/149‑1 (MQ=255)
                                                                                                         tATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCcccacc                                     >  4:139161/1‑148 (MQ=255)
                                                                                                          aTTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCcccaccc                                    >  4:248685/1‑148 (MQ=255)
                                                                                                           tttGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCa                                                                                <  2:308059/103‑1 (MQ=255)
                                                                                                           tttGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCa                                                                                >  1:308059/1‑103 (MQ=255)
                                                                                                              gAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAagag                               >  2:214743/1‑149 (MQ=255)
                                                                                                               aaCAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCAc                                                                                      >  8:206882/1‑93 (MQ=255)
                                                                                                               aaCAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAgg                              >  6:19916/1‑149 (MQ=255)
                                                                                                               aaCAAGGCTTGAAAGCTTGCACAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTCAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAagag                               <  1:216278/148‑1 (MQ=255)
                                                                                                                aCAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAg                                                                                                                             <  8:273498/53‑1 (MQ=255)
                                                                                                                aCAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAg                                                                                                                             >  7:273498/1‑53 (MQ=255)
                                                                                                                 cAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAggg                             >  2:140485/1‑148 (MQ=255)
                                                                                                                    ggCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTaa                          >  7:203844/1‑148 (MQ=255)
                                                                                                                       ttGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCATAGAACAAGGTCGTAGATACTACAAGc                                           <  1:191236/128‑1 (MQ=255)
                                                                                                                       ttGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGc                                           >  2:191236/1‑128 (MQ=255)
                                                                                                                         gAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTaa                          <  1:210331/143‑1 (MQ=255)
                                                                                                                         gAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTaa                          >  2:210331/1‑143 (MQ=255)
                                                                                                                         gAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAAc                    >  1:228656/1‑149 (MQ=255)
                                                                                                                              cTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACaaaa                     <  6:78641/143‑1 (MQ=255)
                                                                                                                              cTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACaaaa                     >  5:78641/1‑143 (MQ=255)
                                                                                                                              cTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCttt               <  6:88203/149‑1 (MQ=255)
                                                                                                                                 gCCCAGGCCACATGGCGGCCTGGGCCCAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGcc                                          >  6:40584/1‑119 (MQ=255)
                                                                                                                                 gCCCAGGCCACATGGCGGCCTGGGCCCAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGcc                                          <  5:40584/119‑1 (MQ=255)
                                                                                                                                   ccAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATa          <  3:139161/149‑1 (MQ=255)
                                                                                                                                    cAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAg         >  6:249532/1‑149 (MQ=255)
                                                                                                                                       gccgccATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGAcc      >  2:20117/1‑149 (MQ=255)
                                                                                                                                         cgccATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCtt    >  4:290257/1‑149 (MQ=255)
                                                                                                                                          gccATGTGGGCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTg   >  8:256955/1‑149 (MQ=255)
                                                                                                                                           ccATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTgg  <  3:290257/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGG  >  NC_002947/5552215‑5552502

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: