Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,386,750 C→G 5.5% L11V (CTC→GTC)  tyrR → aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,386,7500CG 5.5% 512.3 / 13.5 144L11V (CTC→GTC) tyrRaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase
Reads supporting (aligned to +/- strand):  ref base C (66/70);  new base G (4/4);  total (70/74)
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.93e-01

Only 100 of 148 total aligned reads displayed.

GCCTCGGGGATTTCCGTCGTCAGCTTATCGGTCAGGTGAAAGAAACGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGAT  >  NC_000913/1386493‑1387010
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                    
gcctcgGGGATTTCCGTCGTCAGCTTATCGGTCAGGTGAAAGAAACGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCTTCTTTNTGTCAATGATTATTGACAGAACACTTCCTGTTTGCCAAATTGTGTAATATCATCATACAAATTGTTCTTTTTTTNGGTAAAGGTTCCCATTCTTCTGGAAGGTTTNNTTGAATCCCAACCCGGTCttaaacgcgcc                                                                                                                                                                                                                                                     >  3:77770/1‑265 (MQ=255)
                        ttATCGGTCAGGTGAAAGAAAAGCTGAAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAag                                                                                                                                                                                                                              <  1:57567/274‑1 (MQ=255)
                          aTCGGTCAGGTGAAAGAAACGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGCgg                                                                                                                                                                                                                           >  2:341873/1‑275 (MQ=255)
                                      aaagaaaCGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTAcg                                                                                                                                                                                                               <  1:77770/275‑1 (MQ=255)
                                      aaagaaaCGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGGGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGAGTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTAcg                                                                                                                                                                                                               <  2:167476/275‑1 (MQ=255)
                                              gCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGGCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTTACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTg                                                                                                                                                                                                       <  4:510151/275‑1 (MQ=255)
                                                  cAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTTAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGAtt                                                                                                                                                                                                                   >  4:448208/1‑259 (MQ=255)
                                                      aaGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTAcg                                                                                                                                                                                                               <  2:448209/259‑1 (MQ=255)
                                                                                               tGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAACTAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGACGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTCCTCGCTCTACTCGTGCTACGCGGCCTTGATTTCCGCGGTATTGAGATTGATCCCATTGTGCGAATCTCCCGCAATTTTGCTTACCTGGAGtt                                                                                                                                                      >  4:135583/1‑275 (MQ=255)
                                                                                                     gATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGa                                                                                                                                                                                                <  1:58960/227‑1 (MQ=255)
                                                                                                                  cGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAATCCTTCCTGATATCCAAAAAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCgg                                                                                                                                                                                                                                                                  <  8:681429/148‑1 (MQ=255)
                                                                                                                                      cTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCg                                                                                                                                                                                                                                              <  8:390564/148‑1 (MQ=255)
                                                                                                                                       tttGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGa                                                                                                                                                                                                                                             <  7:242899/148‑1 (MQ=255)
                                                                                                                                            gtCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTAc                                                                                                                                                                                                                                        >  7:753270/1‑148 (MQ=255)
                                                                                                                                            gtCAATGATTGTTGACAGAAACCTTCCTGCGATCCAAATAGTGTCACATCATCATATTAATTGTGCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCggg                                                                                                         <  2:135584/275‑1 (MQ=255)
                                                                                                                                              cAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTa                                                                                                                                                                                                                                                  >  6:643593/1‑136 (MQ=255)
                                                                                                                                                aTGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGAGTCGGTCt                                                                                                                                                                                                                                                               >  5:375840/1‑121 (MQ=255)
                                                                                                                                                 tGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCgg                                                                                                                                                                                                                                                                  >  5:311375/1‑117 (MQ=255)
                                                                                                                                                 tGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTg                                                                                                                                                                                                                                   >  8:964414/1‑148 (MQ=255)
                                                                                                                                                    ttgttgACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTa                                                                                                                                                                                                                                <  8:753270/148‑1 (MQ=255)
                                                                                                                                                    ttgttgACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTAcc                                                                                                  <  2:7552/274‑1 (MQ=255)
                                                                                                                                                       ttgACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAaga                                                                                                                                                                                                                             >  8:974281/1‑148 (MQ=255)
                                                                                                                                                             gAAACCTTCCTGATAGCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTAAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTa                                                                                        <  2:213041/275‑1 (MQ=255)
                                                                                                                                                              aaaCCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTg                                                                                                                                                                                                                      <  7:974281/148‑1 (MQ=255)
                                                                                                                                                                aCCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGAt                                                                                                                                                                                                                    <  7:1070540/148‑1 (MQ=255)
                                                                                                                                                                   ttCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTa                                                                                                                                                                                                                                                  >  5:361116/1‑115 (MQ=255)
                                                                                                                                                                    tCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGAtt                                                                                                                                                                                                                   <  7:471323/145‑1 (MQ=255)
                                                                                                                                                                    tCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGAtt                                                                                                                                                                                                                   >  8:471323/1‑145 (MQ=255)
                                                                                                                                                                    tCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTccc                                                                                 >  2:360162/1‑275 (MQ=255)
                                                                                                                                                                        gCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCcg                                                                                                                                                                                                                                                         <  7:957781/103‑1 (MQ=255)
                                                                                                                                                                           aTCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTAcg                                                                                                                                                                                                               >  7:205615/1‑142 (MQ=255)
                                                                                                                                                                             ccAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACt                                                                                                                                                                                                                                                >  5:998304/1‑107 (MQ=255)
                                                                                                                                                                             ccAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTg                                                                                                                                                                                                       >  5:985457/1‑148 (MQ=255)
                                                                                                                                                                                aaTAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGt                                                                                                                                                                                                                                                                 <  7:370205/87‑1 (MQ=255)
                                                                                                                                                                                aaTAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAag                                                                                                                                                                                                                              >  5:228060/1‑122 (MQ=255)
                                                                                                                                                                                aaTAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTgaga                                                                                                                                                                                                    >  5:645070/1‑148 (MQ=255)
                                                                                                                                                                                               tcatATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGTCCGAAATACGCCGTCTTGCCGGTTTTCCCGATGTGCGTACTGTCCCTTGGTTGCTTTCCGAACGTGAGCATCt                                                      >  2:59612/1‑275 (MQ=255)
                                                                                                                                                                                                catATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGaa                                                                                                                                                                                                                                                     >  8:290991/1‑83 (MQ=255)
                                                                                                                                                                                                catATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGaa                                                                                                                                                                                                                                                     <  7:290991/83‑1 (MQ=255)
                                                                                                                                                                                                catATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGa                                                                                                                                                                                                                                                      <  7:1101907/82‑1 (MQ=255)
                                                                                                                                                                                                 atatTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGa                                                                                                                                                                                   <  6:645070/148‑1 (MQ=255)
                                                                                                                                                                                                  tatTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGa                                                                                                                                                                                                >  7:576120/1‑134 (MQ=255)
                                                                                                                                                                                                  tatTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGa                                                                                                                                                                                                <  8:576120/134‑1 (MQ=255)
                                                                                                                                                                                                     tAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGGCTGGAAGTCTTTTTGGACAACCGACTCGGGCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGtt                                                <  4:17226/275‑1 (MQ=255)
                                                                                                                                                                                                      aaTTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTa                                                                                                                                                                              <  6:985457/148‑1 (MQ=255)
                                                                                                                                                                                                        ttGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACt                                                                                                                                                                                                                                                                     <  5:57361/59‑1 (MQ=255)
                                                                                                                                                                                                        ttGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTg                                                                                                                                                                                                       <  7:708015/121‑1 (MQ=255)
                                                                                                                                                                                                         tGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTgaga                                                                                                                                                                                                    >  5:854003/1‑123 (MQ=255)
                                                                                                                                                                                                            tCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCa                                                                                                                                                                                           <  5:1079533/129‑1 (MQ=255)
                                                                                                                                                                                                              ttttttCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTg                                                                                                                                                                                                                                   <  7:982634/87‑1 (MQ=255)
                                                                                                                                                                                                               tttttCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTa                                                                                                                                                                              >  6:372816/1‑139 (MQ=255)
                                                                                                                                                                                                               tttttCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTa                                                                                                                                                                              <  5:372816/139‑1 (MQ=255)
                                                                                                                                                                                                                      ggTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGa                                                                                                                                                                                   >  7:669626/1‑127 (MQ=255)
                                                                                                                                                                                                                      ggTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGa                                                                                                                                                                                   <  8:669626/127‑1 (MQ=255)
                                                                                                                                                                                                                        tGAAGGTTCCCATGGGTCTTGAATGCTTTTTTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGTTCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCg                             <  4:341872/275‑1 (MQ=255)
                                                                                                                                                                                                                        tGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCg                             <  2:211803/275‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGGTTCCCATGCGTCTGGAAGTCTGTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCa                                                                                                                                                                                           <  8:728033/116‑1 (MQ=255)
                                                                                                                                                                                                                              ttCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATAGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGtt                                                                                                                                                      >  5:325203/1‑148 (MQ=255)
                                                                                                                                                                                                                                  cATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCa                                                                                                                                                                                                                         >  5:466648/1‑77 (MQ=255)
                                                                                                                                                                                                                                  cATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTgag                                                                                                                                                  <  8:26003/148‑1 (MQ=255)
                                                                                                                                                                                                                                  cATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTgag                                                                                                                                                  <  6:325203/148‑1 (MQ=255)
                                                                                                                                                                                                                                  cATGCGTCTGAAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCa                                                                                                                                                                                                                         <  6:466648/77‑1 (MQ=255)
                                                                                                                                                                                                                                      cGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGtt                                                                                                                                              <  8:223566/148‑1 (MQ=255)
                                                                                                                                                                                                                                        tCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGa                                                                                                                                                                                                <  5:28894/96‑1 (MQ=255)
                                                                                                                                                                                                                                        tCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGa                                                                                                                                                                                                >  6:28894/1‑96 (MQ=255)
                                                                                                                                                                                                                                         cTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCa                                                                                                                                                                         <  6:479342/118‑1 (MQ=255)
                                                                                                                                                                                                                                         cTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTca                                                                                                                                           >  5:759958/1‑148 (MQ=255)
                                                                                                                                                                                                                                          tGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGACTTACTCGATCTACTCGTGCTACGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGCAATACGCCTTATTGCGGGTTTTAACGATGTTCGTACTGTCCCGTGGATGCCTTCCGAACGTGACCTTCTGGCGTTTGGCGCGTTACTGGGGGCGTTGCCTTACCCCGTGccc           >  4:425407/1‑273 (MQ=255)
                                                                                                                                                                                                                                          tGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCt                                                                                                                                                                 >  4:281375/1‑125 (MQ=255)
                                                                                                                                                                                                                                          tGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCCCGTTACTGGAGGCGTTGCCTGAACCtggtctc           >  2:18511/1‑270 (MQ=255)
                                                                                                                                                                                                                                            gAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTgg                                                                                                                                                          <  6:354462/130‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTg                                                                                                                                                                   >  7:1079148/1‑120 (MQ=255)
                                                                                                                                                                                                                                             aaGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTg                                                                                                                                                                   <  8:1079148/120‑1 (MQ=255)
                                                                                                                                                                                                                                              aGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCAt                                                                                                                                                                                          <  8:108654/96‑1 (MQ=255)
                                                                                                                                                                                                                                              aGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTAcc                                                                                                                                                                            >  6:296858/1‑110 (MQ=255)
                                                                                                                                                                                                                                              aGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAAc                                                                                                                                                             <  2:281376/125‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCg    >  2:1467/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  ttttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGa                                                                                                                                  <  6:759958/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGa                                                                                                                                  <  6:374309/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGa   >  2:345085/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   tttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATcc                                                                                                                                                                                             >  5:654623/1‑88 (MQ=255)
                                                                                                                                                                                                                                                   tttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCaa                                                                                                                                                                        <  8:173776/109‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCACATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGAt                                                                                                                                 >  7:788757/1‑148 (MQ=255)
                                                                                                                                                                                                                                                       tgAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAAtttt                                                                                                                                                                    <  6:672561/109‑1 (MQ=255)
                                                                                                                                                                                                                                                        gAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTc                                                                                                                                            >  8:242079/1‑132 (MQ=255)
                                                                                                                                                                                                                                                        gAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTc                                                                                                                                            <  7:242079/132‑1 (MQ=255)
                                                                                                                                                                                                                                                         caGACCGAGTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGa                                                                                                                                  >  7:881123/2‑141 (MQ=255)
                                                                                                                                                                                                                                                         caGACCGAGTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGa                                                                                                                                  <  8:881123/140‑1 (MQ=255)
                                                                                                                                                                                                                                                         caGAACGAGTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTg                                                                                                                                                                                        >  6:129875/2‑87 (MQ=255)
                                                                                                                                                                                                                                                         aaGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGaa                                                                                                                                                                                  <  6:542521/93‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGaa                                                                                                                                                                                  >  5:542521/1‑93 (MQ=255)
                                                                                                                                                                                                                                                          aGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTggg                                                                                                                                                                                      >  7:935197/1‑88 (MQ=255)
                                                                                                                                                                                                                                                          aGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTggg                                                                                                                                                                                      <  8:935197/88‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccGAGTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACt                                                                                                                                                            <  6:131577/111‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccGAGTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACt                                                                                                                                                            >  5:131577/1‑111 (MQ=255)
                                                                                                                                                                                                                                                               gACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACg                                                                                                                     >  8:393266/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               gACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACg                                                                                                                     >  7:212915/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 cTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGt                                                                                                                                                       >  6:920720/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                 cTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGttt                                                                                                                                             <  8:150138/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGttt                                                                                                                                             >  7:150138/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                    
GCCTCGGGGATTTCCGTCGTCAGCTTATCGGTCAGGTGAAAGAAACGCTGCAAAAAGGCAAAACGCCCAGCGAAAAATAATGCAATATCGGGTGCTGACCGGATATCTTTACGCCGAAGTGCCCGTTTTTCCGTCTTTGTGTCAATGATTGTTGACAGAAACCTTCCTGCTATCCAAATAGTGTCATATCATCATATTAATTGTTCTTTTTTCAGGTGAAGGTTCCCATGCGTCTGGAAGTCTTTTGTGAAGACCGACTCGGTCTGACCCGCGAATTACTCGATCTACTCGTGCTAAGAGGCATTGATTTACGCGGTATTGAGATTGATCCCATTGGGCGAATCTACCTCAATTTTGCTGAACTGGAGTTTGAGAGTTTCAGCAGTCTGATGGCCGAAATACGCCGTATTGCGGGTGTTACCGATGTGCGTACTGTCCCGTGGATGCCTTCCGAACGTGAGCATCTGGCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGAT  >  NC_000913/1386493‑1387010

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