Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 707,926 G→T 6.6% intergenic (+169/‑408) glnS → / → chiP glutamyl‑tRNA synthetase/chitoporin, uptake of chitosugars

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913707,9260GT 6.6% 209.7 / 4.4 61intergenic (+169/‑408)glnS/chiPglutamyl‑tRNA synthetase/chitoporin, uptake of chitosugars
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base T (2/2);  total (29/32)
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.98e-01

GGGCGTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTAACCG  >  NC_000913/707782‑708069
                                                                                                                                                |                                                                                                                                               
gggCGTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACACTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                             >  3:390688/1‑149 (MQ=255)
  gCGTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                           <  4:163679/149‑1 (MQ=25)
    gTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                         >  1:2347/1‑149 (MQ=255)
    gTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                         >  1:578246/1‑149 (MQ=255)
     tttACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                        >  4:129208/1‑149 (MQ=255)
           tCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAAc                                                                                                                                  >  3:32049/1‑149 (MQ=255)
             gCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAg                                                                                                                                >  3:167581/1‑149 (MQ=255)
               aTCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCa                                                                                                                              >  1:214177/1‑149 (MQ=255)
                 ccGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCaaa                                                                                                                            <  3:227166/149‑1 (MQ=255)
                    gCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACaa                                                                                                                         <  2:696927/149‑1 (MQ=255)
                    gCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACaa                                                                                                                         <  1:703933/149‑1 (MQ=255)
                      aaGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATc                                                                                                                       <  1:537559/149‑1 (MQ=255)
                      aaGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATc                                                                                                                       <  2:213416/149‑1 (MQ=255)
                      aaGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATc                                                                                                                       <  4:727437/149‑1 (MQ=255)
                        ggAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCa                                                                                                                     <  2:721903/149‑1 (MQ=255)
                          aaCAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCaaa                                                                                                                   >  1:265608/1‑149 (MQ=255)
                                 aaaCAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAAcgccgc                                                                                                            >  2:116003/1‑149 (MQ=255)
                                  aaCAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCcgcg                                                                                                           >  4:641698/1‑149 (MQ=255)
                                   aCAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAAcgccgctg                                                                                                          >  4:741566/1‑147 (MQ=255)
                                        ttCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAgc                                                                                                     <  1:170744/149‑1 (MQ=255)
                                         tCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAgcg                                                                                                    <  3:80593/149‑1 (MQ=255)
                                           aaaCGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAgcggc                                                                                                  >  3:224500/1‑149 (MQ=255)
                                           aaaCGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAgcggc                                                                                                  >  1:597534/1‑149 (MQ=255)
                                                cTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGtttt                                                                                             <  3:741566/149‑1 (MQ=18)
                                                 tAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGttttg                                                                                            <  2:2347/149‑2 (MQ=21)
                                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGtttttt                                                                                           <  4:457371/149‑1 (MQ=255)
                                                   aaTTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTc                                                                                          >  4:240354/1‑149 (MQ=255)
                                                     ttGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTg                                                                                        >  2:351402/1‑149 (MQ=255)
                                                     ttGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTg                                                                                        <  1:209202/149‑1 (MQ=255)
                                                        ccTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCtt                                                                                     <  3:66095/149‑1 (MQ=255)
                                                                     gCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGaatt                                                                        <  3:417457/149‑1 (MQ=255)
                                                                          tCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTc                                                                   <  2:718899/149‑1 (MQ=255)
                                                                                         ctctGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTAt                                                    <  2:578246/149‑1 (MQ=255)
                                                                                         ctctGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTAt                                                    <  3:611733/149‑1 (MQ=255)
                                                                                          tctGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATc                                                   <  3:395900/149‑1 (MQ=255)
                                                                                          tctGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATc                                                   <  3:641698/149‑1 (MQ=255)
                                                                                               aaTATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCataa                                              >  3:81763/1‑149 (MQ=255)
                                                                                                atatatTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCataaa                                             >  1:531454/1‑149 (MQ=255)
                                                                                                     ttGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCataaaataaa                                        >  3:535084/1‑149 (MQ=255)
                                                                                                     ttGAGTTTGCGCGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCataaaataaa                                        >  1:403721/1‑149 (MQ=255)
                                                                                                      tGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAaaataaat                                       <  2:206795/149‑1 (MQ=255)
                                                                                                         gTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTa                                    >  4:101751/1‑149 (MQ=255)
                                                                                                            tGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTaaaa                                 <  2:597534/149‑1 (MQ=255)
                                                                                                               gTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATaa                              <  4:615656/149‑1 (MQ=255)
                                                                                                                tGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAAc                             <  2:695172/149‑1 (MQ=255)
                                                                                                                tGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAAc                             <  4:535084/149‑1 (MQ=255)
                                                                                                                   ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACaaa                          <  3:240354/149‑1 (MQ=255)
                                                                                                                           gcgggATAAGGCGTTCACGCCTCACCCGGCAAGAAACAGCAAACAATCCAAAACGACGCGTTCAGCGGTGTTTTTTCTGCTTTTCGTCGAGAATTAATTCCGCTTAGCAATTGATCCATGAAAGAAATTTAAAATAACAAAACATAAtt                  >  4:470896/4‑149 (MQ=255)
                                                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAAt                   >  2:24118/1‑148 (MQ=255)
                                                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAAt                   <  1:24118/148‑1 (MQ=255)
                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAGCGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATACATAAAATAAATTTAAAATAACAAAACATAAtt                  >  4:501178/1‑149 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATta                 >  4:295425/1‑149 (MQ=255)
                                                                                                                              ggATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATtaaa               <  4:81763/149‑1 (MQ=255)
                                                                                                                                  aaGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATtaaataaa           <  1:116003/149‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCaaaa                                                                                                                  <  1:383070/44‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCaaaa                                                                                                                  >  2:383070/1‑44 (MQ=255)
                                                                                                                                       gTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTa      >  2:768661/1‑149 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTaa     <  3:259740/149‑1 (MQ=255)
                                                                                                                                          cACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTAAcc   >  4:375645/1‑149 (MQ=255)
                                                                                                                                           aCGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTAACCg  >  3:641016/1‑149 (MQ=255)
                                                                                                                                           aCGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAAGTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTAACCg  >  3:515823/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGGCGTTTACGTCGCATCCGGCAAGGAACAGACAAACAGTTTCAAACGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTTTTCTGCTTTTCTTCGCGAATTAATTCCGCTTCGCAATTTATCCATAAAATAAATTTAAAATAACAAAACATAATTAAATAAAATGTAACCG  >  NC_000913/707782‑708069

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: ‑