Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,877 G→A 5.5% D137N (GAC→AAC)  soxR → redox‑sensitive transcriptional activator of soxS; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,8770GA 5.5% 238.9 / 3.0 73D137N (GAC→AAC) soxRredox‑sensitive transcriptional activator of soxS; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base A (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACAC  >  NC_000913/4277739‑4278024
                                                                                                                                          |                                                                                                                                                   
tGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAg                                                                                                                                            <  6:242360/148‑1 (MQ=255)
         cTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGaagaag                                                                                                                                   <  4:117331/148‑1 (MQ=255)
            tCGTCCCAATGGCGAGAAGAGTTTGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCTTTGCGTAACCCGGGCGACCGCtt                                                                                                                                              <  5:96407/134‑1 (MQ=255)
             cGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTTCATTTCGCGCAGTCATTGCCCGTTCCGTAAACCGGGAGACCGCTTa                                                                                                                                             >  6:96407/1‑134 (MQ=255)
             cGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTa                                                                                                                                <  4:77233/147‑1 (MQ=255)
               tCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCAACCGCTTAGGAGAAGAAGGTACCg                                                                                                                             >  1:145003/1‑148 (MQ=255)
               tCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCAACCGCTTAGGAGAAGAAGGTACCg                                                                                                                             >  1:149452/1‑148 (MQ=255)
                cccAATGGCGAGAAGAGTTGGATCGGCGGATTCATAGCTTAGTGGCGCTGCGTGCCGAACTGGACGGATGTATGGGTTGTGGCTGCCTTTCGCGAAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGTAGAAGAGGGTCCCgg                                                                                                                            >  2:260618/1‑148 (MQ=255)
                cccAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCgg                                                                                                                            >  1:285905/1‑148 (MQ=255)
                cccAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCgg                                                                                                                            >  5:57756/1‑148 (MQ=255)
                  cAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGgcg                                                                                                                          >  7:288664/1‑148 (MQ=255)
                    aTGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCa                                                                                                                        >  8:137788/1‑148 (MQ=255)
                     tGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCCTGCCGAACTGGACGGATGTAATGGTTGTGCCTGCCTATCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCAc                                                                                                                       >  2:129063/1‑148 (MQ=255)
                      ggCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACg                                                                                                                      >  4:166196/1‑148 (MQ=255)
                        cGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCt                                                                                                                    >  7:156521/1‑148 (MQ=255)
                          agaagaGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTg                                                                                                                  <  3:166196/148‑1 (MQ=255)
                              gagTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTgg                                                                                                              <  6:83124/148‑1 (MQ=255)
                                  tGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCAGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTgg                                                                                                              >  3:61458/1‑144 (MQ=255)
                                    gATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGAt                                                                                                         <  8:9430/147‑1 (MQ=255)
                                    gATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGAt                                                                                                         >  7:185342/1‑147 (MQ=255)
                                    gATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATg                                                                                                        <  7:222293/148‑1 (MQ=255)
                                     aTCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGa                                                                                                       <  2:261111/148‑1 (MQ=255)
                                        ggCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTACGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc                                                                                                     <  5:296627/147‑1 (MQ=255)
                                          cgcATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACaa                                                                                                   >  7:54813/1‑147 (MQ=255)
                                          cgcATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCAACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACaaa                                                                                                  <  2:145003/148‑1 (MQ=255)
                                          cgcATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCAACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACaaa                                                                                                  <  2:149452/148‑1 (MQ=255)
                                          cgcATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTCTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCATAGGAGAAGAAGCTACCGGCGCACGCTTGCTGGAAGATGAACaa                                                                                                   >  7:60183/1‑147 (MQ=255)
                                            cATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAAc                                                                                                >  6:230644/1‑148 (MQ=255)
                                               tCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTa                                                                                              >  7:191386/1‑147 (MQ=255)
                                                 aTACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCAGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTaaa                                                                                            <  4:61458/147‑1 (MQ=255)
                                                  tACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAgc                                                                                          >  6:23710/1‑148 (MQ=255)
                                                      ttAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCca                                                                                      >  8:299583/1‑148 (MQ=255)
                                                                  cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                          <  7:137788/148‑1 (MQ=255)
                                                                  cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                          <  8:288664/148‑1 (MQ=255)
                                                                  cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                                                            <  8:185342/146‑1 (MQ=255)
                                                                       cGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgttt                                                                     <  2:285905/148‑1 (MQ=255)
                                                                         aaCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgt                                                                   <  6:57756/148‑1 (MQ=255)
                                                                         aaCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgt                                                                   >  5:259260/1‑148 (MQ=255)
                                                                           cggcaCGGAGGTATTCGTGGGGGCTGCCGTGCGCGCAGTGATTGCCCGTTGCGAACCCCGGGCGACCGCTTCGGAGAAGAAGCCACCGGCGCCCGCTTGCTCGCAGAGCAACAAAACTAACCCGCCCCATGGGCGCTTTAgtttgttt                                                                 <  4:128779/144‑1 (MQ=255)
                                                                                  gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtc                                                          <  5:230644/148‑1 (MQ=255)
                                                                                  gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTAACGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtc                                                          >  8:184463/1‑148 (MQ=255)
                                                                                        ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtc                                                    >  8:183532/1‑148 (MQ=255)
                                                                                            ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCt                                                                                                                                               >  5:252569/1‑53 (MQ=255)
                                                                                            ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                   >  8:316315/1‑145 (MQ=255)
                                                                                            ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                   <  7:316315/145‑1 (MQ=255)
                                                                                            ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtc                                                >  4:44854/1‑148 (MQ=255)
                                                                                                    gCCTTTCGCGCAGTGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATc                                        <  8:156521/148‑1 (MQ=255)
                                                                                                     ccTTTCGCCCAGTGATTGCCCGATGCGTAAACCGGGCGACCGCTTAGGACAAGAATGTAACGGCTAACGCTTCCTGTAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTACGCACTTTCTCATTCTCTCTATcc                                       >  4:253987/1‑148 (MQ=255)
                                                                                                        ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgt                                                           >  5:86516/1‑125 (MQ=255)
                                                                                                        ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgt                                                           <  7:66329/125‑1 (MQ=255)
                                                                                                        ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgt                                                           >  8:66329/1‑125 (MQ=255)
                                                                                                        ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgt                                                           <  6:86516/125‑1 (MQ=255)
                                                                                                        ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACACTGGCTGGAAGATGAACAAAACTAAAGCGCCCCAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                    >  6:43153/1‑148 (MQ=255)
                                                                                                         tCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                   >  7:223923/1‑148 (MQ=255)
                                                                                                          cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                   <  7:183532/147‑1 (MQ=255)
                                                                                                            cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGt                                <  5:23710/148‑1 (MQ=255)
                                                                                                              cAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                                                            <  4:293775/102‑1 (MQ=255)
                                                                                                              cAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                                                            >  3:293775/1‑102 (MQ=255)
                                                                                                                 tGATTGCCGGTTGCGTAACCCGGTCGACCGCTTAGGAGAAGAAGGTACCGGCGCACTCTTGCTGGAAGATGTAGAAAACTAAAGCGGCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCGGCTGGTAcaca                           <  8:60183/148‑1 (MQ=255)
                                                                                                                 tGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAcaca                           <  8:54813/148‑1 (MQ=255)
                                                                                                                  gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACag                          >  2:309553/1‑148 (MQ=255)
                                                                                                                        ccGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGt                    <  7:299583/148‑1 (MQ=255)
                                                                                                                         cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGtt                   >  1:109260/1‑148 (MQ=255)
                                                                                                                            tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTacac                            <  6:148494/136‑1 (MQ=255)
                                                                                                                            tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTacac                            >  5:148494/1‑136 (MQ=255)
                                                                                                                            tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTGTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcc                >  5:13346/1‑148 (MQ=255)
                                                                                                                            tGCGTAACCCGGGCGACCGCATAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCAAAAGGGCGCTTTAGATTGTTTTCCGGTATTTGTCTTTCTCTCTATCCAGCTGGTACACAGGAGGCTTTcc                >  4:4475/1‑148 (MQ=255)
                                                                                                                                tAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATTAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCg            <  1:129063/148‑1 (MQ=255)
                                                                                                                                tAACCCGGGCGACCGCTTAGGAGAAGAACGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCg            <  3:4475/148‑1 (MQ=255)
                                                                                                                                 aaCCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc               >  3:257130/1‑144 (MQ=255)
                                                                                                                                 aaCCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc               <  4:257130/144‑1 (MQ=255)
                                                                                                                                  aCCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGAc          >  4:280164/1‑148 (MQ=255)
                                                                                                                                    ccGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGt        >  1:131087/1‑148 (MQ=255)
                                                                                                                                          gACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacac  <  2:109260/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACAC  >  NC_000913/4277739‑4278024

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