Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,757,963 (T)7→8 intergenic (‑136/+54) yiaY ← / ← selB L‑threonine dehydrogenase/selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,9561.T95.5% 2156.7 / 78.8 508intergenic (‑129/+61)yiaY/selBL‑threonine dehydrogenase/selenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base . (9/14);  new base T (240/245);  total (249/259)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 548 total aligned reads displayed.

CGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATA  >  NC_000913/3757856‑3758051
                                                                                                     |                                                                                                
 gggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑t                                                                                                >  2:7268527/1‑101 (MQ=255)
     aGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑ttttt                                                                                            <  1:7268123/101‑1 (MQ=255)
      gCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑tttttt                                                                                           >  1:7268531/1‑101 (MQ=255)
       cTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑ttttttt                                                                                          >  2:7268535/1‑101 (MQ=255)
        gtACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑tttttttt                                                                                         >  1:7268543/2‑100 (MQ=255)
          aCTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTtaa                                                                                       >  2:7268548/1‑101 (MQ=255)
          aCTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTtaa                                                                                       <  2:7267774/101‑1 (MQ=255)
            tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTtaaat                                                                                     <  1:7267890/101‑3 (MQ=255)
            tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTtaaat                                                                                     <  2:7267951/101‑3 (MQ=255)
             gAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAAtt                                                                                     >  1:7268554/1‑101 (MQ=255)
               ggATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCa                                                                                   >  1:7268562/1‑101 (MQ=255)
               ggATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCa                                                                                   >  2:7268559/1‑101 (MQ=255)
                gATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCAt                                                                                  <  1:7267830/101‑1 (MQ=255)
                gATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCAt                                                                                  <  2:7267669/101‑1 (MQ=255)
                 aTTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATa                                                                                 >  1:7268572/1‑101 (MQ=255)
                  ttttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATaa                                                                                <  1:7267975/101‑1 (MQ=255)
                   tttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAg                                                                               >  2:7268575/1‑101 (MQ=255)
                    ttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGa                                                                              <  1:7268142/101‑1 (MQ=255)
                     tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGaa                                                                             >  1:7268582/1‑101 (MQ=255)
                      cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAAt                                                                            <  2:7267972/101‑1 (MQ=255)
                      cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAAt                                                                            >  2:7268584/1‑101 (MQ=255)
                        ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAAttt                                                                          >  1:7268593/1‑101 (MQ=255)
                          cGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTg                                                                        >  2:7268598/1‑101 (MQ=255)
                           gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGt                                                                       >  1:7268603/1‑101 (MQ=255)
                            aaaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTa                                                                      <  1:7268066/101‑1 (MQ=255)
                             aaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTAt                                                                     >  2:7268608/1‑101 (MQ=255)
                              aaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTAtt                                                                    >  1:7268616/1‑101 (MQ=255)
                                 atGATCTGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAat                                                                 >  2:7268621/1‑101 (MQ=255)
                                 atGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTTAAATTCATAAGAATTTTGTATTaa                                                                  <  1:7267993/101‑1 (MQ=255)
                                  tGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAatt                                                                <  2:7267824/101‑1 (MQ=255)
                                   gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAattt                                                               >  1:7268627/1‑101 (MQ=255)
                                   gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAattt                                                               >  2:7268629/1‑101 (MQ=255)
                                   gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAattt                                                               <  1:7267769/101‑1 (MQ=255)
                                    aTCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAattta                                                              >  1:7268633/1‑101 (MQ=255)
                                     tCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAatttat                                                             <  2:7268103/101‑1 (MQ=255)
                                       aGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAatttattt                                                           >  2:7268641/1‑101 (MQ=255)
                                       aGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAatttattt                                                           <  2:7267874/101‑1 (MQ=255)
                                         tattatTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTa                                                         >  1:7268647/1‑101 (MQ=255)
                                          attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                        <  1:7268239/101‑1 (MQ=255)
                                          attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                        <  1:7268023/101‑1 (MQ=255)
                                          attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                        >  2:7268651/1‑101 (MQ=255)
                                          attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                        <  1:7267469/101‑1 (MQ=255)
                                           ttattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                       <  1:7268086/101‑1 (MQ=255)
                                           ttattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                       <  1:7267915/101‑1 (MQ=255)
                                           ttattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                       <  1:7267747/101‑1 (MQ=255)
                                            tattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaaa                                                      >  1:7268657/1‑101 (MQ=255)
                                             attGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                     >  2:7268664/1‑101 (MQ=255)
                                              ttGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACt                                                    <  2:7268060/101‑1 (MQ=255)
                                                gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtta                                                  >  2:7268673/1‑101 (MQ=255)
                                                gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtta                                                  >  1:7268676/1‑101 (MQ=255)
                                                 aTTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaa                                                 <  2:7267722/101‑1 (MQ=255)
                                                 aTTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaa                                                 <  2:7267914/101‑1 (MQ=255)
                                                  ttATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                <  2:7268148/101‑1 (MQ=255)
                                                  ttATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                <  2:7268353/101‑1 (MQ=255)
                                                  ttATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                >  2:7268682/1‑101 (MQ=255)
                                                   tATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAattta                                                              <  2:7268077/86‑1 (MQ=255)
                                                    aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaatta                                              <  1:7267814/101‑1 (MQ=255)
                                                     tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaattaa                                             >  2:7268691/1‑101 (MQ=255)
                                                      cAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTtaattaat                                            >  1:7268693/1‑101 (MQ=255)
                                                       aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATc                                           >  1:7268701/1‑101 (MQ=255)
                                                        aaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCa                                          <  1:7267694/101‑1 (MQ=255)
                                                          aTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCAtt                                        >  2:7268708/1‑101 (MQ=255)
                                                           ttAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTc                                      >  1:7268709/1‑101 (MQ=255)
                                                              aTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                    >  1:7268718/1‑101 (MQ=255)
                                                               tCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCtt                                   <  1:7267792/101‑1 (MQ=255)
                                                               tCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCtt                                   <  1:7268051/101‑1 (MQ=255)
                                                                 taataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAt                                 <  2:7268252/101‑1 (MQ=255)
                                                                  aataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                <  2:7268032/101‑1 (MQ=255)
                                                                  aataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                >  2:7268721/1‑101 (MQ=255)
                                                                    taaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtttt                              <  1:7268338/101‑1 (MQ=255)
                                                                    taaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtttt                              <  2:7268313/101‑1 (MQ=255)
                                                                     aaaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttttt                             <  2:7268190/101‑1 (MQ=255)
                                                                      aaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTc                            >  1:7268732/1‑101 (MQ=255)
                                                                      aaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTc                            >  2:7268730/1‑101 (MQ=255)
                                                                        aaGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                          >  1:7268741/1‑101 (MQ=255)
                                                                        aaGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                          >  2:7268740/1‑101 (MQ=255)
                                                                         aGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                         <  2:7267847/101‑1 (MQ=255)
                                                                          gCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                        <  1:7267864/101‑1 (MQ=255)
                                                                          gCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                        <  2:7267609/101‑1 (MQ=255)
                                                                          gCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                        <  2:7268171/101‑1 (MQ=255)
                                                                           cacaTTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttt                               >  2:7268752/1‑93 (MQ=255)
                                                                            acaTTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGaaa                      >  1:7268753/1‑101 (MQ=255)
                                                                              aTTCTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaata                    >  2:7268758/1‑101 (MQ=255)
                                                                                tGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataat                  <  1:7268097/101‑1 (MQ=255)
                                                                                tGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataat                  >  1:7268764/1‑101 (MQ=255)
                                                                                    ttCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATg              <  1:7268668/101‑1 (MQ=255)
                                                                                      cATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCg            <  1:7268408/101‑1 (MQ=255)
                                                                                      cATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCg            <  1:7268283/101‑1 (MQ=255)
                                                                                       aTTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGt           >  2:7268774/1‑101 (MQ=255)
                                                                                       aTTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGt           >  1:7268775/1‑101 (MQ=255)
                                                                                       aTTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGt           <  2:7268519/101‑1 (MQ=255)
                                                                                         tAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTcg         >  1:7268782/1‑101 (MQ=255)
                                                                                         tAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTcg         <  1:7268167/101‑1 (MQ=255)
                                                                                          aaaCAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTcgc        <  2:7267996/101‑1 (MQ=255)
                                                                                             cAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCa     >  2:7268788/1‑101 (MQ=255)
                                                                                              aaTGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa    >  1:7268792/1‑101 (MQ=255)
                                                                                              aaTGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa    <  1:7268216/101‑1 (MQ=255)
                                                                                              aaTGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa    <  1:7267948/101‑1 (MQ=255)
                                                                                                tgtgCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATa  <  2:7268120/101‑1 (MQ=255)
                                                                                                tgtgCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATa  <  1:7268194/101‑1 (MQ=255)
                                                                                                     |                                                                                                
CGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGC‑‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATA  >  NC_000913/3757856‑3758051

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