Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,561:1 +G 8.5% coding (372/540 nt) moaB‑I → molybdenum cofactor biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,5611.G 8.5% 222.2 / 14.6 59coding (372/540 nt)moaB‑Imolybdenum cofactor biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base . (28/26);  new base G (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GATTACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGG  >  NC_002947/2421414‑2421709
                                                                                                                                                    |                                                                                                                                                     
gATTACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGccg                                                                                                                                                       >  6:34388/1‑148 (MQ=255)
gATTACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCt                                                                                                                                                      >  2:375498/1‑149 (MQ=255)
 aTTACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGGTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGccgt                                                                                                                                                      >  2:302963/1‑147 (MQ=255)
  ttACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTg                                                                                                                                                    <  1:147475/149‑1 (MQ=255)
     ccGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCTCCTCCACCGTGCAGTCGCGTGCCTTGGcc                                                                                                                                                 <  4:54324/149‑1 (MQ=255)
     ccGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGcc                                                                                                                                                 >  6:210384/1‑149 (MQ=255)
      cggcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCg                                                                                                                                                <  6:14606/149‑1 (MQ=255)
      cggcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCg                                                                                                                                                <  8:189661/149‑1 (MQ=255)
        gcggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCggg                                                                                                                                              >  6:85495/1‑149 (MQ=255)
         cggTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGc                                                                                                                                             <  6:179096/149‑1 (MQ=255)
           gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                           <  4:199906/149‑1 (MQ=255)
            tACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTgg                                                                                                                                          <  6:21826/149‑1 (MQ=255)
                     tACTGGTCGCGACAGCACCCCCGCGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggc                                                                                                                                   <  5:34388/148‑4 (MQ=255)
                        tGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGc                                                                                                                              >  3:114647/1‑149 (MQ=255)
                          gTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCAAGGCACGcga                                                                                                                            <  1:302963/149‑3 (MQ=255)
                            cgcgACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGgt                                                                                                                          <  1:344991/149‑1 (MQ=255)
                             gcgACAGCACCCCCGAGGCGGTGGACTGCCTGCTGGACAATCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGAACggttgtt                                                                                                                         >  7:57737/1‑143 (MQ=255)
                               gACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCTAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGtgtg                                                                                                                       >  5:374312/1‑149 (MQ=255)
                                     aCCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcct                                                                                                                 >  7:5216/1‑149 (MQ=255)
                                       ccccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcctgc                                                                                                               >  1:326330/1‑149 (MQ=255)
                                        cccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaag                                                                                                                                     >  4:158030/1‑126 (MQ=255)
                                         ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCa                                                                                                             <  8:55818/149‑1 (MQ=255)
                                         ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAAGGGCCGCTGGTGTGctgcctgc                                                                                                               >  3:195895/1‑147 (MQ=255)
                                                ggTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTcca                                                                                                      >  5:213650/1‑149 (MQ=255)
                                                   ggCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGGCCGCTGGTGTGGTGCCTGCCCGGCTccacc                                                                                                    >  1:383089/1‑148 (MQ=255)
                                                      ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                  <  1:209379/116‑5 (MQ=38)
                                                      ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                  >  2:209379/1‑112 (MQ=38)
                                                        gcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCt                                                                                              >  4:231679/1‑149 (MQ=255)
                                                            gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                          >  3:108052/1‑149 (MQ=255)
                                                            gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                          >  4:98209/1‑149 (MQ=255)
                                                                gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                           >  1:27396/1‑96 (MQ=255)
                                                                gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                           <  2:27396/96‑1 (MQ=255)
                                                                gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCATGCTCCACCAACGCTGTGCGTAc                                                                                      >  2:239925/1‑149 (MQ=255)
                                                                gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                      >  4:81055/1‑149 (MQ=255)
                                                                gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                      >  1:231193/1‑149 (MQ=255)
                                                                     gcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGtt                                                                                 >  3:311270/1‑149 (MQ=255)
                                                                         gTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGa                                                                             >  8:222836/1‑149 (MQ=255)
                                                                            gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGGCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                           <  5:85495/148‑1 (MQ=255)
                                                                            gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                           <  7:222836/148‑1 (MQ=255)
                                                                                gCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATc                                                                      <  2:326330/149‑1 (MQ=255)
                                                                                  tttGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCt                                                                    >  5:93215/1‑149 (MQ=255)
                                                                                      gTGAGCTGTCCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGc                                                                 <  6:213650/148‑1 (MQ=255)
                                                                                       tGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCg                                                               <  4:108052/149‑1 (MQ=255)
                                                                                         aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGag                                                             <  5:210384/149‑1 (MQ=255)
                                                                                         aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGa                                                              <  4:124599/148‑1 (MQ=255)
                                                                                              ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGagcagc                                                         <  2:385297/148‑1 (MQ=255)
                                                                                              ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTGTGCCGGGCGGGCCAAACGCACCCGGGTGTGCTTCCAGCCCCGCTCCACCAACGCTGTGGGTTCCGGGTGGGGTTGCCTGCTGGCCGCGCCGCt                                                        >  1:300476/1‑149 (MQ=255)
                                                                                                ccGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTc                                                       <  3:231679/148‑1 (MQ=255)
                                                                                                 cGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGa                                                     >  3:67840/1‑149 (MQ=255)
                                                                                                   gCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATg                                                  <  6:374312/149‑1 (MQ=255)
                                                                                                      gATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGTCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                               <  4:311270/149‑1 (MQ=255)
                                                                                                        tttCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGt                                             <  2:130834/149‑1 (MQ=255)
                                                                                                          tCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTc                                            <  4:114647/148‑1 (MQ=255)
                                                                                                           ccGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCAt                                          >  3:315297/1‑149 (MQ=255)
                                                                                                                 ccGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCAGCCATGCTCCACCAACCCTTTGCGTACCGGTTGGGATTGCATCCAGGCCGACAAGCACGAAGCCCGTCATCGccca                                    >  7:232975/1‑148 (MQ=255)
                                                                                                                          gCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                           <  4:195895/149‑1 (MQ=255)
                                                                                                                              cTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTgg                       <  1:239925/149‑1 (MQ=255)
                                                                                                                                 cACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGgcg                     >  6:327271/1‑148 (MQ=255)
                                                                                                                                     gTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCt                >  5:35296/1‑149 (MQ=255)
                                                                                                                                        cAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGGTCCACCAACGCTGTGGGTTCCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTg               >  6:287393/1‑147 (MQ=255)
                                                                                                                                                tGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgcaggcag     >  1:359623/1‑149 (MQ=255)
                                                                                                                                                 gCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc    <  3:81055/149‑1 (MQ=255)
                                                                                                                                                  ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATCGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc    <  1:280829/149‑1 (MQ=255)
                                                                                                                                                  ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg   <  2:231193/149‑1 (MQ=255)
                                                                                                                                                  ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg   <  5:327271/149‑1 (MQ=255)
                                                                                                                                                   cTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg  >  6:134092/1‑149 (MQ=255)
                                                                                                                                                   cTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg  >  3:116310/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTACCGGCGGTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGG  >  NC_002947/2421414‑2421709

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