Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,095,164 G→A 6.0% G23S (GGC→AGC)  yggS → UPF0001 family protein, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,095,1640GA 6.0% 221.0 / 3.5 67G23S (GGC→AGC) yggSUPF0001 family protein, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAAC  >  NC_000913/3095026‑3095311
                                                                                                                                          |                                                                                                                                                   
gtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTcc                                                                                                                                           <  4:647490/149‑1 (MQ=255)
 tAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCa                                                                                                                                          <  4:444171/149‑1 (MQ=255)
   gATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCaga                                                                                                                                        <  2:443100/149‑1 (MQ=255)
     tCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGCAGcc                                                                                                                                                  <  2:71902/137‑1 (MQ=255)
     tCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGCAGcc                                                                                                                                                  >  1:71902/1‑137 (MQ=255)
        gagaCGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGaaa                                                                                                                                   <  2:460026/149‑1 (MQ=255)
           aCGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTa                                                                                                                                <  1:144990/149‑1 (MQ=255)
             gTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACg                                                                                                                              >  3:353951/1‑149 (MQ=255)
             gTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACg                                                                                                                              >  3:599128/1‑149 (MQ=255)
             gTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACg                                                                                                                              >  2:272669/1‑149 (MQ=255)
                 tGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgctgc                                                                                                                          <  4:353951/149‑1 (MQ=255)
                  gCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgctgct                                                                                                                         >  3:666508/1‑149 (MQ=255)
                    tttACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgctgc                                                                                                                          >  1:458100/1‑146 (MQ=255)
                    tttACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgctgc                                                                                                                          <  2:458100/146‑1 (MQ=255)
                     ttACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGc                                                                                                                      >  4:86197/1‑149 (MQ=255)
                     ttACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGc                                                                                                                      >  2:392468/1‑149 (MQ=255)
                                 ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTTCGCTGCTTGCAGTCAGTaaaac                                                                                                          >  3:423618/1‑149 (MQ=255)
                                 ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTaaaac                                                                                                          >  4:189517/1‑149 (MQ=255)
                                 ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTaaaac                                                                                                          >  1:531466/1‑149 (MQ=255)
                                     aaTTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTaaaacaaaa                                                                                                      <  1:585692/149‑1 (MQ=255)
                                         tcttcCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTg                                                                                                  >  4:137226/1‑149 (MQ=255)
                                         tcttcCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTg                                                                                                  >  3:274160/1‑149 (MQ=255)
                                                  ttCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCa                                                                                         >  3:531848/1‑149 (MQ=255)
                                                    cATACTAAGATCCTCGTAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATc                                                                                       <  3:137226/149‑1 (MQ=255)
                                                      tACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgc                                                                                                                             >  4:41907/1‑109 (MQ=255)
                                                      tACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACgc                                                                                                                             <  3:41907/109‑1 (MQ=255)
                                                      tACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGc                                                                                     >  4:576594/1‑149 (MQ=255)
                                                       aCTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGCAGCCgctg                                                                                                                                              <  3:309029/91‑4 (MQ=39)
                                                       aCTAAGATACTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGCAGCCgctg                                                                                                                                              >  4:309029/1‑88 (MQ=39)
                                                         tacgATCCTCGGCAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTa                                                                                                                                >  2:450238/4‑103 (MQ=255)
                                                         tAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTa                                                                                                                                <  1:450238/103‑1 (MQ=255)
                                                         tAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGa                                                                                  >  3:12829/1‑149 (MQ=255)
                                                          aaGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGATGGCGCTCGCAggt                                                                                 <  4:666508/149‑3 (MQ=255)
                                                                 tCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTg                                                                          <  3:584710/149‑1 (MQ=255)
                                                                   ggAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAgcgc                                                                                           <  3:268506/130‑1 (MQ=255)
                                                                   ggAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAgcgc                                                                                           >  4:268506/1‑130 (MQ=255)
                                                                       aaTGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCg                                                                    >  4:153470/1‑149 (MQ=255)
                                                                                  ttGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCaa                                                         >  2:229633/1‑149 (MQ=255)
                                                                                  ttGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCaa                                                         >  2:500418/1‑149 (MQ=255)
                                                                                    gcgcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAAtt                                                       <  2:125083/149‑1 (MQ=255)
                                                                                    gcgcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAAtt                                                       <  1:272669/149‑1 (MQ=255)
                                                                                      gcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGaa                                                                                 <  1:240144/121‑1 (MQ=255)
                                                                                      gcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGaa                                                                                 >  2:240144/1‑121 (MQ=255)
                                                                                        aTAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGt                                                   >  2:82066/1‑149 (MQ=255)
                                                                                        aTAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGACAGCGTCAATTTGGt                                                   >  4:372852/1‑149 (MQ=255)
                                                                                             cTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGaaaa                                              <  3:153470/149‑1 (MQ=255)
                                                                                                     ggTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTc                                      >  2:558485/1‑149 (MQ=255)
                                                                                                       tCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTATAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAg                                    <  4:599128/149‑1 (MQ=255)
                                                                                                           ggACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAg                                <  4:274160/149‑1 (MQ=255)
                                                                                                            gACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgg                               <  3:576594/149‑1 (MQ=255)
                                                                                                             aCAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGTTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCATCAATTTGGTGAAAACTACGTTCAGGAAggg                              <  4:423618/149‑1 (MQ=255)
                                                                                                             aCAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAggg                              <  4:12829/149‑1 (MQ=255)
                                                                                                             aCAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAggg                              >  2:250259/1‑149 (MQ=255)
                                                                                                             aCAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGACAGCGTCAATTTGGTGAAAACTACGTTCAGGAAggg                              <  3:372852/149‑1 (MQ=255)
                                                                                                              cAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggg                             >  3:89914/1‑149 (MQ=255)
                                                                                                               aaaaTCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGt                            >  4:31042/1‑149 (MQ=255)
                                                                                                                aaaTCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTa                           >  4:112048/1‑149 (MQ=255)
                                                                                                                 aaTCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATc                                                                                       >  2:569933/1‑88 (MQ=255)
                                                                                                                 aaTCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATc                                                                                       <  1:569933/88‑1 (MQ=255)
                                                                                                                     tcAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTgg                                                    >  4:334619/1‑119 (MQ=255)
                                                                                                                     tcAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTgg                                                    <  3:334619/119‑1 (MQ=255)
                                                                                                                       aGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATaaaa                    <  3:189517/149‑1 (MQ=255)
                                                                                                                          ggCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTc                 <  3:86197/149‑1 (MQ=255)
                                                                                                                           gCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCg                <  1:82066/149‑1 (MQ=255)
                                                                                                                                aaCGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACt           >  2:533408/1‑149 (MQ=255)
                                                                                                                                        gCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGaa   <  2:531466/149‑1 (MQ=255)
                                                                                                                                         cGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAAc  <  4:568615/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAAC  >  NC_000913/3095026‑3095311

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