Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,032 A→G 5.5% intergenic (+121/+68) PP_5440 → / ← PP_t72 hypothetical protein/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,0320AG 5.5% 207.5 / 2.2 72intergenic (+121/+68)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base G (2/2);  total (33/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 = 6.92e-01

TCGTTCCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343890‑344178
                                                                                                                                               |                                                                                                                                                   
tCGTTCCTTGTCGTACTATTAAGGTTTGGCCCGAATTATGTTAATGGCTGCGGGTTTTTCTATTTCATGGTTAGCCTTAAGGTTGCTTCTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAAc                                                                                                                                                 <  1:93233/148‑1 (MQ=255)
      cATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCAc                                                                                                                                           >  4:138445/1‑147 (MQ=255)
       aTGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAg                                                                                                                                        >  4:162985/1‑149 (MQ=255)
          tCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCg                                                                                                                                    >  1:204497/1‑149 (MQ=255)
            ggACGTGTACGGTGTGGCCCGCAGTCTGCTCATGGCTGGGGGTTTCTCGTTTGCATGGCTAGCCTTCTGTTTGGGTTGTCAGGTCGCGTGTTCTTCGACGCCGCCTGAAGTAAAAAATCACAGACAACAAAAACCCCGCACTAGGCgg                                                                                                                                   <  4:222211/148‑1 (MQ=255)
            ggACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCgg                                                                                                                                   <  2:204497/148‑1 (MQ=255)
             gACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGt                                                                                                                                 >  6:248441/1‑149 (MQ=255)
                gAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGtttt                                                                                                                              >  3:257358/1‑148 (MQ=255)
                  gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            <  1:315482/149‑1 (MQ=255)
                  gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            >  5:213270/1‑149 (MQ=255)
                  gTAACGGTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            <  6:51130/149‑1 (MQ=255)
                   tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                           >  8:8778/1‑149 (MQ=255)
                    aaGGTTTGGCCCGCAGTCTGCTGATGGCTGCGGGTTTCTATATTTCATTGATATAATTGAGTTTTCGTTGTTAGTTCTAGTTTTCTTCGATGATTCCTTAATTTAAAAATTACATTCAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                           <  2:132402/148‑1 (MQ=255)
                    aaGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                           <  3:162985/148‑1 (MQ=255)
                     aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGTGTTTCttgtt                                                                                                                         <  7:14435/149‑1 (MQ=255)
                     aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            <  4:42020/146‑1 (MQ=255)
                      ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGTTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttgttg                                                                                                                        <  4:257358/149‑1 (MQ=39)
                       ctttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGGCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTtgttgt                                                                                                                       <  3:266324/148‑1 (MQ=255)
                        tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTc                                                                                                                      >  4:224486/1‑149 (MQ=255)
                        tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTAAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTc                                                                                                                      >  5:160280/1‑149 (MQ=255)
                          tGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTTTGTTGTct                                                                                                                     <  4:162598/148‑2 (MQ=38)
                            gCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCtt                                                                                                                  >  8:280544/1‑149 (MQ=255)
                             cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTc                                                                                                                 <  3:209614/149‑1 (MQ=255)
                              ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCg                                                                                                                >  3:64590/1‑149 (MQ=255)
                               cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCgg                                                                                                               <  2:95992/149‑1 (MQ=255)
                                gCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCggg                                                                                                              <  4:64590/149‑1 (MQ=255)
                                 cAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGt                                                                                                             <  7:8778/149‑1 (MQ=255)
                                      tGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACtt                                                                                                        >  7:119976/1‑149 (MQ=255)
                                        cTCCTGGCTGCGGGGTTCTCGATTGCCTGGCTAGCCTTCAGGTTGCGTTGTCAGGTCTCGTGTTCTACGACGCCGCCTGAAGTAAAAACTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGAAACTTTg                                                                                                      <  8:304703/149‑1 (MQ=255)
                                        cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTg                                                                                                      <  6:227904/149‑1 (MQ=255)
                                        cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTg                                                                                                      >  1:59089/1‑149 (MQ=255)
                                           aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCa                                                                                                    <  1:110964/148‑1 (MQ=255)
                                           aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCa                                                                                                    <  2:59089/148‑1 (MQ=255)
                                               cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                               >  4:308137/1‑149 (MQ=255)
                                               cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                               <  5:248441/149‑1 (MQ=255)
                                                   gggTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAg                                                                                            >  3:314141/1‑148 (MQ=255)
                                                         ctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCGTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGATGTCGCTTCGGGTAACTTTGCAACCACCAGagcctg                                                                                      >  6:94108/1‑143 (MQ=255)
                                                         ctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCAGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAACCtt                                                                                      >  8:122183/1‑148 (MQ=255)
                                                          tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGa                                                                                    <  3:138445/149‑1 (MQ=255)
                                                          tcGATTGCATGGCTAGCCTTCACGTTGCGTTGTGAGGACTCGTGTTCATCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACAACGCGCGTTTCTTCTTGTCGCTTCTGCAAACTTTGCAAGCACCAGAAGCTTg                                                                                     >  5:49040/1‑148 (MQ=255)
                                                            gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  7:296405/148‑1 (MQ=255)
                                                            gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  4:314141/148‑1 (MQ=255)
                                                                gcttggcttggctTGAGGTTGCGTTGTGAGGTGTCGTGTTCTTCGACGGCGCCTGATGTTAAAAATGACAGACAACAAAAAACGCGCTCTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCATCCACCAGAAGGTTGTAggtg                                                                               <  6:49040/145‑1 (MQ=255)
                                                                     gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                          <  2:59579/148‑1 (MQ=255)
                                                                       tAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                          >  7:44262/1‑146 (MQ=255)
                                                                        aGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  5:144142/136‑1 (MQ=255)
                                                                        aGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   >  6:144142/1‑136 (MQ=255)
                                                                            ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                  <  8:119976/149‑1 (MQ=255)
                                                                               aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCCCAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  6:213270/148‑1 (MQ=255)
                                                                               aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  3:224486/148‑1 (MQ=255)
                                                                                 gTTGCGTTGTGAGGTCTCGTGTTCTTCGAAGCCGCCTGAAGTTAAAAATCACATAGAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTTCCCAGAGACGGAGTc                                                             <  4:116111/149‑1 (MQ=255)
                                                                                           tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                     >  1:98620/1‑148 (MQ=255)
                                                                                             aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGac                                                  <  3:308137/149‑1 (MQ=255)
                                                                                                  tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                             >  6:12741/1‑149 (MQ=255)
                                                                                                         cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGTTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                       <  6:160280/148‑1 (MQ=255)
                                                                                                         cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                       >  4:87521/1‑148 (MQ=255)
                                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGa                                                           <  4:303138/127‑1 (MQ=255)
                                                                                                          ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGa                                                           >  3:303138/1‑127 (MQ=255)
                                                                                                                   gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGc                            >  5:312041/1‑149 (MQ=255)
                                                                                                                              aaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  8:183173/83‑1 (MQ=255)
                                                                                                                              aaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   >  7:183173/1‑83 (MQ=255)
                                                                                                                              aaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACt                  >  7:119210/1‑148 (MQ=255)
                                                                                                                               aaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  6:111349/149‑1 (MQ=255)
                                                                                                                               aaaTCAAAGACAACAAAAAACACGCACTAGGCCGGATTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGAACAGAGACCCAGTCGAACCGCCGAAACGAGCATTTTCAATCCTCTGCTCAACCGACTGa                >  8:99236/1‑149 (MQ=255)
                                                                                                                                    acaGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGCCACGAGGATTTTCAATCCTCTTCTCTACCGACTGAGATa            >  2:283294/1‑148 (MQ=255)
                                                                                                                                      aGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGTTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGac                                                  <  7:109773/108‑1 (MQ=255)
                                                                                                                                      aGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGac                                                  >  8:109773/1‑108 (MQ=255)
                                                                                                                                      aGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCt         <  1:305584/149‑1 (MQ=255)
                                                                                                                                          aacaaGAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc     >  8:94647/1‑149 (MQ=255)
                                                                                                                                          aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                      <  5:198964/132‑1 (MQ=255)
                                                                                                                                          aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                      >  6:198964/1‑132 (MQ=255)
                                                                                                                                          aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc     >  4:58227/1‑149 (MQ=255)
                                                                                                                                              aaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac  >  8:209621/1‑148 (MQ=255)
                                                                                                                                               aaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac  >  8:8746/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCGTTCCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343890‑344178

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: