Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,746 C→A 5.3% intergenic (+89/‑22) 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,7460CA 5.3% 259.2 / 4.7 76intergenic (+89/‑22)ycbF/pyrDputative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base A (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCGGCAGAAAGTG  >  NC_000913/1004603‑1004893
                                                                                                                                               |                                                                                                                                                   
ggcgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCgg                                                                                                                                                <  2:749390/149‑1 (MQ=255)
  cgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGa                                                                                                                                              <  2:390247/149‑1 (MQ=255)
   gATAAATGATTACGGCGGGTTGAGTGCCAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                                                                                                                                             >  3:53805/1‑149 (MQ=255)
   gATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                                                                             <  2:696395/149‑3 (MQ=255)
    aTAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                                                                            <  3:277080/149‑4 (MQ=255)
     tAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTTCACACCGGGAATc                                                                                                                                           >  1:391281/1‑149 (MQ=255)
     tAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATc                                                                                                                                           >  2:373562/1‑149 (MQ=255)
       aaTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCa                                                                                                                                         >  1:410617/1‑149 (MQ=255)
        aTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAg                                                                                                                                        <  3:613946/149‑1 (MQ=255)
         tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAgg                                                                                                                                       >  3:158876/1‑149 (MQ=255)
          gATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGga                                                                                                                                      <  4:262095/149‑1 (MQ=255)
            ttACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                                                                                    <  1:216724/149‑1 (MQ=255)
            ttACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                                                                                    >  3:303239/1‑149 (MQ=255)
            ttACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGATCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTTCACACCTGGAAACCAGgaga                                                                                                                                    >  2:137019/1‑149 (MQ=255)
                  cggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                                                                              <  1:693620/149‑1 (MQ=255)
                  cggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                                                                              <  3:248026/149‑1 (MQ=255)
                  cggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                                                                              <  1:304454/149‑1 (MQ=255)
                     gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGta                                                                                                                           <  1:279911/149‑1 (MQ=255)
                         aGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtactac                                                                                                                       <  4:289669/149‑1 (MQ=255)
                           tGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAccc                                                                                                                     >  1:435245/1‑149 (MQ=255)
                            gCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcccc                                                                                                                    >  1:159592/1‑149 (MQ=255)
                               aagaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttc                                                                                                                 <  3:61103/149‑1 (MQ=255)
                                 gaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                                                                               <  1:331631/149‑1 (MQ=255)
                                 gaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgt                                                                                                               >  3:225174/1‑149 (MQ=255)
                                     gAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTtcgttcgt                                                                                                           >  4:489543/1‑149 (MQ=255)
                                       gCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaa                                                                                                         >  2:180594/1‑149 (MQ=255)
                                        cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTAGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAGTCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaaa                                                                                                        >  2:219015/1‑149 (MQ=255)
                                          aaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACTCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGc                                                                                                      <  4:101111/149‑1 (MQ=255)
                                            aTCTGCCCTGAAACAGGTTCGGAGAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTTCCGTTTGGTCCGAACGCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGTAGAGTTCATGTACTACCCCTTGGTTCGTAAAGccc                                                                                                    <  4:476259/149‑1 (MQ=255)
                                            aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGccc                                                                                                    <  4:240213/149‑1 (MQ=255)
                                             tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCt                                                                                                   >  3:136328/1‑149 (MQ=255)
                                               tGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCttt                                                                                                 >  1:447499/1‑149 (MQ=255)
                                                     gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGc                                                                                           <  4:6603/149‑1 (MQ=255)
                                                      aaaCAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATTCGCCCCCGGTTTTCACACCGGGAATCCAGGACAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTGCAGCt                                                                                          <  4:481760/149‑1 (MQ=255)
                                                          aGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGAt                                                                                      >  3:663342/1‑149 (MQ=255)
                                                             ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                                                                              <  4:289954/106‑1 (MQ=255)
                                                             ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                                                                                                                              >  3:289954/1‑106 (MQ=255)
                                                             ttCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCa                                                                                   >  4:156992/1‑149 (MQ=255)
                                                                     aCGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcgc                                                                           <  1:379140/149‑1 (MQ=255)
                                                                           gCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa                                                                     <  2:606431/149‑1 (MQ=255)
                                                                           gCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGa                                                                     <  1:658733/149‑1 (MQ=255)
                                                                             gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGt                                                                   >  2:84289/1‑149 (MQ=255)
                                                                              ttttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt                                                                  <  4:158876/149‑1 (MQ=255)
                                                                              gtttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCt                                                                          <  4:284391/140‑1 (MQ=255)
                                                                              gtttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCt                                                                          >  3:284391/2‑141 (MQ=255)
                                                                                ttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTa                                                                <  3:616691/149‑1 (MQ=255)
                                                                                      gCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGATATGCACACCGGGAATCCAGAAGAGTTCATGTACTACCCCTACGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTAATTTTc                                                          <  1:759712/149‑1 (MQ=255)
                                                                                          cAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCa                                                      <  1:84289/149‑1 (MQ=255)
                                                                                             gTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAAtt                                                   >  4:399753/1‑149 (MQ=255)
                                                                                                aaTTCCCTTTTGGTCCGAACTCGCATATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCCTGAGTTTACTTTTCAGCAATTACg                                                >  3:182846/1‑149 (MQ=255)
                                                                                                aaTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACg                                                >  2:453611/1‑149 (MQ=255)
                                                                                                 aTTCCCTTTTGGTCAGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGc                                               >  3:53947/1‑149 (MQ=255)
                                                                                                     ccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTa                                           <  1:217822/149‑1 (MQ=255)
                                                                                                           ggTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAg                                     <  1:137019/149‑1 (MQ=255)
                                                                                                             tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGa                                   >  4:46409/1‑149 (MQ=255)
                                                                                                             tCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGa                                   >  2:740685/1‑149 (MQ=255)
                                                                                                                      gCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGttt                          >  3:123430/1‑149 (MQ=255)
                                                                                                                      gCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGttt                          >  1:215525/1‑149 (MQ=255)
                                                                                                                      gCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGttt                          <  3:247560/149‑1 (MQ=255)
                                                                                                                         caTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGaa                       <  3:489543/149‑1 (MQ=255)
                                                                                                                         caTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGaa                       <  1:206351/149‑1 (MQ=255)
                                                                                                                         caTAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGaa                       <  2:435245/149‑1 (MQ=255)
                                                                                                                              taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACt                  >  1:542211/1‑149 (MQ=255)
                                                                                                                               aCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTg                 <  4:123430/149‑1 (MQ=255)
                                                                                                                                gattcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTgg                >  2:443516/5‑149 (MQ=255)
                                                                                                                                cGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTgg                <  4:406027/149‑1 (MQ=255)
                                                                                                                                 gCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGt               <  4:400910/149‑1 (MQ=255)
                                                                                                                                  cccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTg              <  2:657799/149‑1 (MQ=255)
                                                                                                                                   tcccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGc             >  3:428353/2‑149 (MQ=255)
                                                                                                                                   ccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGc             <  3:494224/149‑1 (MQ=255)
                                                                                                                                    cccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCCGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCg            >  3:515898/1‑149 (MQ=255)
                                                                                                                                    cccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt                                                                 >  3:236354/1‑96 (MQ=255)
                                                                                                                                    cccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGttt                                                                 <  4:236354/96‑1 (MQ=255)
                                                                                                                                    cccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCg            >  4:363533/1‑149 (MQ=255)
                                                                                                                                     ccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCgg           >  3:612589/1‑149 (MQ=255)
                                                                                                                                     ccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCgg           >  2:372039/1‑149 (MQ=255)
                                                                                                                                              acacCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCGGCAGAAAGTg  >  2:76790/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGCGGCAGAAAGTG  >  NC_000913/1004603‑1004893

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