Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,355,046 Δ1 bp 9.1% coding (246/681 nt) modB → molybdate ABC transporter permease ModB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,355,0460G. 9.1% 206.5 / 14.6 66coding (246/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base . (2/4);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

TCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGC  >  NC_002947/4354907‑4355188
                                                                                                                                           |                                                                                                                                               
tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTg                                                                                                                                       >  4:182259/1‑149 (MQ=255)
  gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                      <  7:288610/148‑1 (MQ=255)
  gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTggg                                                                                                                                     <  4:273639/149‑1 (MQ=255)
   ctcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCTCTGGTACTGCCACCGACGGTTATCGGTTTTTACCTGTTGATCGACCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggt                                                                                                                                    <  7:213731/147‑2 (MQ=255)
   cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggt                                                                                                                                    <  2:160403/149‑2 (MQ=255)
   cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGGCAGGCCACCCAGGCCTAGCGCCTgggt                                                                                                                                    <  1:231513/149‑2 (MQ=255)
      gcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgt                                                                                                                                 >  7:160882/1‑149 (MQ=255)
        aCGCGCTCTTGGCTGAGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgtgt                                                                                                                               >  7:54530/1‑149 (MQ=255)
            gcTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggtg                                                                                                                           >  5:139411/1‑149 (MQ=255)
                        gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                              >  2:100436/1‑134 (MQ=255)
                        gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCCGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                              <  1:100436/134‑1 (MQ=255)
                             ccGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCt                                                                                                          >  5:66994/1‑149 (MQ=255)
                                gtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGt                                                                                                       >  3:199600/1‑149 (MQ=255)
                                  gggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGa                                                                                                     <  8:160882/149‑1 (MQ=255)
                                   ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAt                                                                                                    <  6:91333/149‑1 (MQ=255)
                                   cccgatctggtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGt                                                                                                       <  7:72300/139‑1 (MQ=255)
                                        ggtggtGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                               <  3:182259/149‑1 (MQ=255)
                                            gtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGt                                                                                                       >  8:72300/1‑137 (MQ=255)
                                             tggCACTGCCTCTGGTACTGCGTCCGTCGGTGATGGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCTCAGGCCTGGTGATTGGCTCGACg                                                                                          <  4:199600/149‑1 (MQ=255)
                                             tggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt                                                                                                                            <  7:109704/115‑1 (MQ=255)
                                             tggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt                                                                                                                            >  8:109704/1‑115 (MQ=255)
                                              ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACgg                                                                                         >  7:117095/1‑149 (MQ=255)
                                               gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                        <  6:139411/149‑1 (MQ=255)
                                                cACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                       <  5:190308/149‑1 (MQ=255)
                                                 aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                       >  1:15034/1‑148 (MQ=255)
                                                     ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTc                                                                                  <  8:117095/149‑1 (MQ=255)
                                                                  cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                     >  8:246115/1‑149 (MQ=255)
                                                                  cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                     <  2:15034/149‑1 (MQ=255)
                                                                    ccGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                   >  6:163944/1‑149 (MQ=255)
                                                                     cGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                  <  7:246115/149‑1 (MQ=255)
                                                                      gACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGc                                                                 <  6:66994/149‑1 (MQ=255)
                                                                             aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCt                                                          >  6:230992/1‑149 (MQ=255)
                                                                              tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTg                                                         >  4:44484/1‑149 (MQ=255)
                                                                                 gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                               >  5:26846/1‑140 (MQ=255)
                                                                                 gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                               <  6:26846/140‑1 (MQ=255)
                                                                                 gTTTTTACCTGTTGATAGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGc                                                                                                            >  6:98948/1‑95 (MQ=255)
                                                                                 gTTTTTACCTGTTGATAGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGc                                                                                                            <  5:98948/95‑1 (MQ=255)
                                                                                         cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCtt                                              >  3:203006/1‑149 (MQ=255)
                                                                                         cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACACAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACACCTGCAGAATGCCtt                                              >  5:240869/1‑149 (MQ=255)
                                                                                            ttGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTgg                                           >  7:264707/1‑149 (MQ=255)
                                                                                               aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcg                                         <  3:44484/148‑1 (MQ=255)
                                                                                                  gCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGAGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGACGTTCGTGGTGCAACCCCTGCAGAATGTCTTTGGCgaga                                      >  6:10715/1‑148 (MQ=255)
                                                                                                    ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCACCCCCTGCAGAATGCCTTTGGCGCGATCg                                   >  5:260055/1‑149 (MQ=255)
                                                                                                       tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGcc                                                >  2:276990/1‑133 (MQ=255)
                                                                                                       tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGcc                                                <  1:276990/133‑1 (MQ=255)
                                                                                                        ggCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCa                               >  4:74338/1‑149 (MQ=255)
                                                                                                          cccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgc                             <  5:230992/149‑1 (MQ=255)
                                                                                                           ccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcg                            >  7:93440/1‑149 (MQ=255)
                                                                                                                  gctggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTg                     >  3:252766/1‑149 (MQ=255)
                                                                                                                   ctggctTGGTCAGGCCACCCAGGCGCTATGCCTGGGGAGTGTGGTGTTCAGCTTCACAGGCCTGTTGATTGGCTCGAGGGTGTACTCCATGCCGTTCGTGGTGCAAGCCCTGCAGAATGTCTTTGGCGCGATCGGCCAGCGCCCCCTgg                    <  5:10715/149‑1 (MQ=255)
                                                                                                                    tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                             >  7:35744/1‑123 (MQ=255)
                                                                                                                    tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                             <  8:35744/123‑1 (MQ=255)
                                                                                                                     ggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                  <  5:88979/149‑1 (MQ=255)
                                                                                                                      gctTGGTCAGGCCAGCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAg                 <  8:54530/149‑1 (MQ=255)
                                                                                                                      gctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                  <  5:163944/148‑1 (MQ=255)
                                                                                                                        ttGGTCAGTCCACCCATGCGCTAGGCCTGGGCATTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACTGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                             <  2:212562/119‑1 (MQ=255)
                                                                                                                        ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                             >  1:212562/1‑119 (MQ=255)
                                                                                                                           gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgcc            <  3:74338/149‑1 (MQ=255)
                                                                                                                           gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcg              >  8:263799/1‑147 (MQ=255)
                                                                                                                            tCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg           >  1:142051/1‑149 (MQ=255)
                                                                                                                             cAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccgc          <  8:264707/149‑1 (MQ=255)
                                                                                                                                gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGCCGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCTGGAAGTCGCCGCTAct      >  7:237246/1‑149 (MQ=255)
                                                                                                                                actgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAc       >  2:10435/5‑149 (MQ=255)
                                                                                                                                  tgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCAGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGACTGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActc     >  1:214159/3‑149 (MQ=255)
                                                                                                                                        ggCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGc  >  6:216872/1‑146 (MQ=255)
                                                                                                                                        ggCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGc  <  5:216872/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
TCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGC  >  NC_002947/4354907‑4355188

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: