Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,741 T→G 16.2% intergenic (+84/‑27) ycbF → / → pyrD putative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,004,7410TG16.2% 209.6 / 40.0 93intergenic (+84/‑27)ycbF/pyrDputative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base T (36/42);  new base G (7/8);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCA  >  NC_000913/1004642‑1004838
                                                                                                   |                                                                                                 
gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                                                                                                   >  1:838380/1‑100 (MQ=255)
gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                                                                                                   >  2:926768/1‑100 (MQ=255)
 cAAAATCTGCCCTGAACCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                  <  1:1517604/100‑1 (MQ=255)
 cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTTTTTTTGTTTATTTTTCTTTTTTTTTTTTCTCGCACATAATACGCCCCCGGttt                                                                                                  >  2:2537689/1‑100 (MQ=255)
  aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTg                                                                                                 <  1:186676/100‑1 (MQ=255)
  aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCAGGTTTg                                                                                                 >  2:1877157/1‑100 (MQ=255)
   aaaTCTGCCCCGAAACAGGTTCGGAGAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTg                                                                                                 <  1:527455/99‑1 (MQ=255)
     aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTg                                                                                                 <  2:1814043/97‑1 (MQ=255)
          cccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCgg                                                                                         <  1:2439018/100‑1 (MQ=255)
           ccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCgg                                                                                         >  1:1737317/1‑99 (MQ=255)
            cTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGa                                                                                       >  2:1877261/1‑100 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                                                                                      <  1:233210/100‑1 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                      <  2:472376/100‑3 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAAt                                                                                     >  2:9670/1‑100 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAAt                                                                                     >  2:1191304/1‑100 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  1:77895/100‑4 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  1:1929475/100‑4 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  2:2435350/100‑4 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  1:652608/100‑4 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  2:1630030/100‑4 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                     <  1:1351410/100‑4 (MQ=255)
                aaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATcc                                                                                   <  1:965736/100‑1 (MQ=255)
                 aCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCa                                                                                  <  1:2507197/100‑1 (MQ=255)
                  cAGGTTCGGAGAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGCACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATcc                                                                                   <  1:1773644/98‑1 (MQ=255)
                  cAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATc                                                                                    <  2:731030/97‑1 (MQ=255)
                   aGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAgg                                                                                <  2:2544095/100‑1 (MQ=255)
                   aGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggcgtat                                                                                 <  1:2550800/99‑8 (MQ=38)
                    ggTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGga                                                                               >  2:80030/1‑100 (MQ=255)
                     gTTCGGAAAACGTTTGCGTTTTTTTTGTTGTTGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACTCCCCCGGTTTGCACACCGGGAATCCAGga                                                                               >  1:2128847/1‑99 (MQ=255)
                     gTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCTAGgag                                                                              <  1:2270143/100‑1 (MQ=255)
                           aaaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAAATCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa                                                                        >  1:2555262/1‑100 (MQ=255)
                           aaaaCGTTTGCGTTATTTTTACCACAGGCAAATTCCCTTTTGGTAAGAACTCGCCCATAATACGCCCCAGGATTGAAAACCGGGAATCCaaaaaagtt                                                                          >  1:2473911/1‑90 (MQ=255)
                            aaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa                                                                        <  2:1865430/99‑1 (MQ=255)
                              aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGt                                                                     <  1:591650/100‑1 (MQ=255)
                                gTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac                                                                   <  1:1589112/100‑1 (MQ=255)
                                    gCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcc                                                               <  2:619161/100‑1 (MQ=255)
                                     cGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAccc                                                              >  2:2499233/1‑100 (MQ=255)
                                      gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTTCGAACTCGCACATAATACGCCCCCGCTTTGCACACCGGGAATCCAGGAGTGTTCATGTACTAcccc                                                             >  1:1401815/1‑100 (MQ=255)
                                      gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcccc                                                             >  2:233101/1‑100 (MQ=255)
                                       ttttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCt                                                            >  2:883988/1‑100 (MQ=255)
                                          tttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcg                                                         <  1:601073/100‑1 (MQ=255)
                                            tttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgtt                                                       <  2:1679107/100‑1 (MQ=255)
                                             ttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttc                                                      <  2:2135032/100‑1 (MQ=255)
                                              tGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                     <  2:128834/100‑1 (MQ=255)
                                              tGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                     >  1:51660/1‑100 (MQ=255)
                                              tGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACAGAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                     <  2:1045732/100‑1 (MQ=255)
                                               gCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTtcgttcgt                                                    >  1:1167127/1‑100 (MQ=255)
                                                ccGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                   <  1:1319432/100‑1 (MQ=255)
                                                  gCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaa                                                  >  1:254260/1‑99 (MQ=255)
                                                   cAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAg                                                <  1:554980/100‑1 (MQ=255)
                                                    aGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAg                                                <  1:371499/99‑1 (MQ=255)
                                                    aGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGc                                               <  2:1552806/100‑1 (MQ=255)
                                                    aGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGc                                               <  2:2136484/100‑1 (MQ=255)
                                                     ggTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGcc                                              >  1:2491508/1‑100 (MQ=255)
                                                          aTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGTTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCtttt                                         <  1:219102/100‑1 (MQ=255)
                                                            tCCCTTTAGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACATCGGGAATCCAGGAGAGTTCAGGTACTACCCCTTCGTTCGTAAAGCCCTTTTcc                                       <  2:588383/100‑1 (MQ=255)
                                                             cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                      >  2:1251245/1‑100 (MQ=255)
                                                              ccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAg                                     >  1:723609/1‑100 (MQ=255)
                                                                 tttGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTc                                  <  1:326173/100‑1 (MQ=255)
                                                                    ggTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACAGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGAt                               >  1:980505/1‑100 (MQ=255)
                                                                     gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATc                              >  2:416407/1‑100 (MQ=255)
                                                                     gTCCGAACTCGCACATAATACACCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATc                              >  2:81334/1‑100 (MQ=255)
                                                                      tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                             >  1:27059/1‑100 (MQ=255)
                                                                      tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                             <  2:380136/100‑1 (MQ=255)
                                                                      tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                             <  2:426078/100‑1 (MQ=255)
                                                                       ccGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                             >  2:2351460/1‑99 (MQ=255)
                                                                        cGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCag                           >  2:66193/1‑100 (MQ=255)
                                                                        cGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCag                           >  2:73515/1‑100 (MQ=255)
                                                                           aCTCGCACATAACACGCCCCCGGTTTGCACACCGGGAAGCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCagag                         <  2:793767/99‑1 (MQ=255)
                                                                            cTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcg                       <  2:1002303/100‑1 (MQ=255)
                                                                              cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                     >  1:340940/1‑100 (MQ=255)
                                                                              cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                     <  1:2097262/100‑1 (MQ=255)
                                                                                  caTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCTAGAGCGCGCTCa                 <  2:1702957/100‑1 (MQ=255)
                                                                                   ctcctggattcccGGTGTGCAAACCGTGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCAt                >  2:1577103/11‑100 (MQ=38)
                                                                                    taataCGCCCCCGTTTTCCACACCGGGAATCCAGGAGATTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATg               <  1:1340849/100‑1 (MQ=255)
                                                                                    taataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATg               >  1:1360227/1‑100 (MQ=255)
                                                                                     cctggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa              >  1:2053865/9‑100 (MQ=38)
                                                                                     aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa              >  2:1557246/1‑100 (MQ=255)
                                                                                     aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa              <  2:1043764/100‑1 (MQ=255)
                                                                                      ctggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg             >  1:1288822/8‑100 (MQ=255)
                                                                                      ataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg             >  2:360402/1‑100 (MQ=255)
                                                                                      ataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg             <  1:1177238/100‑1 (MQ=255)
                                                                                      ataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg             <  1:360400/100‑1 (MQ=255)
                                                                                        aCGCCCCCGGTTTGCACAGCGGGAATCCAGGAGAGTTCATGTACTACCCCGTCGTTCGTAAAGACCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt           <  2:455193/100‑1 (MQ=255)
                                                                                        aCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt           <  2:455143/100‑1 (MQ=255)
                                                                                         ngCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt          >  1:1544830/2‑100 (MQ=255)
                                                                                         cGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGCGttt          >  1:76588/1‑100 (MQ=255)
                                                                                          attcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa         >  1:115964/4‑100 (MQ=255)
                                                                                          attcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa         >  2:1532774/4‑100 (MQ=255)
                                                                                          attcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGACCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa         >  1:1796841/4‑100 (MQ=255)
                                                                                           ttcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc        >  1:1679109/3‑100 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc        >  1:973623/1‑100 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc        >  2:1114631/1‑100 (MQ=255)
                                                                                            ccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACt       >  2:1113259/1‑100 (MQ=255)
                                                                                                ggTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc   <  1:1919613/100‑1 (MQ=255)
                                                                                                ggTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc   >  2:2226886/1‑100 (MQ=255)
                                                                                                 gTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCa  <  2:1425009/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCA  >  NC_000913/1004642‑1004838

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