Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,355,050:1 +C 9.0% coding (250/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,0501.C 9.0% 252.1 / 17.6 67coding (250/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base . (33/28);  new base C (2/4);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGCCAG  >  NC_002947/4354907‑4355195
                                                                                                                                                |                                                                                                                                                  
tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTg                                                                                                                                               >  4:182259/1‑149 (MQ=255)
  gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTgg                                                                                                                                              <  7:288610/148‑1 (MQ=255)
  gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTggg                                                                                                                                             <  4:273639/149‑1 (MQ=255)
   ctcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCTCTGGTACTGCCACCGACGGTTATCGGTTTTTACCTGTTGATCGACCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggt                                                                                                                                            <  7:213731/147‑2 (MQ=255)
   cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggt                                                                                                                                            <  2:160403/149‑2 (MQ=255)
   cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGGCAGGCCACCCAGGC‑CTAGCGCCTgggt                                                                                                                                            <  1:231513/149‑2 (MQ=255)
      gcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAgt                                                                                                                                         >  7:160882/1‑149 (MQ=255)
        aCGCGCTCTTGGCTGAGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAgtgt                                                                                                                                       >  7:54530/1‑149 (MQ=255)
            gcTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTgtggtg                                                                                                                                   >  5:139411/1‑149 (MQ=255)
                        gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGtgtg                                                                                                                                      >  2:100436/1‑134 (MQ=255)
                        gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCCGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGtgtg                                                                                                                                      <  1:100436/134‑1 (MQ=255)
                             ccGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑t                                                                                                                 >  5:66994/1‑149 (MQ=255)
                                gtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGt                                                                                                              >  3:199600/1‑149 (MQ=255)
                                  gggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGa                                                                                                            <  8:160882/149‑1 (MQ=255)
                                   ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGAt                                                                                                           <  6:91333/149‑1 (MQ=255)
                                   cccgatctggtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGt                                                                                                              <  7:72300/139‑1 (MQ=255)
                                        ggtggtGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCt                                                                                                      <  3:182259/149‑1 (MQ=255)
                                            gtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGt                                                                                                              >  8:72300/1‑137 (MQ=255)
                                             tggCACTGCCTCTGGTACTGCGTCCGTCGGTGATGGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCTCAGGCC‑TGGTGATTGGCTCGACg                                                                                                 <  4:199600/149‑1 (MQ=255)
                                             tggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTgtggt                                                                                                                                    >  8:109704/1‑115 (MQ=255)
                                             tggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTgtggt                                                                                                                                    <  7:109704/115‑1 (MQ=255)
                                              ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACgg                                                                                                >  7:117095/1‑149 (MQ=255)
                                               gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGgt                                                                                               <  6:139411/149‑1 (MQ=255)
                                                cACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGgtg                                                                                              <  5:190308/149‑1 (MQ=255)
                                                 aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGgtg                                                                                              >  1:15034/1‑148 (MQ=255)
                                                     ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTc                                                                                         <  8:117095/149‑1 (MQ=255)
                                                                  cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                            <  2:15034/149‑1 (MQ=255)
                                                                  cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                            >  8:246115/1‑149 (MQ=255)
                                                                    ccGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                          >  6:163944/1‑149 (MQ=255)
                                                                     cGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                         <  7:246115/149‑1 (MQ=255)
                                                                      gACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGc                                                                        <  6:66994/149‑1 (MQ=255)
                                                                             aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCt                                                                 >  6:230992/1‑149 (MQ=255)
                                                                              tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTg                                                                >  4:44484/1‑149 (MQ=255)
                                                                                 gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                                      <  6:26846/140‑1 (MQ=255)
                                                                                 gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                                      >  5:26846/1‑140 (MQ=255)
                                                                                 gTTTTTACCTGTTGATAGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGc                                                                                                                    <  5:98948/95‑1 (MQ=255)
                                                                                 gTTTTTACCTGTTGATAGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGc                                                                                                                    >  6:98948/1‑95 (MQ=255)
                                                                                         cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCtt                                                     >  3:203006/1‑149 (MQ=255)
                                                                                         cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACACAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACACCTGCAGAATGCCtt                                                     >  5:240869/1‑149 (MQ=255)
                                                                                            ttGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTgg                                                  >  7:264707/1‑149 (MQ=255)
                                                                                               aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcg                                                <  3:44484/148‑1 (MQ=255)
                                                                                                  gCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGAGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGACGTTCGTGGTGCAACCCCTGCAGAATGTCTTTGGCgaga                                             >  6:10715/1‑148 (MQ=255)
                                                                                                    ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCACCCCCTGCAGAATGCCTTTGGCGCGATCg                                          >  5:260055/1‑149 (MQ=255)
                                                                                                       tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGcc                                                       <  1:276990/133‑1 (MQ=255)
                                                                                                       tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGcc                                                       >  2:276990/1‑133 (MQ=255)
                                                                                                        ggCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCa                                      >  4:74338/1‑149 (MQ=255)
                                                                                                          cccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgc                                    <  5:230992/149‑1 (MQ=255)
                                                                                                           ccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcg                                   >  7:93440/1‑149 (MQ=255)
                                                                                                                  gctggctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTg                            >  3:252766/1‑149 (MQ=255)
                                                                                                                   ctggctTGGTCAGGCCACCCAGGCGCTAT‑GCCTGGGGAGTGTGGTGTTCAGCTTCACAGGCC‑TGTTGATTGGCTCGAGGGTGTACTCCATGCCGTTCGTGGTGCAAGCCCTGCAGAATGTCTTTGGCGCGATCGGCCAGCGCCCCCTgg                           <  5:10715/149‑1 (MQ=255)
                                                                                                                    tggctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                    >  7:35744/1‑123 (MQ=255)
                                                                                                                    tggctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                    <  8:35744/123‑1 (MQ=255)
                                                                                                                     ggctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                         <  5:88979/149‑1 (MQ=255)
                                                                                                                      gctTGGTCAGGCCAGCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAg                        <  8:54530/149‑1 (MQ=255)
                                                                                                                      gctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                         <  5:163944/148‑1 (MQ=255)
                                                                                                                        ttGGTCAGTCCACCCATGCGCTAG‑GCCTGGGCATTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACTGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                    <  2:212562/119‑1 (MQ=255)
                                                                                                                        ttGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                    >  1:212562/1‑119 (MQ=255)
                                                                                                                           gTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgcc                   <  3:74338/149‑1 (MQ=255)
                                                                                                                           gTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcg                     >  8:263799/1‑147 (MQ=255)
                                                                                                                            tCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg                  >  1:142051/1‑149 (MQ=255)
                                                                                                                             cAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccgc                 <  8:264707/149‑1 (MQ=255)
                                                                                                                                gCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGCCGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCG‑CCCCTGGAAGTCGCCGCTAct             >  7:237246/1‑149 (MQ=255)
                                                                                                                                actgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAc              >  2:10435/5‑149 (MQ=255)
                                                                                                                                  tgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑AGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGACTGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActc            >  1:214159/3‑149 (MQ=255)
                                                                                                                                        ggCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGc         >  6:216872/1‑146 (MQ=255)
                                                                                                                                        ggCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGc         <  5:216872/146‑1 (MQ=255)
                                                                                                                                             tAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGagccag  >  1:216115/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCC‑TGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGCCAG  >  NC_002947/4354907‑4355195

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: ‑