Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,006,594 T→C 8.3% intergenic (+6/+93) ydcJ → / ← PP_5261 metalloprotein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,006,5940TC 8.3% 210.7 / 6.3 72intergenic (+6/+93)ydcJ/PP_5261metalloprotein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base C (2/4);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCGAATACGGCAGCAACGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAA  >  NC_002947/6006445‑6006669
                                                                                                                                                     |                                                                           
gCCGAATACGGCAGCAACGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGt                                                                             <  2:1918339/150‑1 (MQ=255)
  cGAATACGGCAGCAACGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTAC                                                                           >  2:2076851‑M1/1‑150 (MQ=255)
     aTACGGCAGCAACGCCAACCGCGAAGCGTTCGAAGCGGCGGTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGG                                                                        >  1:159038‑M1/1‑150 (MQ=255)
            atcaacGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTC                                                                 <  1:350371‑M1/148‑1 (MQ=255)
               aaCGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCG                                                              <  2:1621473‑M1/150‑1 (MQ=255)
               aaCGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCG                                                              <  2:1297878‑M1/150‑1 (MQ=255)
                  gCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGC                                                           <  1:1145836‑M1/150‑1 (MQ=255)
                    cAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTC                                                         <  2:1021497‑M1/150‑1 (MQ=255)
                     aaCCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCG                                                              <  2:662357‑M1/144‑1 (MQ=255)
                     aaCCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCG                                                              >  1:662357‑M1/1‑144 (MQ=255)
                       ccGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCC                                                      <  2:505088‑M1/150‑1 (MQ=255)
                       ccGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACC                                                                          <  1:1052828‑M1/130‑1 (MQ=255)
                       ccGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACC                                                                          >  2:1052828‑M1/1‑130 (MQ=255)
                       ccGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCC                                                      >  2:1693729‑M1/1‑150 (MQ=255)
                        cgcgAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCC                                                     <  2:1250008‑M1/150‑1 (MQ=255)
                              gCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCC                                               <  1:2348277‑M1/150‑1 (MQ=255)
                               cGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                              <  1:820200‑M1/150‑1 (MQ=255)
                               cGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                              <  1:1724733‑M1/150‑1 (MQ=255)
                                  tCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAG                                           <  2:988589‑M1/150‑1 (MQ=255)
                                  tCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAG                                           <  1:521908‑M1/150‑1 (MQ=255)
                                        cggcgCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAg                                     <  1:2076851/150‑1 (MQ=255)
                                         ggcgCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                    <  1:2235840/150‑1 (MQ=255)
                                         ggcgCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                    <  1:1389113/150‑1 (MQ=255)
                                          gcgcTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAggg                                   >  1:1431849/1‑150 (MQ=255)
                                              tGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCGTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGG                                                             <  1:977412‑M1/120‑1 (MQ=255)
                                              tGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGG                                                             >  2:977412‑M1/1‑120 (MQ=255)
                                               ggggTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctg                              <  1:584695‑M1/150‑12 (MQ=255)
                                                 ggTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgta                            <  2:1841726‑M1/150‑14 (MQ=255)
                                                      cAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgt                       >  1:1144650‑M2/1‑144 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCTGCTCCCACAGGTcgtatgccctgtacctgtg                      <  1:1239554‑M1/150‑20 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtg                      >  2:1161985‑M2/1‑143 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtg                      >  2:159038‑M2/1‑143 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtg                      >  2:248756‑M2/1‑143 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCC                                                >  1:1684233‑M1/1‑124 (MQ=255)
                                                       aGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCC                                                <  2:1684233‑M1/124‑1 (MQ=255)
                                                        ggTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgg                     <  2:2061440‑M1/150‑21 (MQ=255)
                                                        ggTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgg                     >  1:464688‑M2/1‑142 (MQ=255)
                                                        ggTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgg                     >  2:1712290‑M2/1‑142 (MQ=255)
                                                        ggTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgg                     <  1:1693729‑M1/150‑21 (MQ=255)
                                                         gTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                    <  2:894469/134‑1 (MQ=255)
                                                         gTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                    >  1:894469/1‑134 (MQ=255)
                                                              ggATTAGCTGTCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgg               >  1:2279685‑M2/1‑136 (MQ=255)
                                                              ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgg               >  1:367001‑M2/1‑136 (MQ=255)
                                                              ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtggga                   >  1:875772‑M2/1‑136 (MQ=255)
                                                              ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtggga                   <  2:875772‑M2/146‑11 (MQ=255)
                                                              ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcgg               <  1:345125‑M1/150‑27 (MQ=255)
                                                               gATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcggg              >  1:771127‑M1/1‑123 (MQ=255)
                                                               gATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcggg              <  2:771127‑M1/150‑28 (MQ=255)
                                                                aTGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGC                                                       >  2:210465‑M1/1‑108 (MQ=255)
                                                                aTGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGC                                                       <  1:210465‑M1/108‑1 (MQ=255)
                                                                  gAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcga           >  2:1431849‑M2/1‑132 (MQ=255)
                                                                         cgcTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCT                                                                     >  2:274585‑M1/1‑85 (MQ=255)
                                                                         cgcTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCT                                                                     <  1:274585‑M1/85‑1 (MQ=255)
                                                                                  cgcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcgg               <  2:122740‑M1/130‑27 (MQ=255)
                                                                                  cgcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcgg               >  1:122740‑M1/1‑104 (MQ=255)
                                                                                   gcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCC‑‑‑‑‑‑TACAGGGCATACGacctgtggga                   <  2:1072266‑M2/119‑11 (MQ=255)
                                                                                   gcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCC‑‑‑‑‑‑TACAGGGCATACGacctgtggga                   >  1:1072266‑M2/1‑109 (MQ=255)
                                                                                   gcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaa                                                >  1:2239454‑M1/1‑93 (MQ=255)
                                                                                   gcAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaa                                                <  2:2239454‑M1/96‑4 (MQ=255)
                                                                                              gcCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcgagcccgcgaa  <  2:473656‑M2/131‑28 (MQ=255)
                                                                                              gcCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcgagcccgcgaa  >  1:473656‑M2/1‑104 (MQ=255)
                                                                                                  gTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgt                       <  1:357203‑M1/106‑19 (MQ=255)
                                                                                                  gTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGCACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgt                       >  2:357203‑M1/1‑88 (MQ=255)
                                                                                                       tGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagc                 >  1:1133338‑M1/1‑83 (MQ=255)
                                                                                                       tGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagc                 <  2:1133338‑M1/107‑25 (MQ=255)
                                                                                                                    acaGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggc             >  2:512353‑M2/1‑82 (MQ=255)
                                                                                                                    acaGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggc             <  1:512353‑M2/98‑17 (MQ=255)
                                                                                                                                 cTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                              <  2:187152‑M1/52‑1 (MQ=255)
                                                                                                                                 cTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                              >  1:187152‑M1/1‑52 (MQ=255)
                                                                                                                                  tGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcgc           >  2:1655025‑M2/1‑68 (MQ=255)
                                                                                                                                  tGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcga           <  1:1655025‑M2/86‑19 (MQ=255)
                                                                                                                                           tGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAG                                           <  2:1675709‑M1/45‑1 (MQ=255)
                                                                                                                                           tGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAG                                           >  1:1675709‑M1/1‑45 (MQ=255)
                                                                                                                                                     |                                                                           
GCCGAATACGGCAGCAACGCCAACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAA  >  NC_002947/6006445‑6006669

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