Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 425,349 A→T 6.8% intergenic (‑40/+17) PP_0349 ← / ← PP_0350 membrane protein/ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947425,3490AT 6.8% 214.9 / 5.3 74intergenic (‑40/+17)PP_0349/PP_0350membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

AGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCA  >  NC_002947/425204‑425496
                                                                                                                                                  |                                                                                                                                                   
aGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggc                                                                                                                                                  <  4:295749/149‑1 (MQ=255)
aGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggc                                                                                                                                                  <  2:388244/149‑1 (MQ=255)
 gATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcg                                                                                                                                                 <  7:159712/149‑1 (MQ=255)
   tAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAgg                                                                                                                                                   <  1:208685/145‑1 (MQ=255)
   tAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAgg                                                                                                                                                   >  2:208685/1‑145 (MQ=255)
      ggggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGcccc                                                                                                                                            <  3:170150/149‑1 (MQ=255)
       gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccg                                                                                                                                           <  1:72106/149‑2 (MQ=255)
         gCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGGCAGGCGGCCCCAGc                                                                                                                                         >  2:391338/1‑149 (MQ=255)
            ggCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctca                                                                                                                                      <  5:347516/149‑1 (MQ=255)
              ccGGTGACTGCCAGCAACAGCATTTAGAGGGTACAGACCAGGCTGCTCCAGGTGTGGTGCAAGGACCAGCGGTGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCTAGGCGGccc                                                                                                                                             <  8:154554/140‑1 (MQ=255)
              ccGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGAGTGTTGTAGGTTGTAACCCGGGGCCGCCAGGGGGCCGCacatccgc                                                                                                                                    >  7:154554/1‑142 (MQ=255)
                ggTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGTAAGGCGGCCCCAGctcact                                                                                                                                    >  8:216506/1‑147 (MQ=255)
                      tGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCGCTCACTCAGAAcgt                                                                                                                            >  4:52872/1‑147 (MQ=255)
                        ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGGGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                            <  7:183325/147‑1 (MQ=255)
                          aGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGTGGCCCCAGCTCACTCAGAACTTGt                                                                                                                          >  4:252432/1‑147 (MQ=255)
                                   agGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTg                                                                                                               <  2:396536/149‑1 (MQ=255)
                                   agGAACAGGGTACAGACCAGGCTGCTCCAGGTGTAGATCACGGACCACCGGCTAAAATTATAAAAATTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGGAAGTGGCGCTg                                                                                                               >  4:392299/1‑149 (MQ=255)
                                           ggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAAc                                                                                                       >  1:83989/1‑149 (MQ=255)
                                             tACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGGCCCAGCTCACTCAGAACTTGTAAGTGGGGCTGGCGACAACGc                                                                                                     >  3:401627/1‑149 (MQ=255)
                                               cAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCt                                                                                                    <  3:244248/148‑1 (MQ=255)
                                                  aCCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCCGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                                >  8:122895/1‑149 (MQ=255)
                                                   ccAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                                >  6:337407/1‑148 (MQ=255)
                                                            tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCAGGCGGCCCCCGCGGACGCAGAACGTGTTAGTGGCGCGGGCGACAACGCTGCGCTCAGCGCCgg                                                                                      >  4:22455/1‑148 (MQ=255)
                                                            tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATAACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCg                                                                                       >  5:333986/1‑148 (MQ=255)
                                                              cAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGTGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCgt                                                                                      <  3:252432/147‑1 (MQ=255)
                                                                ggTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCgta                                                                                     <  7:261291/146‑1 (MQ=255)
                                                                      gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGTAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagta                                                                            <  7:216506/149‑1 (MQ=255)
                                                                        tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                                >  4:154523/1‑127 (MQ=255)
                                                                        tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                                <  3:154523/127‑1 (MQ=255)
                                                                        tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                           >  4:41306/1‑148 (MQ=255)
                                                                        tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                           >  1:335590/1‑148 (MQ=255)
                                                                        tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACACAGAACATGTAAGTTGCGCTGGCGAAAACGCTGCGCACATCGCCGTAGTAGCagtat                                                                           >  6:271952/1‑147 (MQ=255)
                                                                         cAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                           >  1:331227/1‑147 (MQ=255)
                                                                          agggACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGtagtag                                                                        <  3:83222/147‑1 (MQ=255)
                                                                          aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAgtagta                                                                         <  2:83989/148‑1 (MQ=255)
                                                                          aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGtagtag                                                                        <  5:337407/149‑1 (MQ=255)
                                                                              aCCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                      <  3:35304/147‑1 (MQ=255)
                                                                              aCCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                      >  4:35304/1‑147 (MQ=255)
                                                                                 aGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTcac                                                                 <  8:285871/149‑1 (MQ=255)
                                                                                  gcggcgGATGGTGGGGCTCTTCATCACGGTACCAGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTca                                                                  >  1:262864/1‑148 (MQ=255)
                                                                                    ggcggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAg                                                              <  1:213678/149‑1 (MQ=255)
                                                                                         aTGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagag                                                         <  3:41306/149‑1 (MQ=255)
                                                                                          tggtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTTAGTGGCGCTGGCGACAACGCTGCGGTCATCGCCGTAGTAGCAGTAGGAGGAGGAACAATTagagc                                                        >  4:235131/1‑148 (MQ=255)
                                                                                            gtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                       <  5:271952/148‑1 (MQ=255)
                                                                                             tggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAgg                                                                       <  1:62127/131‑1 (MQ=255)
                                                                                             tggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAgg                                                                       >  2:62127/1‑131 (MQ=255)
                                                                                               gggCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTaa                                                   >  5:117554/1‑149 (MQ=255)
                                                                                                 gCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagta                                                                            <  7:331248/122‑1 (MQ=255)
                                                                                                  ctctTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATcc                                                <  6:117554/149‑1 (MQ=255)
                                                                                                       catcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGt                                            <  6:333986/148‑1 (MQ=255)
                                                                                                        atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                            >  8:210665/1‑67 (MQ=255)
                                                                                                        atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                            <  7:210665/67‑1 (MQ=255)
                                                                                                           aCGGTACCGATTGTTGTCGGTTGGAACCGGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcag                                        >  1:230332/1‑148 (MQ=255)
                                                                                                             ggTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCACAACTTGTAAGTGGCCCAGGCGACAACGCTGAGCTCATCGCCGGAGTAGCAGAAGTAGGAGTCACAACTAGAGAGGTAATACATGTTcagcag                                     >  2:110849/1‑149 (MQ=255)
                                                                                                               tACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagaga                                                        >  5:159406/1‑128 (MQ=255)
                                                                                                               tACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagaga                                                        <  6:159406/128‑1 (MQ=255)
                                                                                                                  cGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                           >  4:74216/1‑106 (MQ=255)
                                                                                                                  cGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                           <  3:74216/106‑1 (MQ=255)
                                                                                                                    gtttgttgTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggt                               <  1:110849/147‑1 (MQ=255)
                                                                                                                    gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggt                               <  6:117061/149‑1 (MQ=255)
                                                                                                                           gtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTaa                                                                                                                        <  5:222270/53‑1 (MQ=255)
                                                                                                                                   gAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGCTAATCCTTGTTCAGCAGGTTGGTGGCGGGGACCGcga                 <  3:392299/148‑1 (MQ=255)
                                                                                                                                     aCCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgacc               >  8:141121/1‑148 (MQ=255)
                                                                                                                                      gctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccga             >  5:225848/4‑149 (MQ=255)
                                                                                                                                      cccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccga             >  1:376251/1‑149 (MQ=255)
                                                                                                                                       ctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac            >  7:241022/3‑149 (MQ=255)
                                                                                                                                       ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac            >  4:28424/1‑149 (MQ=255)
                                                                                                                                          gggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGAcgcg         <  8:241022/149‑1 (MQ=255)
                                                                                                                                              cGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGGAGCAGTAGGAGGAGTCACAAGGAGCGAGGGACTCCGTGTTCCGCCGGGTGGTGGCGGTGCCCGCGACCGACGCGCCCt     >  5:125202/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagag                                                         <  6:311531/96‑1 (MQ=255)
                                                                                                                                               gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc    >  4:213402/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc    >  7:206930/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc    >  2:284733/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCGGCCCCAGCTCAATCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc    >  2:240199/1‑149 (MQ=255)
                                                                                                                                                 aaGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCa  <  3:22455/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCA  >  NC_002947/425204‑425496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: