Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,109,196 Δ1 bp 5.9% coding (363/990 nt) sbp → sulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,109,1960T. 5.9% 228.1 / 6.5 68coding (363/990 nt)sbpsulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base . (2/2);  total (34/34)
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.55e-01

AACGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATCACAACAACAA  >  NC_000913/4109050‑4109339
                                                                                                                                                  |                                                                                                                                               
aaCGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGTCGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTaa                                                                                                                                               <  1:7896/149‑1 (MQ=255)
aaCGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGt                                                                                                                                                 >  3:213587/1‑147 (MQ=255)
aaCGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGt                                                                                                                                                 >  3:74337/1‑147 (MQ=255)
aaCGGTATTGAAGCTGATGTTGTCACGCTGGCTATGGCCTATGACGTGTACGCAATTGCGGCACGCGTGCGGATTGATAAACAGTGGATCAAACGTCTGCCGGATAACACCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGt                                                                                                                                                 >  6:279874/1‑147 (MQ=255)
 aCGGTATTGAAGCTGATGTTGTCACGCTGGCTCTTGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGAACACTTCCACCATTGTTTTCCTGGTACGTAAg                                                                                                                                              >  1:354441/1‑149 (MQ=255)
  cGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAgg                                                                                                                                             >  3:263997/1‑149 (MQ=255)
       ttGAAGCGGTTGTTGTGACCGCGGCTCTGGCCTATTACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAATATTGGATCAAACGTCGTCCGGATAAATCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGaaa                                                                                                                                         <  6:62352/148‑1 (MQ=255)
        tGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATc                                                                                                                                       >  1:334537/1‑149 (MQ=255)
          aaGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATcc                                                                                                                                      <  4:42082/148‑1 (MQ=255)
           aGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAGTCCGa                                                                                                                                    <  4:263997/149‑1 (MQ=255)
                  gttgtCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCCCCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAg                                                                                                                                  >  8:107670/1‑144 (MQ=255)
                          gCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACt                                                                                                                     >  8:22350/1‑149 (MQ=255)
                           cTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                                                                                    <  4:74337/149‑1 (MQ=255)
                            tGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                                                                                    >  7:360395/1‑148 (MQ=255)
                                   ggCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                                                                            >  1:217581/1‑149 (MQ=255)
                                    gCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                                                                            <  5:78083/148‑1 (MQ=255)
                                    gCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                                                                            <  8:71597/148‑1 (MQ=255)
                                    gCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCt                                                                                                           <  2:354441/149‑1 (MQ=255)
                                      cTATGACGTGGACGCAATTGCGGAACGCGGGCTGATGGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                                                                         <  4:278455/149‑1 (MQ=255)
                                      cTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAACGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                                                                         <  6:279186/149‑1 (MQ=255)
                                       tATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAt                                                                                                        >  6:55252/1‑149 (MQ=255)
                                            cGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAc                                                                                                   <  8:360395/149‑1 (MQ=255)
                                                     aaTTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                                                                          >  7:329235/1‑149 (MQ=255)
                                                       ttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATcc                                                                                                                           >  1:157736/1‑114 (MQ=255)
                                                       ttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATcc                                                                                                                           <  2:157736/114‑1 (MQ=255)
                                                        tGCGGAACGCGGGCGGATTGATAAAGAGTGGATCATACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAACGGGTGTTTCg                                                                                        <  8:105202/148‑1 (MQ=255)
                                                          cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGt                                                                                      <  7:22350/148‑1 (MQ=255)
                                                            gAACGCGGGCGGATTGATTAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTACGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGAt                                                                                   <  8:369590/149‑1 (MQ=255)
                                                                      ggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAAt                                                                         >  4:155227/1‑149 (MQ=255)
                                                                        aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc                                                                       >  5:342796/1‑149 (MQ=255)
                                                                        aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc                                                                       >  3:310808/1‑148 (MQ=255)
                                                                        aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc                                                                       >  3:305600/1‑148 (MQ=255)
                                                                            aTAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                                                                         >  7:363525/1‑111 (MQ=255)
                                                                            aTAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                                                                         <  8:363525/111‑1 (MQ=255)
                                                                              aaaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                                                                          >  2:308606/1‑124 (MQ=255)
                                                                              aaaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGttt                                                                                          <  1:308606/124‑1 (MQ=255)
                                                                               aaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAg                                                                 <  3:155227/148‑1 (MQ=255)
                                                                                agagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGct                                                               >  7:349792/1‑149 (MQ=255)
                                                                                  agTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctct                                                             <  8:329235/149‑1 (MQ=255)
                                                                                   gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg                                                            <  5:55252/149‑1 (MQ=255)
                                                                                   gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg                                                            <  2:334537/149‑1 (MQ=255)
                                                                                    tGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTTTTTCGGTGATCACGCCTAATCCGAAAAGctctca                                                           >  5:276123/1‑147 (MQ=255)
                                                                                      gATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGCACGATCt                                                                                                           >  2:366439/1‑99 (MQ=255)
                                                                                      gATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctct                                                             >  2:91946/1‑145 (MQ=255)
                                                                                      gATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctct                                                             <  1:91946/145‑1 (MQ=255)
                                                                                      gATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtggtg                                                         <  4:213587/149‑1 (MQ=255)
                                                                                          aaaCGTCTGCCGGATAACTCCGCACCGTACACTTCAACCATTGTTTACCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGAATTGAACGATCTGATTAAACCGCGTGTTTCGGTGATCCCGCCTAATCCGAAAAGCTCtggaggaaa                                                     >  2:307379/1‑146 (MQ=255)
                                                                                                 tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                                                                                >  3:72153/1‑115 (MQ=255)
                                                                                                 tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                                                                                <  4:72153/115‑1 (MQ=255)
                                                                                                      gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTAc                                        <  4:305600/149‑1 (MQ=255)
                                                                                                      gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTAc                                        <  4:310808/149‑1 (MQ=255)
                                                                                                           cTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTgg                                    >  7:277506/1‑149 (MQ=255)
                                                                                                            tCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTgg                                    >  6:244737/1‑148 (MQ=255)
                                                                                                                 aCCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg                                                            >  7:188025/1‑119 (MQ=255)
                                                                                                                 aCCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg                                                            <  8:188025/119‑1 (MQ=255)
                                                                                                                          ttCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTAcg                     >  5:220624/1‑149 (MQ=255)
                                                                                                                          ttCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTAcg                     >  8:382260/1‑149 (MQ=255)
                                                                                                                            ccaccaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                                                                                <  1:184435/88‑1 (MQ=255)
                                                                                                                            ccaccaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                                                                                >  2:184435/1‑88 (MQ=255)
                                                                                                                            ccaccaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTAcg                     <  7:107670/147‑1 (MQ=255)
                                                                                                                             caccaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACgcgc                  >  1:114634/1‑149 (MQ=255)
                                                                                                                              accaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCt                 >  8:388230/1‑149 (MQ=255)
                                                                                                                              accaTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCt                 <  2:217581/149‑1 (MQ=255)
                                                                                                                                caTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTgg                            >  5:315948/1‑136 (MQ=255)
                                                                                                                                caTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTgg                            <  6:315948/136‑1 (MQ=255)
                                                                                                                                caTTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTg                >  1:107078/1‑148 (MQ=255)
                                                                                                                                  ttGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCAt             <  6:342796/149‑1 (MQ=255)
                                                                                                                                    gTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATca           >  8:336902/1‑149 (MQ=255)
                                                                                                                                     ttttCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCGGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATcac          <  8:349792/149‑1 (MQ=255)
                                                                                                                                        tCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATCacaa        <  7:382260/148‑1 (MQ=255)
                                                                                                                                             gTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATCAcaacaacaa  >  4:308378/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AACGGTATTGAAGCTGATGTTGTCACGCTGGCTCTGGCCTATGACGTGGACGCAATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGGTGGCGCGCGCTGGAACTACCTGGCAGCCTGGGGCTACGCGCTGCATCACAACAACAA  >  NC_000913/4109050‑4109339

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