Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,207,557 C→T 5.8% N19N (AAC→AAT)  stfP → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,5570CT 5.8% 223.4 / 2.3 69N19N (AAC→AAT) stfPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (28/37);  new base T (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CGTATCTGGGGTGACACCGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTAC  >  NC_000913/1207411‑1207703
                                                                                                                                                  |                                                                                                                                                  
cGTATCTGGGGTGACACCGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACgg                                                                                                                                                  <  3:242277/149‑1 (MQ=255)
   aTCTGGGGTGACACCGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATTGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGttt                                                                                                                                               >  3:186099/1‑149 (MQ=255)
        gggTGACACCGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTAccc                                                                                                                                          <  8:294187/149‑1 (MQ=255)
                ccGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAAc                                                                                                                                  <  7:323736/149‑1 (MQ=255)
                             tgtgCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCa                                                                                                                     <  1:23692/149‑1 (MQ=255)
                                gCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAATGGTTTTACCCGTGGTAACCCCCAGACCGGCACGc                                                                                                                  >  1:161054/1‑149 (MQ=255)
                                gCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAATGGTTTTACCCGTGGTAACCCCCAGACCGGCACGc                                                                                                                  >  1:40946/1‑149 (MQ=255)
                                    aaCAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCt                                                                                                                <  6:270683/147‑1 (MQ=255)
                                     aCAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGc                                                                                                              <  7:81014/148‑1 (MQ=255)
                                     aCAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGc                                                                                                              <  4:96777/148‑1 (MQ=255)
                                          ctctGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGa                                                                                                        <  7:395151/149‑1 (MQ=255)
                                          ctctGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGa                                                                                                        >  8:256561/1‑149 (MQ=255)
                                               cccGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAATGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTg                                                                                                    <  2:40946/148‑1 (MQ=255)
                                               cccGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAATGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTg                                                                                                    <  2:161054/148‑1 (MQ=255)
                                               cccGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTgg                                                                                                   <  5:423678/149‑1 (MQ=255)
                                               cccGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTg                                                                                                    >  2:108248/1‑148 (MQ=255)
                                               cccGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTg                                                                                                    >  7:436833/1‑148 (MQ=255)
                                                 cGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgat                                                                                                 <  1:105388/149‑1 (MQ=255)
                                                 cGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgat                                                                                                 <  3:330705/149‑1 (MQ=255)
                                                   tCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatga                                                                                               <  3:197926/149‑1 (MQ=255)
                                                   tCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatga                                                                                               <  1:206073/149‑1 (MQ=255)
                                                    cGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatgat                                                                                              >  2:29416/1‑149 (MQ=255)
                                                     ccATACCTACTTAATTTTTAAATATCAGCAGTACTACATGCATCGTATACACACGAAAACCCCGAAGAACGATAAGTTCGGCGCTGGTAAGAACTGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatga                                                                                               <  4:265798/146‑1 (MQ=25)
                                                      cATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGAtgatga                                                                                            <  5:219564/149‑1 (MQ=255)
                                                       aTACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGac                                                                                           >  6:233647/1‑149 (MQ=255)
                                                           cTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGactact                                                                                       >  3:185907/1‑149 (MQ=255)
                                                              cGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACttt                                                                                     <  6:8314/148‑1 (MQ=255)
                                                              cGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACtt                                                                                      <  1:81199/147‑1 (MQ=255)
                                                               gTAATTTTTAAATATCCGGAGT‑ATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGa                                                                                   >  2:3855/1‑148 (MQ=255)
                                                                tAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTg                                                                                                                                       <  1:152923/96‑1 (MQ=255)
                                                                tAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTg                                                                                                                                       >  2:152923/1‑96 (MQ=255)
                                                                tAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGAc                                                                                  <  5:233647/149‑1 (MQ=255)
                                                                   tttttAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATg                                                                               >  1:35761/1‑149 (MQ=255)
                                                                     tttAAATATCCGGAGT‑ATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGtt                                                                             <  1:3855/148‑1 (MQ=255)
                                                                          atatCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCagg                                                                        <  1:108248/149‑1 (MQ=255)
                                                                           tatCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCagga                                                                       <  1:244621/149‑1 (MQ=255)
                                                                                   gTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACttt                                                               <  2:35761/149‑1 (MQ=255)
                                                                                    tAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGAt                                                                                                       <  8:25572/108‑1 (MQ=255)
                                                                                    tAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGAt                                                                                                       >  7:25572/1‑108 (MQ=255)
                                                                                    tAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTg                                                              <  3:30380/149‑1 (MQ=255)
                                                                                      aTCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCa                                                            <  4:191873/149‑1 (MQ=255)
                                                                                      aTCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTCGACATGTTGCAGGAGGAACTTTGCa                                                            <  1:81264/149‑1 (MQ=255)
                                                                                        ccATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGc                                                          >  2:426505/1‑149 (MQ=255)
                                                                                          aTGCATAGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCTGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCgt                                                        >  7:179719/1‑149 (MQ=255)
                                                                                           tGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGATCTTTGCAGCgtg                                                       <  1:426505/149‑1 (MQ=255)
                                                                                                 gTATAGACACGAAACCCGCGCAGCAGGATAAGTGCGGCGCGGGTAAGAACGGCTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGa                                                  <  6:272982/148‑1 (MQ=255)
                                                                                                    taGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCa                                              >  1:116877/1‑149 (MQ=255)
                                                                                                         acGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCAcc                                                                                                          >  2:263901/1‑84 (MQ=255)
                                                                                                         acGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCAcc                                                                                                          <  1:263901/84‑1 (MQ=255)
                                                                                                         acGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCgg                                         >  8:50463/1‑149 (MQ=255)
                                                                                                         acGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCgg                                         >  5:311798/1‑149 (MQ=255)
                                                                                                         acGAAAACCGCGCAGAAGGAGAATTTCGGCGCGGGTAAGAACGGTTTTCCCCGTGGTTACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGAGATGTTGCAGGAGGACCGATGCAGCGAGGTCGCGGCATCCgc                                         >  7:312525/1‑148 (MQ=21)
                                                                                                           gAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTg                                       >  7:386418/1‑149 (MQ=255)
                                                                                                              aaCCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCTGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCa                                    >  6:248906/1‑149 (MQ=255)
                                                                                                                   cgcAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTg                               <  4:185907/149‑1 (MQ=255)
                                                                                                                        aaGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGgag                            >  5:428844/1‑147 (MQ=255)
                                                                                                                                gTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGCCCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAg                                                 >  8:416981/1‑118 (MQ=255)
                                                                                                                                gTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatgat                                                                                              >  6:95008/1‑73 (MQ=255)
                                                                                                                                gTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGgatgat                                                                                              <  5:95008/73‑1 (MQ=255)
                                                                                                                                     gcgcGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGAc             <  6:428844/149‑1 (MQ=255)
                                                                                                                                         gggTAAGACCGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGt         >  5:366927/1‑149 (MQ=255)
                                                                                                                                         gggTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGt         <  1:29416/149‑1 (MQ=255)
                                                                                                                                           gTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGtta       >  7:271280/1‑149 (MQ=255)
                                                                                                                                               gAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGgaga                           <  4:365471/125‑1 (MQ=255)
                                                                                                                                               gAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGgaga                           >  3:365471/1‑125 (MQ=255)
                                                                                                                                               gAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGttactta   <  7:256561/149‑1 (MQ=255)
                                                                                                                                                aaCGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAg                                                 <  4:144621/102‑1 (MQ=255)
                                                                                                                                                aaCGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAg                                                 >  3:144621/1‑102 (MQ=255)
                                                                                                                                                aaCGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGttacttac  >  7:331993/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTATCTGGGGTGACACCGTTGTCGAGTGTGTGCTTAACAAACTCTGCCCGTCGCATACCTACGTAATTTTTAAATATCCGGAGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTAC  >  NC_000913/1207411‑1207703

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