Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,741 T→G 15.3% 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.3% 386.8 / 66.3 162intergenic (+84/‑27)ycbF/pyrDputative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base T (53/84);  new base G (14/11);  total (67/95)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-02

Only 100 of 175 total aligned reads displayed.

GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGC  >  NC_000913/1004642‑1004840
                                                                                                     |                                                                                                   
 cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                    <  2:1991692/100‑1 (MQ=255)
 cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                    <  1:744473/100‑1 (MQ=255)
  aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                    >  2:1231686/1‑100 (MQ=255)
  aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                    >  2:3555602/1‑99 (MQ=255)
   aaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                  <  2:3384172/100‑1 (MQ=255)
   aaaTCTGCCCCGAAACAGGTTCGGAAAACGTTTGCGGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                                                                                                    <  1:1295220/99‑1 (MQ=255)
    aaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGca                                                                                                 <  2:588157/100‑1 (MQ=255)
     aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGTTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcac                                                                                                >  1:1708715/1‑100 (MQ=255)
     aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                                                                                                  <  2:2633089/98‑1 (MQ=255)
       cTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACAcc                                                                                             <  2:1046012/100‑1 (MQ=255)
         gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCg                                                                                            >  1:2424377/1‑100 (MQ=255)
          cccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTTTTTTTTTTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCacac                                                                                              >  1:2176001/1‑97 (MQ=255)
          cccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgg                                                                                           >  1:1912734/1‑100 (MQ=255)
           ccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCggg                                                                                          >  1:66654/1‑100 (MQ=255)
           ccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCggg                                                                                          <  1:1021928/100‑1 (MQ=255)
             tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCgg                                                                                           <  1:1298459/97‑1 (MQ=255)
              gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                       <  2:88973/100‑4 (MQ=255)
               aaaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggcg                                                                                      <  2:1912741/100‑5 (MQ=255)
               aaaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggcg                                                                                      <  1:1367596/100‑5 (MQ=255)
                aaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGGAGGTCAATTTCCTTTGGGTCCGAACGCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                       <  1:1833316/98‑4 (MQ=255)
                  cAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCGCCCGGTTTGCACACCGGGAATCCAg                                                                                   >  1:1601590/1‑100 (MQ=255)
                  cAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAg                                                                                   <  2:1964070/100‑1 (MQ=255)
                     gTTCGGAAAACGTTTGCGTTTTTTTTTTTTTTTTTTAATTCCCTTTTGGCCCGAACTCTCACATATTACGCCCCCGGTTTGCCCACCGGGAATc                                                                                      >  1:3480473/1‑94 (MQ=255)
                     gTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgag                                                                                <  2:2540763/100‑1 (MQ=255)
                     gTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgag                                                                                <  1:1489290/100‑1 (MQ=255)
                      ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                               >  2:1128202/1‑100 (MQ=255)
                      ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                               <  2:2965733/100‑1 (MQ=255)
                           aaaaCGTTTGCGTTTTTTTTGCTGCAGGTCAATTACCTGTTGGTCCGAACTCTCACATAATACGCACCCGCTTTGCACACCGCGAATCCAGGAGCGTTCa                                                                          >  1:2635616/1‑100 (MQ=255)
                           aaaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa                                                                          <  2:1619840/100‑1 (MQ=255)
                           aaaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa                                                                          <  1:4113126/100‑1 (MQ=255)
                           aaaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa                                                                          <  2:3821404/100‑1 (MQ=255)
                            aaaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                         >  2:1766678/1‑100 (MQ=255)
                             aaCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATg                                                                        <  1:1104252/100‑1 (MQ=255)
                              aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGTTCCCGGTTTGCACATCGGGAATCTATGAGAGTTTAttt                                                                       <  1:360295/100‑3 (MQ=255)
                              aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGt                                                                       <  1:4299986/100‑1 (MQ=255)
                               cGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGta                                                                      <  1:3874104/100‑1 (MQ=255)
                               cGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGta                                                                      <  1:774086/100‑1 (MQ=255)
                                   tGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtactac                                                                  <  2:1541955/100‑1 (MQ=255)
                                    gCGTTTTTTTTTTTTCATGTTAATTCCATTTTGGTCCGAACTCGCACATAATATGCGCTCGGTTTTCACATCGGGAATCCAAGAGAGTTCATGttct                                                                    >  1:2868408/1‑97 (MQ=255)
                                    gCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcc                                                                 <  2:1296530/100‑1 (MQ=255)
                                       gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcccc                                                               >  2:3768140/1‑100 (MQ=255)
                                         ttttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCt                                                              <  2:3150413/100‑1 (MQ=255)
                                             ttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                          >  2:2784993/1‑100 (MQ=255)
                                             ttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                          <  2:1572836/100‑1 (MQ=255)
                                              tttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgtt                                                         <  1:1425552/100‑1 (MQ=255)
                                                  ccGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                     >  1:1444674/1‑100 (MQ=255)
                                                  ccGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                     >  2:2612705/1‑100 (MQ=255)
                                                   cGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaa                                                    >  2:4272612/1‑100 (MQ=255)
                                                   cGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                     >  1:1000090/1‑99 (MQ=255)
                                                    gCAGGTCAATTCGCTTTTGGTCCGAACGCGCACATAATACGGCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttcg                                                       <  1:1942426/96‑1 (MQ=255)
                                                    gCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaaa                                                   >  1:4028642/1‑100 (MQ=255)
                                                       ggTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGcc                                                <  2:3599274/100‑1 (MQ=255)
                                                        gTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGccc                                               <  2:3919272/100‑1 (MQ=255)
                                                         tCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCt                                              <  2:2640407/100‑1 (MQ=255)
                                                           aaTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCttt                                            >  2:84704/1‑100 (MQ=255)
                                                             ttCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCtttt                                           <  2:1331326/99‑1 (MQ=255)
                                                              tCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCttttac                                         <  2:2981200/100‑3 (MQ=255)
                                                               cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                        >  1:1423283/1‑100 (MQ=255)
                                                               cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                        <  1:1636083/100‑1 (MQ=255)
                                                               cccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                        <  2:1621897/100‑1 (MQ=255)
                                                                  ttttTGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                     <  1:1128198/100‑1 (MQ=255)
                                                                  ttttGGTCCGAACTCTCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                     <  2:4335878/100‑1 (MQ=255)
                                                                   tttGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTCCTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTc                                    >  1:2326703/1‑100 (MQ=255)
                                                                   tttGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                     <  1:598104/99‑1 (MQ=255)
                                                                    ttGGTCCGAACTCGCGCATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                     <  2:2012554/98‑1 (MQ=255)
                                                                        tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                               <  2:2176008/100‑1 (MQ=255)
                                                                        tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                               <  2:3357057/100‑1 (MQ=255)
                                                                        tCCGAACTCGCACATAATACACCCCCGGTTTGCACACCGAGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATcc                               <  1:2934753/100‑1 (MQ=255)
                                                                         ccGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCa                              >  2:565730/1‑100 (MQ=255)
                                                                         ccGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCa                              <  2:3620569/100‑1 (MQ=255)
                                                                          cGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCag                             >  1:443908/1‑100 (MQ=255)
                                                                           gAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                            >  2:247702/1‑100 (MQ=255)
                                                                           gAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAAGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                            <  2:20717/100‑1 (MQ=255)
                                                                            aaCTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCagag                           <  1:2808395/100‑1 (MQ=255)
                                                                            aaCTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCagag                           <  1:1360189/100‑1 (MQ=255)
                                                                                cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                       <  1:4029216/100‑1 (MQ=255)
                                                                                cGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcg                       <  2:226027/100‑1 (MQ=255)
                                                                                 gCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcgc                      <  1:804318/100‑1 (MQ=255)
                                                                                  cacaTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTACAGCCCTTTTCCAGCTCGATCCAgcgcgcgc                      >  1:3528811/1‑99 (MQ=255)
                                                                                  cacaTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCt                     >  1:328141/1‑100 (MQ=255)
                                                                                  cacaTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCATAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCt                     <  2:3737293/100‑1 (MQ=255)
                                                                                      taataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATg                 <  2:507403/100‑1 (MQ=255)
                                                                                        ataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAg               >  2:1687735/1‑100 (MQ=255)
                                                                                         tggattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGt              >  1:3775991/7‑100 (MQ=255)
                                                                                           gattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt            >  2:4171246/5‑100 (MQ=255)
                                                                                           gattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt            >  1:3973009/5‑100 (MQ=255)
                                                                                           cGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt            <  2:3663596/100‑1 (MQ=255)
                                                                                            gCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa           >  2:3846962/1‑100 (MQ=255)
                                                                                            attcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa           >  2:3693822/4‑100 (MQ=255)
                                                                                             ttcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc          >  2:2574627/3‑100 (MQ=255)
                                                                                             ttcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTAATACCCCATCGTTCGTAAAGCCCTTTTCCAGCTAGATCCAGAGAGCGCTCATGCGtt             >  1:3433922/3‑97 (MQ=255)
                                                                                             cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAc          >  2:2401428/1‑100 (MQ=255)
                                                                                              tcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGACCTTTTCCAGCTCGATCGAGGGCGCGCTCCTGAGTTTACt         >  2:1667754/2‑100 (MQ=255)
                                                                                               cccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACtt        >  2:2335909/1‑100 (MQ=255)
                                                                                               cccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACtt        >  1:694405/1‑100 (MQ=255)
                                                                                                  ggTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc     <  1:3754234/100‑1 (MQ=255)
                                                                                                   gTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc     <  1:1422687/99‑1 (MQ=255)
                                                                                                    tttGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAg   <  2:3809582/100‑1 (MQ=255)
                                                                                                    tttGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAg   <  2:2860559/100‑1 (MQ=255)
                                                                                                     ttGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGc  <  1:2310073/100‑1 (MQ=255)
                                                                                                     |                                                                                                   
GCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGC  >  NC_000913/1004642‑1004840

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