Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,997,902 A→G intergenic (+3/‑81) yigB → / → uvrD 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase/ssDNA translocase and dsDNA helicase ‑ DNA helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,997,9020AG100.0% 255.5 / NA 69intergenic (+3/‑81)yigB/uvrD5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase/ssDNA translocase and dsDNA helicase ‑ DNA helicase II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/47);  total (22/47)

GGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGA  >  NC_000913/3997765‑3998029
                                                                                                                                         |                                                                                                                               
ggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGc                                                                                                                               >  1:963495/1‑140 (MQ=255)
ggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGc                                                                                                                               >  2:1062857/1‑140 (MQ=255)
 gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGTAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  2:872428/140‑1 (MQ=255)
 gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  1:960819/140‑1 (MQ=255)
 gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  1:448659/140‑1 (MQ=255)
 gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  1:1188722/140‑1 (MQ=255)
 gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  2:1032470/140‑1 (MQ=255)
 gtgggGCAATTCGCAGCGGAAGGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  2:265931/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:1001963/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             >  1:522516/1‑140 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:499920/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:294177/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:1062857/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:1115127/140‑1 (MQ=255)
  tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaa                                                                                                                             <  1:272982/140‑1 (MQ=255)
   ggggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCaaa                                                                                                                            <  1:87737/140‑1 (MQ=255)
    gggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAAt                                                                                                                           <  2:465667/140‑1 (MQ=255)
    gggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAAt                                                                                                                           <  1:709499/140‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCACACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              <  2:974301/136‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCa                                                                                                                              >  1:974301/1‑136 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          <  2:209278/140‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          <  1:950349/140‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          >  1:822136/1‑140 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          <  1:531372/140‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          >  1:515234/1‑140 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          <  2:1170111/140‑1 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          >  1:222683/1‑140 (MQ=255)
     ggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATc                                                                                                                          <  1:243536/140‑1 (MQ=255)
      gCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCt                                                                                                                         <  2:522516/140‑1 (MQ=255)
      gCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCt                                                                                                                         <  2:479183/140‑1 (MQ=255)
      gCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCt                                                                                                                         <  2:732560/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  1:916355/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  2:845955/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  2:848420/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  2:1141193/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  2:963495/140‑1 (MQ=255)
       cAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  2:161528/140‑1 (MQ=255)
                    aaTGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAg                                                                                                           <  1:86304/140‑1 (MQ=255)
                      tGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTTACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        >  2:1127384/1‑125 (MQ=255)
                      tGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTg                                                                                                                        <  1:1127384/125‑1 (MQ=255)
                           gCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTg                                                                                                    >  2:1010835/1‑140 (MQ=255)
                               gttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCgg                                                                                                <  2:822136/140‑1 (MQ=255)
                                     tCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCg                                                                                          <  2:228693/140‑1 (MQ=255)
                                         aCCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTgc                                                                                      >  1:607858/1‑140 (MQ=255)
                                              aaaaTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTtc                                                                                 <  2:889195/140‑1 (MQ=255)
                                                  tGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCg                                                                             <  2:236051/140‑1 (MQ=255)
                                                        tCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAAc                                                                       <  2:607858/140‑1 (MQ=255)
                                                            aTGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAt                                                                   >  2:148177/1‑140 (MQ=255)
                                                                    cTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAcgc                                                           >  2:101472/1‑140 (MQ=255)
                                                                      gggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACgcgg                                                         >  2:564252/1‑140 (MQ=255)
                                                                        gACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACgcggcg                                                       <  1:101472/140‑1 (MQ=255)
                                                                            gCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGc                                                   >  1:1148927/1‑140 (MQ=255)
                                                                                 ttACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATg                                              >  2:830588/1‑140 (MQ=255)
                                                                                   aCTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGa                                            <  2:1148927/140‑1 (MQ=255)
                                                                                   aCTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGa                                            <  1:276572/140‑1 (MQ=255)
                                                                                   aCTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTCTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGa                                            <  1:725959/140‑1 (MQ=255)
                                                                                       ccGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGtt                                        >  1:155151/1‑140 (MQ=255)
                                                                                       ccGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGtt                                        >  1:455426/1‑140 (MQ=255)
                                                                                                     ttCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGa                          <  1:973877/140‑1 (MQ=255)
                                                                                                      tCCCGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGAc                         <  1:564252/140‑1 (MQ=255)
                                                                                                       cccGGTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACa                        <  2:455426/140‑1 (MQ=255)
                                                                                                          ggTTGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGcc                     >  1:691781/1‑140 (MQ=255)
                                                                                                             tGGCATCTCTGACCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTa                  <  2:447405/140‑1 (MQ=255)
                                                                                                                        aCCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt                    >  1:1027029/1‑127 (MQ=255)
                                                                                                                        aCCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt                    <  2:1027029/127‑1 (MQ=255)
                                                                                                                        aCCTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAg       >  2:423470/1‑140 (MQ=255)
                                                                                                                         ccTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt                    <  1:217624/126‑1 (MQ=255)
                                                                                                                         ccTCGCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt                    >  2:217624/1‑126 (MQ=255)
                                                                                                                             gCTGATATAATCGGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGa  >  1:966320/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                               
GGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGA  >  NC_000913/3997765‑3998029

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