Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,757,942 C→A 5.4% intergenic (‑115/+75) yiaY ← / ← selB L‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,9420CA 5.4% 619.7 / 14.8 193intergenic (‑115/+75)yiaY/selBL‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base C (94/88);  new base A (7/4);  total (101/92)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 195 total aligned reads displayed.

ACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCG  >  NC_000913/3757796‑3758088
                                                                                                                                                   |                                                                                                                                                   
   gAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTaa                                                                                                                                                <  1:1067160/149‑1 (MQ=255)
    aattaaTAAAGAACGTTGAAGCTGCCATAACACACTAATGAAATTTTAAAAAAGGGCGTGCAGATTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTaaa                                                                                                                                              <  1:178643/145‑1 (MQ=255)
    aaGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTaaa                                                                                                                                              >  1:1118678/1‑149 (MQ=255)
      ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACa                                                                                                                                            >  2:1514461/1‑149 (MQ=255)
       gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                           >  2:852822/1‑149 (MQ=255)
          tAAAGAACGTTGATGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                        >  2:649394/1‑149 (MQ=255)
          tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                        <  2:1121334/149‑1 (MQ=255)
          tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                        <  2:467693/149‑1 (MQ=255)
            aaGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGc                                                                                                                                      >  2:354727/1‑149 (MQ=255)
            aaGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGc                                                                                                                                      >  1:834232/1‑149 (MQ=255)
             aGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTTTAATGAATACAATGTGCt                                                                                                                                     <  1:1248363/149‑1 (MQ=255)
              taaCGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAATGGGCGGGCAGCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCtt                                                                                                                                    <  1:263243/148‑1 (MQ=255)
              gAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAACAATGTGCttt                                                                                                                                   >  1:1574794/1‑149 (MQ=255)
                 cGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttt                                                                                                                                 <  1:979437/149‑1 (MQ=255)
                  gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCtttttt                                                                                                                                >  1:69267/1‑149 (MQ=255)
                    agtaagctgcccgcccttttttacactttcaggagtgtgttatggcagcttACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                        >  1:1148906‑M2/52‑140 (MQ=255)
                      taagctgcccgcccttttttacactttcaggagtgtgttatggcagcttACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATATAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaa                                                                                                                             >  2:1272768‑M2/50‑149 (MQ=255)
                      AAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaaa                                                                                                                            >  1:1034429‑M1/1‑149 (MQ=255)
                      AAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaaa                                                                                                                            >  1:25001‑M1/1‑149 (MQ=255)
                         CTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTc                                                                                                                         <  2:1134137‑M1/149‑1 (MQ=255)
                         CTGCCATAACACACTCCTGAAAGTGTAAAAAACGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCAATAAACAAAGAGATTTATTTAAatta                                                                                                                         <  2:1596213‑M1/149‑2 (MQ=255)
                           GCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCAt                                                                                                                       >  2:871639‑M1/1‑149 (MQ=255)
                               tAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGa                                                                                                                   >  1:1537905/1‑149 (MQ=255)
                                   acacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAttt                                                                                                               >  2:759528/1‑149 (MQ=255)
                                    cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAtttt                                                                                                              <  2:272377/149‑1 (MQ=255)
                                    aacTCCTGAAAGTGTAAAAAAGTGCTGGCATCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAtttt                                                                                                              <  1:132107/148‑1 (MQ=255)
                                       tCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTa                                                                                                           <  2:1166400/149‑1 (MQ=255)
                                       tCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTa                                                                                                           >  1:879027/1‑149 (MQ=255)
                                           gAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTaa                                                                                                       <  2:1525269/149‑1 (MQ=255)
                                              aGTGTAAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatt                                                                                                     <  1:361779/149‑1 (MQ=255)
                                               gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATAGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAattta                                                                                                   <  2:949615/149‑1 (MQ=255)
                                                 gtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTTTTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttatt                                                                                                 <  2:1416379/149‑1 (MQ=255)
                                                 gtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttatt                                                                                                 <  2:1455371/149‑1 (MQ=255)
                                                  tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttattt                                                                                                >  2:493387/1‑149 (MQ=255)
                                                     aaaaaTGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTa                                                                                              <  2:1530871/149‑1 (MQ=255)
                                                      aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTc                                                                                                                                                     <  1:487572/94‑1 (MQ=255)
                                                        aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaaa                                                                                           <  2:32554/149‑1 (MQ=255)
                                                         aGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                          >  2:470171/1‑149 (MQ=255)
                                                             cgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaa                                                                                      <  1:1328914/149‑1 (MQ=255)
                                                                 cagcttACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaattaa                                                                                  <  2:262739‑M2/143‑1 (MQ=255)
                                                                   gcttACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATc                                                                                >  1:1387256‑M2/5‑149 (MQ=255)
                                                                     ttACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTAATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCAt                                                                              >  2:1484594‑M2/3‑149 (MQ=255)
                                                                         tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTcc                                                                          >  2:217974/1‑149 (MQ=255)
                                                                          gAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                                                         <  2:1537905/149‑1 (MQ=255)
                                                                            tgATTTTGATTCGAAAATATGATGAGTATTATTGATTATTAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTa                                                                       <  1:660193/148‑1 (MQ=255)
                                                                                  tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTc                                                                 <  2:604008/149‑1 (MQ=255)
                                                                                   cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                <  2:484225/149‑1 (MQ=255)
                                                                                     ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaaa                                                                                                                            >  2:544248/1‑87 (MQ=255)
                                                                                     ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaaa                                                                                                                            <  1:544248/87‑1 (MQ=255)
                                                                                     ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                              >  1:1191949/1‑149 (MQ=255)
                                                                                     ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                              >  2:1307170/1‑149 (MQ=255)
                                                                                     ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCAgg                                                              >  1:1045504/1‑149 (MQ=255)
                                                                                      tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAATCAATGTGCTTTTTTTAAATACATAAGAATTTTGTATTAATTTATTTTAAAACTTATTTAATCATTTCCTTATTTTTCCGGa                                                             >  2:225149/1‑149 (MQ=255)
                                                                                      tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAATTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                                                             <  1:295072/149‑1 (MQ=255)
                                                                                         aaaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAATAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCGGAaat                                                          >  1:1474528/1‑149 (MQ=255)
                                                                                         aaaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGTAaat                                                          >  2:207210/1‑149 (MQ=255)
                                                                                               tGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataat                                                    <  1:581218/149‑1 (MQ=255)
                                                                                               tGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataat                                                    >  1:1077334/1‑149 (MQ=255)
                                                                                                gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATg                                                   <  1:14481/149‑1 (MQ=255)
                                                                                                 aTCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAaataataatta                                                  >  2:758206/1‑147 (MQ=255)
                                                                                                 aTCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGc                                                  >  2:291954/1‑149 (MQ=255)
                                                                                                      tattatTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTcgc                                             >  2:1251864/1‑149 (MQ=255)
                                                                                                      tattatTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGAAAATAATAATGCGTcgc                                             >  2:942060/1‑149 (MQ=255)
                                                                                                       attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTcgcg                                            >  2:1115120/1‑149 (MQ=255)
                                                                                                       attattGATTATAAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGAGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGGcgcg                                            <  1:1408950/149‑1 (MQ=255)
                                                                                                            tGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATa                                       <  2:359240/149‑1 (MQ=255)
                                                                                                             gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaa                                      <  2:1080875/149‑1 (MQ=255)
                                                                                                              aTTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                     <  1:58033/149‑1 (MQ=255)
                                                                                                               ttATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAt                                    <  2:1440762/149‑1 (MQ=255)
                                                                                                               ttATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAt                                    <  1:448739/149‑1 (MQ=255)
                                                                                                                 ataaaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGa                                  <  1:758206/146‑1 (MQ=255)
                                                                                                                 aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGa                                  <  1:1514461/149‑1 (MQ=255)
                                                                                                                 aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGa                                  <  1:852822/149‑1 (MQ=255)
                                                                                                                  tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                 >  1:571075/1‑149 (MQ=255)
                                                                                                                  tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                 >  1:742390/1‑149 (MQ=255)
                                                                                                                  tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                 <  1:225994/149‑1 (MQ=255)
                                                                                                                    aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                               >  1:1562113/1‑149 (MQ=255)
                                                                                                                     aaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                              >  2:1237221/1‑149 (MQ=255)
                                                                                                                      aaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAtt                             <  2:830040/149‑1 (MQ=255)
                                                                                                                       aTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAttt                            >  2:147167/1‑149 (MQ=255)
                                                                                                                        ttACTCTAATAAAAAGCACATTGTATTCAGTCAACACTGTGCTTTTTTTAAATTCATAAGAATTTTGTATTACTGTCTTTTAACACTTAATTAATCATTTCCTTATTTATCCGCAAATAATAATTCGTCGCGCAATAAATGATCATTTc                           <  2:248285/149‑1 (MQ=255)
                                                                                                                         tAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTcc                          <  1:725536/149‑1 (MQ=255)
                                                                                                                          aaTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCa                         <  2:1490197/149‑1 (MQ=255)
                                                                                                                           aTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCac                        <  2:69267/149‑1 (MQ=255)
                                                                                                                              taataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGAGCAATAAATGATCATTTCCacgac                     >  1:499393/1‑149 (MQ=255)
                                                                                                                                 taaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgac                  >  2:911719/1‑149 (MQ=255)
                                                                                                                                 aaaaaaGCACATTGTTTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGTAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgac                  >  2:585817/2‑149 (MQ=38)
                                                                                                                                   aaaaGCACATTGTTTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACTAcgc                <  2:1353331/149‑1 (MQ=39)
                                                                                                                                    aaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgcg               >  2:1380482/1‑149 (MQ=255)
                                                                                                                                        cacaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTaaa           <  1:529227/149‑1 (MQ=255)
                                                                                                                                        cacaTTGTATTCATTAAACAATGTGCTTTTTTTAAAATCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTaaa           <  2:185945/149‑1 (MQ=255)
                                                                                                                                          caTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGc         <  2:1078251/149‑1 (MQ=255)
                                                                                                                                          caTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGc         <  2:851225/149‑1 (MQ=255)
                                                                                                                                           aTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATTCGTCGCGCAATAAATGATCATTTCAACGAAGACGCGTAAAGca        >  1:329752/1‑148 (MQ=255)
                                                                                                                                           aTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGcc        >  2:1567489/1‑149 (MQ=255)
                                                                                                                                           aTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGcc        >  2:107688/1‑149 (MQ=255)
                                                                                                                                           aTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATAATTTTAAATGAAATAATAATGATTCTAGAAATAAATTATAATTTCAACTACGACGCTTAAAGca        >  1:103775/1‑148 (MQ=255)
                                                                                                                                              ttATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAAt     <  2:1191949/148‑1 (MQ=255)
                                                                                                                                               tATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATAATTTCAACGACTACGCGTAAAGACAATg    >  1:225381/1‑149 (MQ=255)
                                                                                                                                                 ttCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCg  >  2:156088/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCG  >  NC_000913/3757796‑3758088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: