Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,095,151 T→C 7.5% A18A (GCT→GCC)  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,1510TC 7.5% 212.1 / 5.5 66A18A (GCT→GCC) yggSUPF0001 family protein, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/3);  total (33/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 = 1.00e+00

GCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCG  >  NC_000913/3095013‑3095297
                                                                                                                                          |                                                                                                                                                  
gcgcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCACCGGCTGCAACGCGtt                                                                                                                                          <  4:511059/149‑1 (MQ=255)
    tGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGcag                                                                                                                                      <  2:655924/149‑3 (MQ=255)
             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                                                                                                        >  2:392468/1‑149 (MQ=255)
                                  ttACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGc                                                                                                        >  4:86197/1‑149 (MQ=255)
                                              ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTTCGCTGCTTGCAGTCAGTaaaac                                                                                            >  3:423618/1‑149 (MQ=255)
                                              ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTaaaac                                                                                            >  1:531466/1‑149 (MQ=255)
                                              ccACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTaaaac                                                                                            >  4:189517/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                                                                             >  4:268506/1‑130 (MQ=255)
                                                                                ggAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAgcgc                                                                             <  3:268506/130‑1 (MQ=255)
                                                                                    aaTGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCg                                                      >  4:153470/1‑149 (MQ=255)
                                                                                               ttGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCaa                                           >  2:500418/1‑149 (MQ=255)
                                                                                               ttGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCaa                                           >  2:229633/1‑149 (MQ=255)
                                                                                                 gcgcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAAtt                                         <  1:272669/149‑1 (MQ=255)
                                                                                                 gcgcATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAAtt                                         <  2:125083/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)
                                                                                                                                          |                                                                                                                                                  
GCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAAATTACGCTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCG  >  NC_000913/3095013‑3095297

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