Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,956,702 C→T 5.4% A102A (GCG→GCA torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,956,7020CT 5.4% 242.1 / 2.6 74A102A (GCG→GCAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (38/32);  new base T (2/2);  total (40/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 = 8.30e-01

ACACGCCAGAAGAGCGCCAGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCT  >  NC_000913/1956557‑1956844
                                                                                                                                                 |                                                                                                                                              
acacGCCAGAAGAGCGCCAGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTACCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGc                                                                                                                                              >  6:129446/1‑148 (MQ=255)
         aagaGCGCCAGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAg                                                                                                                                     >  1:352517/1‑148 (MQ=255)
            agCGCCAGCCGTAAGATCCGGCCAAAATGGCCGAGGGGCCGTACGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATCGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTgg                                                                                                                                  <  2:294160/148‑1 (MQ=255)
               gcCAGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGAt                                                                                                                                >  6:42838/1‑147 (MQ=255)
                  aGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGt                                                                                                                                          <  2:161467/134‑1 (MQ=255)
                  aGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGt                                                                                                                                          >  1:161467/1‑134 (MQ=255)
                   gCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCg                                                                                                                           <  7:229400/148‑1 (MQ=255)
                   gCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCAGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCg                                                                                                                           <  6:316283/148‑1 (MQ=255)
                     cGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCgag                                                                                                                         <  4:153244/148‑1 (MQ=255)
                       tAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGAt                                                                                                                       >  3:187307/1‑148 (MQ=255)
                        aaGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCgaga                                                                                                                        >  2:189343/1‑146 (MQ=255)
                             ccGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCttt                                                                                                                 >  1:112034/1‑148 (MQ=255)
                              cGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCtttt                                                                                                                <  4:187307/148‑1 (MQ=255)
                               ggCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCttttt                                                                                                               <  2:352517/148‑1 (MQ=255)
                                  aaaaaTGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTtctc                                                                                                            <  3:107384/148‑1 (MQ=255)
                                    aaaTGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGGAGt                                                                                                                                    <  8:94920/122‑1 (MQ=255)
                                    aaaTGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGGAGt                                                                                                                                    >  7:94920/1‑122 (MQ=255)
                                         gCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTTGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCg                                                                                                     >  7:221451/1‑148 (MQ=255)
                                            gATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTTGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGAt                                                                                                  <  8:221451/148‑1 (MQ=255)
                                             aTGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATg                                                                                                 >  1:199232/1‑148 (MQ=255)
                                              tGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGc                                                                                                <  1:399768/148‑1 (MQ=255)
                                              tGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGc                                                                                                >  4:24949/1‑148 (MQ=255)
                                                 tCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTg                                                                                              >  1:127117/1‑147 (MQ=255)
                                                      tAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATAcg                                                                                        >  1:49197/1‑148 (MQ=255)
                                                       aGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATAcg                                                                                        <  2:199232/147‑1 (MQ=255)
                                                            tACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTTGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCgt                                                                                  >  3:247464/1‑148 (MQ=255)
                                                            tACGAATACGATCATGTTGCTCATGAATCAGGGTTAATGCCTGTTCCCCGCTGACCTGTACATAGGTATCTTCGCCACGACCTTGCGGCGGGTGCAGGGGATTTTCGAGATAGCTTTTTCTCACCATCGGATGCTGCATACGCGCCgt                                                                                  >  4:81764/1‑148 (MQ=255)
                                                                aaTACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCgtggtg                                                                              >  1:211173/1‑148 (MQ=255)
                                                                  tACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTgg                                                                                                                                  <  1:218633/94‑1 (MQ=255)
                                                                  tACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTgg                                                                                                                                  >  2:218633/1‑94 (MQ=255)
                                                                        cATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTTTACATAGGTATCTTAGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGAGCCGTGGTGTGTACCTg                                                                      >  3:56340/1‑148 (MQ=255)
                                                                        cATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCga                                                                                                                          >  3:56096/1‑96 (MQ=255)
                                                                        cATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCga                                                                                                                          <  4:56096/96‑1 (MQ=255)
                                                                          tGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGAt                                                                    <  5:42838/148‑1 (MQ=255)
                                                                            ttGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTTGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATcc                                                                  <  4:247464/148‑1 (MQ=255)
                                                                            ttGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATcc                                                                  >  6:315025/1‑148 (MQ=255)
                                                                            ttGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATc                                                                   <  5:129446/147‑1 (MQ=255)
                                                                                tCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCt                                                                                               >  7:106564/1‑115 (MQ=255)
                                                                                tCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCt                                                                                               <  8:106564/115‑1 (MQ=255)
                                                                                 cATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATccgccgc                                                             >  7:23341/1‑148 (MQ=255)
                                                                                     aaTCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGt                                                          <  8:23341/147‑1 (MQ=255)
                                                                                      aTCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATccgccgccg                                                           >  1:81480/1‑145 (MQ=255)
                                                                                        cAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGAc                                                      >  3:29494/1‑148 (MQ=255)
                                                                                        cAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTATAACCATCGGATGATGAATACGCGCCGTGGTGTGTACCTGATCCTACGCCGTAGAc                                                      >  4:15908/1‑148 (MQ=255)
                                                                                         aGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACt                                                     >  4:120545/1‑148 (MQ=255)
                                                                                           ttttAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGt                                                   <  2:112034/148‑1 (MQ=255)
                                                                                           tttaAATGCCTGTTCCCAGCTCCCCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACCCGCCGTGGTGTGTACCTGAACCGCCGCCGAACACTg                                                    >  6:25443/1‑147 (MQ=255)
                                                                                             ttAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTa                                                  <  2:127117/147‑1 (MQ=255)
                                                                                              tAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCCCGACATTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAg                                                >  1:312964/1‑148 (MQ=255)
                                                                                              tAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATa                                                                                          <  2:78896/106‑1 (MQ=255)
                                                                                              tAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATa                                                                                          >  1:78896/1‑106 (MQ=255)
                                                                                                 tGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGaa                                             <  7:58645/148‑1 (MQ=255)
                                                                                                   ccTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAAtt                                           >  3:249832/1‑148 (MQ=255)
                                                                                                   ccTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAAtt                                           >  7:221377/1‑148 (MQ=255)
                                                                                                    gtGTTCACAGTTCAACTGTACTTAGGTATATTCTCGATGATCTTTCGCGTTTTGCTGTGGTTTATATAGATAGATTTTTCTCAGAATCGGATGCTTAATAATCGAATTGTTGTGTACCTTATCCGCCGCCGTAGACTGTAAGGAAtt                                           <  8:221377/146‑1 (MQ=255)
                                                                                                                       aCATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAgcg                       >  6:191287/1‑148 (MQ=255)
                                                                                                                         aTAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGt                                                          <  2:258538/111‑1 (MQ=255)
                                                                                                                         aTAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCAGCCGCCGt                                                          >  1:258538/1‑111 (MQ=255)
                                                                                                                           aGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCTCGcc                   >  4:196055/1‑148 (MQ=255)
                                                                                                                            ggTATCTTGGAGACGACCTTTCTCCTGGTGCACTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCTCCGTGGTGTGTACCGGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGcc                   <  6:82333/147‑1 (MQ=255)
                                                                                                                             gTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGTTGCTGAATACGCGCCGTGGTGTGTACGTGATCCGGCGCCGTAGACTGTAAGGTATTCGGTATGGTTTTCGCCAGCGCGCCTg                 <  1:189343/148‑1 (MQ=255)
                                                                                                                                tCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGtt               <  2:211173/147‑1 (MQ=255)
                                                                                                                                 cTTCGCCACGACCTTTCGCCGCTTAAAGTGGATTATCGAGATAGCAATTTCTCACCATCCGAGGCTGATTACGCGCAGTGGTGTGTACCTGAACCGCCGCCGTAGACTCTAAGGAATTCGGTATGGTTTTCGCAAGCGCGCCTGATg              >  4:321270/1‑147 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCgg                                        <  7:123151/115‑1 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCgg                                        >  8:123151/1‑115 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCt                  <  4:29176/137‑1 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCt                  >  3:29176/1‑137 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGaaa        >  8:263869/1‑147 (MQ=255)
                                                                                                                                       cACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAAc       <  3:348761/148‑1 (MQ=255)
                                                                                                                                          gACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAAt    >  7:374942/1‑148 (MQ=255)
                                                                                                                                           aCCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATc   <  5:191287/148‑1 (MQ=255)
                                                                                                                                            ccTTTCGCCGGTGGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCt  <  3:161798/148‑1 (MQ=255)
                                                                                                                                                tCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTa                                                         <  6:290971/89‑1 (MQ=255)
                                                                                                                                                tCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTa                                                         >  5:290971/1‑89 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACACGCCAGAAGAGCGCCAGCCGTAAGATCCGGCAAAAATGGCCGATGGTCCGTTAGCTTTACGAATACGATCATGTTGCTCATGAATCAGTTTTAATGCCTGTTCCCAGCTCACCTGTACATAGGTATCTTCGCCACGACCTTTCGCCGGTTGCAGTGGATTATCGAGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCT  >  NC_000913/1956557‑1956844

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