Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,708,585 A→G 7.1% W223R (TGG→CGG)  cysI ← sulfite reductase hemoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,708,5850AG 7.1% 280.2 / 9.0 98W223R (TGG→CGG) cysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base A (42/49);  new base G (4/3);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 101 total aligned reads displayed.

CCATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGC  >  NC_002947/2708438‑2708733
                                                                                                                                                   |                                                                                                                                                      
ccATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACg                                                                                                                                                      >  2:2344236/1‑150 (MQ=255)
 cATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCaagc                                                                                                                                                     <  1:695904/150‑4 (MQ=255)
      aggCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGa                                                                                                                                                >  1:1509087/2‑150 (MQ=255)
       ggCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGTTCCTGCCGTGGCAGGac                                                                                                                                               <  2:237256/150‑2 (MQ=255)
       ggCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGac                                                                                                                                               <  1:806159/150‑2 (MQ=255)
             gACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         <  2:1694829/150‑1 (MQ=255)
                   ttCGGCGAACACTTCCGGGCTAAGGGCCTTTGCCAGGATCTTGATCCGAGCCTTGTACTGGTGGACATCGCCTTTGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGaa                                                                                                                                   <  1:43762/150‑1 (MQ=255)
                    tCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTGGCGGGGGCCGTAACGGTTGTAGACGCGCAGGATGGCGTCGAGGTAGCTGGTCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑c                                                                                                                                <  1:318448/150‑1 (MQ=255)
                    tCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑c                                                                                                                                <  2:1173558/150‑1 (MQ=255)
                          acaacTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCCTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACCGCAGGAACTCGCTGATGAA‑‑CGAAccc                                                                                                                          >  2:2402316/4‑150 (MQ=255)
                             acTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAAAGGTTGTACACACGCATGTTAGCGGCTTGGTAGCTGATCATGTCCGTGCACGGCATGGACTCGGTGATGAA‑‑CGAGCccaac                                                                                                                       >  2:1109527/1‑148 (MQ=255)
                               ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         >  2:1211601/1‑132 (MQ=255)
                               ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         <  1:1211601/132‑1 (MQ=255)
                               ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTtgat                                                                                                                                      <  2:290990/135‑1 (MQ=255)
                               ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTtgat                                                                                                                                      >  1:290990/1‑135 (MQ=255)
                                tCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCaccacc                                                                                                                    >  2:343213/1‑150 (MQ=255)
                                  cGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCgg                                                                                                                  <  1:1109527/150‑1 (MQ=255)
                                    gggTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCgggg                                                                                                                >  1:243535/1‑150 (MQ=255)
                                      gTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCg                                                                                                                   >  2:1155832/1‑145 (MQ=255)
                                      gTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCg                                                                                                                   <  1:1155832/145‑1 (MQ=255)
                                        aaGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACg                                                                                                            <  1:2344236/150‑1 (MQ=255)
                                         aGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACgg                                                                                                           <  2:1638243/150‑1 (MQ=255)
                                              cTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAg                                                                                                      >  1:1101279/1‑150 (MQ=255)
                                                tGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATTGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGc                                                                                                    <  1:867195/150‑1 (MQ=255)
                                                  aCCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGcc                                                                                                                           >  1:2295188/1‑125 (MQ=255)
                                                  aCCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGcc                                                                                                                           <  2:2295188/125‑1 (MQ=255)
                                                   ccAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGccac                                                                                                 <  2:68361/150‑1 (MQ=255)
                                                   ccAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGccac                                                                                                 >  1:620641/1‑150 (MQ=255)
                                                    cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGccacc                                                                                                >  1:1317114/1‑150 (MQ=255)
                                                      ggATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAg                                                                                                                                                  <  2:1266188/100‑1 (MQ=255)
                                                      ggATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAg                                                                                                                                                  >  1:1266188/1‑100 (MQ=255)
                                                       gATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTCCTTGCTCAATTACAGCGACGGCAGTAACTCGTTGATAAACCCCAGCCCAATAC‑‑GGGTATGGGCCAGGCCaccacc                                                                                             >  1:1853415/1‑150 (MQ=255)
                                                       gATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCaccacc                                                                                                                    >  1:1448215/1‑127 (MQ=255)
                                                       gATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCaccacc                                                                                                                    <  2:1448215/127‑1 (MQ=255)
                                                            tGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGTTCCTGCCGTGGCAGGACCtgat                                                                                                                                          >  2:2249993/1‑99 (MQ=255)
                                                            tGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACg                                                                                                            >  2:1223107/1‑130 (MQ=255)
                                                            tGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACg                                                                                                            <  1:1223107/130‑1 (MQ=255)
                                                                  gcgcCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCacac                                                                                  >  2:1266810/1‑150 (MQ=255)
                                                                    gcCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACgg                                                                                <  2:1509087/150‑1 (MQ=255)
                                                                      cTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGaa                                                                              <  2:1317114/150‑1 (MQ=255)
                                                                      cTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGaa                                                                              <  1:343213/150‑1 (MQ=255)
                                                                       ttGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGaa                                                                                                                                               >  2:766853/1‑86 (MQ=255)
                                                                       ttGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGaa                                                                                                                                               <  1:766853/86‑1 (MQ=255)
                                                                         gTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                           <  2:1101279/150‑1 (MQ=255)
                                                                           gttggttgtCACGGCGACGGTAAAGATTGTACGCGCGCAGGATGGGGGCGAGGTAGCTGATCAGGGCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCa                                                                         <  2:463973/146‑1 (MQ=255)
                                                                            cTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTAc                                                                                                             <  2:1786093/113‑1 (MQ=255)
                                                                            cTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTAc                                                                                                             >  1:1786093/1‑113 (MQ=255)
                                                                             ttgttgTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATTAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGt                                                                       >  2:763176/1‑150 (MQ=255)
                                                                             ttgttgTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGt                                                                       >  2:596938/1‑150 (MQ=255)
                                                                                ttgtCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCg                                                                                                                   <  2:1548878/103‑1 (MQ=255)
                                                                                ttgtCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCg                                                                                                                   >  1:1548878/1‑103 (MQ=255)
                                                                                ttGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCG                                                                    >  1:1769175‑M1/1‑150 (MQ=255)
                                                                                ttGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGAACGGCCCCGGCCACCACCTACCAGCACACGGAAGCCCAGTTCG                                                                    >  2:445639‑M1/1‑150 (MQ=255)
                                                                                  gTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGCAGCCCAGTTCGCC                                                                  <  2:2333115‑M1/150‑1 (MQ=255)
                                                                                  gTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCC                                                                  <  1:596938‑M1/150‑1 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑c                                                                                                                                <  2:2188500/87‑1 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑c                                                                                                                                >  1:2188500/1‑87 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCggg                                                                                                                 <  1:2175616/102‑1 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCggg                                                                                                                 >  2:2175616/1‑102 (MQ=255)
                                                                                     aCGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         >  2:515200/1‑78 (MQ=255)
                                                                                     aCGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         <  1:515200/78‑1 (MQ=255)
                                                                                       ggCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                           >  2:777688/1‑136 (MQ=255)
                                                                                       ggCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                           <  1:777688/136‑1 (MQ=255)
                                                                                        gCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCAACGGCAGGGACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGc                                                            >  1:1479229/1‑150 (MQ=255)
                                                                                              gTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgcg                                                      >  1:839908/1‑150 (MQ=255)
                                                                                                   ggTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   <  2:1068867/54‑1 (MQ=255)
                                                                                                   ggTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   >  1:1068867/1‑54 (MQ=255)
                                                                                                    gTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCAccgg                                                <  2:1742257‑M2/150‑5 (MQ=255)
                                                                                                    gTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCAccgg                                                <  1:2183640‑M2/150‑5 (MQ=255)
                                                                                                     ttGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCcaccac                                                                                              >  2:1230449/1‑103 (MQ=255)
                                                                                                     ttGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCcaccac                                                                                              <  1:1230449/103‑1 (MQ=255)
                                                                                                       gTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGtgcgcaacgctgccggcga                                             <  1:2025923‑M1/150‑20 (MQ=255)
                                                                                                       gTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTc                                             <  1:1742257/150‑1 (MQ=255)
                                                                                                       gTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTc                                             <  1:2402316/150‑1 (MQ=255)
                                                                                                       gTACACGCGCAGGATGGCGTCGAGGTAGCCGATCAGGTCCTGCCACGTCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTc                                             <  2:127145/150‑1 (MQ=255)
                                                                                                          cacGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAg                                          <  2:243535/150‑1 (MQ=255)
                                                                                                             gcgcAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAg                                                                                                                             >  2:525487/1‑64 (MQ=255)
                                                                                                             gcgcAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAg                                                                                                                             <  1:525487/64‑1 (MQ=255)
                                                                                                                     tGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGAAGGAACTCGTTGATGCA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGt                               <  2:271691/150‑1 (MQ=255)
                                                                                                                        cGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCa                            >  1:947983/1‑150 (MQ=255)
                                                                                                                                 aGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCCGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATgg                   <  1:1266810/150‑1 (MQ=255)
                                                                                                                                    tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAg                                                                                                      >  2:1110505/1‑64 (MQ=255)
                                                                                                                                    tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAg                                                                                                      <  1:1110505/64‑1 (MQ=255)
                                                                                                                                       gcAGGTCCTGCGACGGCAGGAACTGGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGgcgc             >  1:953634/2‑150 (MQ=255)
                                                                                                                                         aGGGCCTGCCACGGCAGGAACTCGTTGACGAA‑‑CGAGCCCACCAACGGGGTACGGACCAGGCCGACACCTACCACCAAACGGACGCGCAGGTGTCCGGCAGGGTTCCGGACGGGCTGAACGCCAATGTCGTGCACTTCGATGGCGGCGCgg           <  2:1479229/150‑1 (MQ=255)
                                                                                                                                          gttCCTGCCGTGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTg                              >  1:1828055/3‑130 (MQ=255)
                                                                                                                                          gttCCTGCCGTGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTg                              <  2:1828055/128‑1 (MQ=255)
                                                                                                                                          gttCCTGCCGTGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGATGAAGTCCAATGTCGTGCACTTCGATGGCGGCGCGGt          >  2:2371189/3‑150 (MQ=255)
                                                                                                                                           ttCCTGCCGTGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTc         >  2:1246749/2‑150 (MQ=255)
                                                                                                                                            tCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGc                                                                                                    <  2:2250500/58‑1 (MQ=255)
                                                                                                                                            tCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGc                                                                                                    <  2:2250003/58‑1 (MQ=255)
                                                                                                                                            tCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGc                                                                                                    >  1:2250003/1‑58 (MQ=255)
                                                                                                                                            tCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGc                                                                                                    >  1:2250500/1‑58 (MQ=255)
                                                                                                                                              cTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc      >  2:838136/1‑150 (MQ=255)
                                                                                                                                              cTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc      <  2:620641/150‑1 (MQ=255)
                                                                                                                                                gCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCt    <  2:2183640/150‑1 (MQ=255)
                                                                                                                                                 ccACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGCTCGCCGGCAGCGTTGCGCACCGGCTCAAGCCCAATGTCGTGCACTTCGATGGCGGCGCGGTGTTCCTg   >  2:30244/1‑150 (MQ=255)
                                                                                                                                                 ccACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCAACCACGTACCAGCCCACGGAAGCCCAGTTCGCCGCCAGCGTTGCGCACCAACTCAAGGCCAACGTCGCGGAACTCGATGGCGGGGCGGTGTTGCTg   >  2:113489/1‑150 (MQ=255)
                                                                                                                                                   aCGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGc  >  1:1849202/1‑149 (MQ=255)
                                                                                                                                                   aCGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGc  <  2:1849202/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAA‑‑CGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGC  >  NC_002947/2708438‑2708733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: