Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,741 T→G 15.5% 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,7410TG15.5% 340.7 / 55.1 136intergenic (+84/‑27)ycbF/pyrDputative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base T (54/60);  new base G (12/9);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 145 total aligned reads displayed.

GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAG  >  NC_000913/1004642‑1004839
                                                                                                   |                                                                                                  
gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                                                                                                    >  1:3535529/1‑100 (MQ=255)
gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                                                                                                    >  1:1598489/1‑100 (MQ=255)
gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTCCCGCAGGGCAAGTCCCTTTTGGGCCGAACGCGCACATAATACGCCCCCGGtt                                                                                                    >  2:1813258/1‑100 (MQ=255)
 cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                   >  1:921730/1‑100 (MQ=255)
  aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTTCCGCAGTTCAATTCCCTTTTGGTACGAACTCGCACATAATACGCCCCCGGttt                                                                                                   >  2:3357418/1‑99 (MQ=255)
   aaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                 <  1:233301/100‑1 (MQ=255)
   aaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                 <  2:3014854/100‑1 (MQ=255)
   aaaTCTGACCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                 >  1:2200864/1‑100 (MQ=255)
    aaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGca                                                                                                >  1:4167265/1‑100 (MQ=255)
     aTCGGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                 <  1:3362809/98‑1 (MQ=255)
      tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcaca                                                                                              >  1:1385974/1‑100 (MQ=255)
      tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcaca                                                                                              <  1:3086394/100‑1 (MQ=255)
         gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCg                                                                                           >  1:3345283/1‑100 (MQ=255)
         gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCg                                                                                           >  2:406178/1‑100 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCTCAGGTCAATTCCCTTTTGGTTCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                                                                                       >  2:3600642/1‑100 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                       <  2:2700128/100‑3 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                       <  2:901672/100‑3 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                       <  1:1913503/100‑3 (MQ=255)
              gagaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTGTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                       <  2:3679773/96‑3 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAAt                                                                                      >  2:2388287/1‑100 (MQ=255)
               aaaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggcg                                                                                     <  2:3232036/100‑5 (MQ=255)
                aaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATcc                                                                                    >  2:177781/1‑100 (MQ=255)
                 aCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCa                                                                                   <  1:1628657/100‑1 (MQ=255)
                 aCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCa                                                                                   <  1:3092502/100‑1 (MQ=255)
                 aCAGGTTCGGAAAACGTTTGCGTGTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATcc                                                                                    <  2:794270/99‑1 (MQ=255)
                    ggTTCGGAAAACGTTTGCGTTTTTTTTTTTTTTTTTCAATTCCCTTTTGGTCCGAACTCTCTCATAATTCGCCCCCGGTTTGCCCACCGTGATTCCAGga                                                                                >  2:3424552/1‑100 (MQ=255)
                    ggTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAgg                                                                                 <  2:3296431/99‑1 (MQ=255)
                    ggTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGga                                                                                >  1:910410/1‑100 (MQ=255)
                      ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                              >  2:3132521/1‑100 (MQ=255)
                      ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                              <  1:1987985/100‑1 (MQ=255)
                         ggAAAACGTTTGCGTTTTTTTTGCTGCAGTTAAATTCCCTTTTGGTCCGAACTCGCACATAATACGCGCCCGGTTTTCACACCGTGAATCCAGGAGAGtt                                                                           >  2:524559/1‑100 (MQ=255)
                         ggAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGtt                                                                           >  1:2767067/1‑100 (MQ=255)
                         ggAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGAGAATCCAGGAGAGtt                                                                           >  1:4137279/1‑100 (MQ=255)
                          gAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTc                                                                          <  1:263052/100‑1 (MQ=255)
                              aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGt                                                                      <  1:520498/100‑1 (MQ=255)
                              aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGt                                                                      <  2:3960422/100‑1 (MQ=255)
                               cGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGta                                                                     <  1:853145/100‑1 (MQ=255)
                                gTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac                                                                    >  1:237922/1‑100 (MQ=255)
                                     cGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcc                                                                >  2:2809161/1‑99 (MQ=255)
                                       ttttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCt                                                             <  2:4159650/100‑1 (MQ=255)
                                        tttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCtt                                                            <  1:1000500/100‑1 (MQ=255)
                                        tttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCtt                                                            <  1:1011558/100‑1 (MQ=255)
                                        tttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCtt                                                            >  1:4062128/1‑100 (MQ=255)
                                        ttgttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCtt                                                            <  2:760217/97‑1 (MQ=255)
                                           ttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                         >  1:2095615/1‑100 (MQ=255)
                                           ttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcg                                                          <  1:2008103/99‑1 (MQ=255)
                                            tttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                         <  2:2074651/99‑1 (MQ=255)
                                              tGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                      <  1:1766028/100‑1 (MQ=255)
                                              tGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                      <  1:2893228/100‑1 (MQ=255)
                                              tGCCGCAGGGCAATGCCCTTTTGTTCCGAACTCGCACATAATCCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttc                                                       <  1:1813262/99‑1 (MQ=255)
                                                ccGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                    >  1:3445445/1‑100 (MQ=255)
                                                 cGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCTTTCGTaa                                                   >  1:3421898/1‑100 (MQ=255)
                                                   cAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAg                                                 <  2:3894025/100‑1 (MQ=255)
                                                    aGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGc                                                <  1:3916522/100‑1 (MQ=255)
                                                      gTCAATTCCCTTTTGGTCCGAACTCGCACATAAGACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGccc                                              <  2:152018/100‑1 (MQ=255)
                                                        cAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCtt                                            >  2:1575439/1‑100 (MQ=255)
                                                          aTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCtttt                                          >  2:1743860/1‑100 (MQ=255)
                                                            tCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTcc                                        <  2:1598485/100‑1 (MQ=255)
                                                            tCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTcc                                        <  2:3378553/100‑1 (MQ=255)
                                                             cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                       >  1:3960134/1‑100 (MQ=255)
                                                             cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                       >  1:505956/1‑100 (MQ=255)
                                                              ccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAg                                      <  1:177781/100‑1 (MQ=255)
                                                               cgtttGGTCCGAACTCGCACGTAATACGCGCCCGGGTTGCACACCTGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTccn                                       <  2:2834938/96‑2 (MQ=255)
                                                               cTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGc                                     <  2:747952/100‑1 (MQ=255)
                                                               cTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGc                                     <  1:3613928/100‑1 (MQ=255)
                                                                 tttGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                    <  1:6241/99‑1 (MQ=255)
                                                                   tGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGa                                 >  2:3019669/1‑100 (MQ=255)
                                                                    ggTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGAt                                >  1:4083987/1‑100 (MQ=255)
                                                                     gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGACTCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCTGCTCGAt                                >  2:3029226/1‑99 (MQ=255)
                                                                     gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCCGCTCGATc                               >  1:1813985/1‑100 (MQ=255)
                                                                      tCCGAACTCGCACATAATACCCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                              >  2:109595/1‑100 (MQ=255)
                                                                       ccGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCa                             >  1:757696/1‑100 (MQ=255)
                                                                         gAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                           <  2:561043/100‑1 (MQ=255)
                                                                          aaCTCGCACATAATGCGCCCCCGGTATGCACAACGGTAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                           <  2:3215388/99‑1 (MQ=255)
                                                                            cccGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcg                        <  1:2165801/98‑1 (MQ=255)
                                                                              cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                      >  1:3287888/1‑100 (MQ=255)
                                                                              cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                      >  2:2318574/1‑100 (MQ=255)
                                                                                    taataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATg                <  2:3313330/100‑1 (MQ=255)
                                                                                     aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa               >  1:3871517/1‑100 (MQ=255)
                                                                                      ataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg              <  2:3795943/100‑1 (MQ=255)
                                                                                       tggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg              >  1:2893263/7‑99 (MQ=255)
                                                                                       tggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg              >  2:2713624/7‑99 (MQ=255)
                                                                                       tggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg              >  2:3969473/7‑99 (MQ=255)
                                                                                       taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg              <  2:1806519/99‑1 (MQ=255)
                                                                                        ggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt            >  2:2136219/6‑100 (MQ=255)
                                                                                        aCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACGCCTGCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGt             <  2:2385702/99‑1 (MQ=255)
                                                                                        aCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt            <  2:481316/100‑1 (MQ=255)
                                                                                         gatctccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt           <  2:1425082/95‑1 (MQ=255)
                                                                                         cGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt           <  1:894126/100‑1 (MQ=255)
                                                                                          gCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa          >  2:2864832/1‑100 (MQ=255)
                                                                                          attcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa          >  1:113442/4‑100 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc         <  1:1701228/100‑1 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc         >  2:2292768/1‑100 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc         <  1:1743864/100‑1 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc         >  1:2361283/1‑100 (MQ=255)
                                                                                           cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc         >  2:1543153/1‑100 (MQ=255)
                                                                                             cccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACtt       >  1:1686631/1‑100 (MQ=255)
                                                                                              ccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTagatc      >  1:1425085/1‑96 (MQ=255)
                                                                                                ggTTTGCACACCGGGAATCCAGGGGAGTGCATGGACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc    <  2:1493665/100‑1 (MQ=255)
                                                                                                 gTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCa   <  2:1999547/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAG  >  NC_000913/1004642‑1004839

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: