Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,766,262 A→T 5.6% F808L (TTT→TTA)  ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,766,2620AT 5.6% 254.2 / 4.8 71F808L (TTT→TTA) ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (2/2);  total (36/35)
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.16e-01

CAGCTTCGCCATACGGTTGAGGAAATCCGCCGTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATTTGTTGTTGTAGC  >  NC_000913/1766115‑1766409
                                                                                                                                                   |                                                                                                                                                   
cAGCTTCGCCATACGGTTGAGGAAATCCGCCGTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACa                                                                                                                                                     >  4:1006310/1‑148 (MQ=255)
       gCCATACGGTTGAGGAAATCCGCCGTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCg                                                                                                                                              <  3:1006310/148‑1 (MQ=255)
         cATACGGTTGAGGAAATCCGCCGTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGcc                                                                                                                                            <  1:1109027/148‑1 (MQ=255)
                             ccgTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTa                                                                                                                           >  2:300575/1‑145 (MQ=255)
                             ccgTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTa                                                                                                                           <  1:300575/145‑1 (MQ=255)
                                tCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTGATATGCTGGGCGTTGCCATTAGACGAAAATGGCAGTAACACAGGCTGGAAACCTAACTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCt                                                                                                                     <  1:1110808/148‑1 (MQ=255)
                                tCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTg                                                                                                                                    >  2:6598/1‑133 (MQ=255)
                                tCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTg                                                                                                                                    <  1:6598/133‑1 (MQ=255)
                                   ttttCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACAGGCAAAATCAGCCACCACTTGCGCATCGTAATAGCTgga                                                                                                                  >  2:1110808/1‑147 (MQ=255)
                                                    aaaCGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACaaa                                                                                                                                                   >  2:989992/1‑98 (MQ=255)
                                                    aaaCGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACaaa                                                                                                                                                   <  1:989992/98‑1 (MQ=255)
                                                               aaaaCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAAt                                                                                                                         >  1:1061711/1‑113 (MQ=255)
                                                               aaaaCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAAt                                                                                                                         <  2:1061711/113‑1 (MQ=255)
                                                               aaaaCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACtgtg                                                                                      >  4:484575/1‑148 (MQ=255)
                                                                  aCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACTAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGt                                                                                                          >  1:1039337/1‑125 (MQ=255)
                                                                  aCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACTAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGt                                                                                                          <  2:1039337/125‑1 (MQ=255)
                                                                  aCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACTAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGt                                                                                                          >  3:828731/1‑125 (MQ=255)
                                                                  aCCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACTAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGt                                                                                                          <  4:828731/125‑1 (MQ=255)
                                                                            tGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTg                                                                                                                    >  1:535649/1‑105 (MQ=255)
                                                                            tGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTg                                                                                                                    <  2:535649/105‑1 (MQ=255)
                                                                             gCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCg                                                                                                                                              <  3:671247/78‑1 (MQ=255)
                                                                             gCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCg                                                                                                                                              >  4:671247/1‑78 (MQ=255)
                                                                                   cTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAgg                                                                                                             >  1:669594/1‑105 (MQ=255)
                                                                                   cTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAgg                                                                                                             <  2:669594/105‑1 (MQ=255)
                                                                                      tGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAac                                                                                            <  3:462040/119‑1 (MQ=255)
                                                                                      tGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAac                                                                                            >  4:462040/1‑119 (MQ=255)
                                                                                           tttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTg                                                                 <  1:519482/141‑1 (MQ=255)
                                                                                           tttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTg                                                                 >  2:519482/1‑141 (MQ=255)
                                                                                            ttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAg                                                                                                              <  3:857291/95‑1 (MQ=255)
                                                                                            ttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAg                                                                                                              >  4:857291/1‑95 (MQ=255)
                                                                                            ttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTc                                                                                                         <  4:515233/100‑1 (MQ=255)
                                                                                            ttGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTc                                                                                                         >  3:515233/1‑100 (MQ=255)
                                                                                                  aaaaTGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTg                                                   >  1:382665/1‑148 (MQ=255)
                                                                                                    aaTGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGtt                                                 <  2:382665/148‑1 (MQ=255)
                                                                                                              acacAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCa                                                        >  3:566125/1‑131 (MQ=255)
                                                                                                              acacAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCa                                                        <  4:566125/131‑1 (MQ=255)
                                                                                                                    gCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAgg                                                                                                             >  3:247999/1‑72 (MQ=255)
                                                                                                                    gCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAgg                                                                                                             <  4:247999/72‑1 (MQ=255)
                                                                                                                       ggAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgc                                                                          >  4:685706/1‑104 (MQ=255)
                                                                                                                       ggAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgc                                                                          <  3:685706/104‑1 (MQ=255)
                                                                                                                        gAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAg                             >  3:450678/1‑148 (MQ=255)
                                                                                                                             cTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATg                        <  3:484575/148‑1 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgc                                                                          >  2:989719/1‑97 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgc                                                                          <  1:989719/97‑1 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCAt                                                                   <  4:1061659/104‑1 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCAt                                                                   >  3:1061659/1‑104 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATg                                      >  2:109854/1‑133 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATg                                      <  1:109854/133‑1 (MQ=255)
                                                                                                                              tAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGc                       >  2:544697/1‑148 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAac                                                                                            <  1:746720/78‑1 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAac                                                                                            >  2:746720/1‑78 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCt                                                                                >  1:546411/1‑90 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCt                                                                                <  4:302703/90‑1 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCt                                                                                <  2:546411/90‑1 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCt                                                                                >  3:302703/1‑90 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgctctgc                                                                     <  1:223944/101‑1 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTtctgctctgc                                                                     >  2:223944/1‑101 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCAt                                       <  3:737210/131‑1 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCAt                                       >  4:737210/1‑131 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTg                                  >  1:166298/1‑136 (MQ=255)
                                                                                                                               aaTTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTg                                  <  2:166298/136‑1 (MQ=255)
                                                                                                                                        aCCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCAttt             >  1:523577/1‑148 (MQ=255)
                                                                                                                                        aCCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCAttt             >  4:476317/1‑148 (MQ=255)
                                                                                                                                          cAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATttgt           <  4:450678/148‑1 (MQ=255)
                                                                                                                                            gACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTc                                                         >  4:701875/1‑100 (MQ=255)
                                                                                                                                            gACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTc                                                         <  3:701875/100‑1 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTga                                                                >  2:102167/1‑92 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTga                                                                <  1:102167/92‑1 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATttgt           <  3:720733/145‑1 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATttgt           >  4:720733/1‑145 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCATGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATttgt           <  1:1111050/145‑1 (MQ=255)
                                                                                                                                             aCGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCATGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATttgt           >  2:1111050/1‑145 (MQ=255)
                                                                                                                                                gACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTc                                                >  4:692302/1‑105 (MQ=255)
                                                                                                                                                gACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTc                                                <  3:692302/105‑1 (MQ=255)
                                                                                                                                                   aaaaTCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATTTGTTGTTGTAGc  >  1:163186/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCTTCGCCATACGGTTGAGGAAATCCGCCGTCTTTTTCGCCGTCTTCGCAAAACGATTAAGAAAACCTTTAATATGCTGGGCTTTGCCATTTGGCGAAAATGGCAGTAACACAGGCTGGAAACCTAATTTTTCGACCAGACGGACAAAATCCGCCACCACTTGCGCATCGTAATAGCTGGTAAAGGGGTCCTGCACCACCAACACTGTGCGCGCTTTCTGCTCTGCATTGAGAGATTCAAGCTGTTCCAGCGTCATGTTTGCCGAGCGATGCCCCACCATTTGTTGTTGTAGC  >  NC_000913/1766115‑1766409

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