Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,355 C→A 9.1% intergenic (+18/+643) 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,3550CA 9.1% 231.1 / 10.9 77intergenic (+18/+643)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base A (3/4);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

CACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGG  >  NC_002947/5552208‑5552502
                                                                                                                                                   |                                                                                                                                                   
cACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCa                                                                                                                                                    >  5:142458/1‑149 (MQ=255)
cACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCa                                                                                                                                                    >  1:54169/1‑149 (MQ=255)
 aCCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCAt                                                                                                                                                   >  5:54970/1‑149 (MQ=255)
       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                                                                                                                                        >  5:297930/1‑144 (MQ=255)
                 acacaGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGc                                                                                                                                        <  6:297930/144‑1 (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                                                                                                                                      <  3:125783/118‑2 (MQ=255)
                                             ggcgCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCACATGGCGGCCTGGGcaa                                                                                                                                      >  4:125783/1‑117 (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)
                                                                                                                                                   |                                                                                                                                                   
CACCGAGGTGTGGGCACACACAGCCAACGTCGAGTTCGCGGCGGCGGCGCCGTATGCGGCGTTGTTGATTGTGGTGTCGGGGTTGCCGGTGTATTTGCTGACTACGCGGATGTATTTGAACAAGGCTTGAAAGCTTGCCCAGGCCGCCATGTGGCCTGGGCACAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGG  >  NC_002947/5552208‑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: