Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,388 T→C 7.3% intergenic (+397/+10) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,3880TC 7.3% 352.6 / 13.0 109intergenic (+397/+10)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/46);  new base C (4/4);  total (59/50)
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.66e-01

Only 100 of 111 total aligned reads displayed.

GCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGC  >  NC_000913/1552249‑1552536
                                                                                                                                           |                                                                                                                                                    
gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                                                                                                                                             >  1:1317396/1‑149 (MQ=255)
 ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                                                                                                                                            <  2:1317396/149‑1 (MQ=255)
  cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                                                                                                                                           >  1:416137/1‑149 (MQ=255)
    gataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                                                                         >  1:1072332/6‑149 (MQ=17)
         ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt                                                                                                                                    >  3:1068882/1‑149 (MQ=255)
           cgttcaCGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa                                                                                                                                  <  2:1072332/144‑1 (MQ=255)
           cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa                                                                                                                                  <  2:416137/149‑1 (MQ=255)
             tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAcc                                                                                                                                <  4:1068882/149‑1 (MQ=255)
               aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGcc                                                                                                                                                   <  1:1404815/128‑1 (MQ=255)
               aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGcc                                                                                                                                                   >  2:1404815/1‑128 (MQ=255)
                 cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGcc                                                                                                                                                   <  3:1167944/126‑1 (MQ=255)
                    cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                                                                                                                                            >  3:1259404/1‑130 (MQ=255)
                    cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt                                                                                                                                    >  4:809369/1‑138 (MQ=255)
                    cTTATCCGGAATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTTCCAtttttt                                                                                                                                            <  4:1259404/130‑1 (MQ=255)
                          cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATc                                                                                                                   <  4:1388119/149‑1 (MQ=255)
                          cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATc                                                                                                                   >  1:679129/1‑149 (MQ=255)
                             cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgcg                                                                                                                >  2:191512/1‑149 (MQ=255)
                              aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        >  2:622179/1‑124 (MQ=255)
                              aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        <  1:622179/124‑1 (MQ=255)
                              aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               <  2:679129/149‑1 (MQ=255)
                               tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               <  3:117379/148‑1 (MQ=255)
                               tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               >  4:117379/1‑148 (MQ=255)
                                  aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa                                                                                                                                  >  2:509086/1‑126 (MQ=255)
                                  aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa                                                                                                                                  <  1:509086/126‑1 (MQ=255)
                                       cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                      >  4:929462/1‑149 (MQ=255)
                                       cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                      <  3:929462/149‑1 (MQ=255)
                                       cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                      >  1:386783/1‑149 (MQ=255)
                                         tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgc                                                                                                                 >  2:1381202/1‑136 (MQ=255)
                                         tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgc                                                                                                                 <  1:1381202/136‑1 (MQ=255)
                                              tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGttttt                                                                                               >  1:15687/1‑149 (MQ=255)
                                                  aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                                                                                                                                                    <  4:600428/92‑1 (MQ=255)
                                                  aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                                                                                                                                                    >  3:600428/1‑92 (MQ=255)
                                                  aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAgcgc                                                                                                                                                    <  1:1006716/92‑1 (MQ=255)
                                                   tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                          <  2:15687/149‑1 (MQ=255)
                                                     gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                         >  4:811119/1‑116 (MQ=255)
                                                     gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                         <  3:811119/116‑1 (MQ=255)
                                                           aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCa                                                                                  >  3:146508/1‑149 (MQ=255)
                                                           aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttt                                                                                                                                               >  3:1101047/1‑88 (MQ=255)
                                                           aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttt                                                                                                                                               <  4:1101047/88‑1 (MQ=255)
                                                                aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGcc                                                                                                                                                   >  1:1230903/1‑79 (MQ=255)
                                                                aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGcc                                                                                                                                                   <  2:1230903/79‑1 (MQ=255)
                                                                 cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTcc                                                                            <  4:146508/149‑1 (MQ=255)
                                                                  gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCa                                                                           >  1:901251/1‑149 (MQ=255)
                                                                    gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                 >  2:21417/1‑93 (MQ=255)
                                                                    gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                 <  1:21417/93‑1 (MQ=255)
                                                                      aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                          >  3:694851/1‑130 (MQ=255)
                                                                      aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                          <  4:694851/130‑1 (MQ=255)
                                                                            gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGc                                                                 >  1:1118542/1‑149 (MQ=255)
                                                                             ttcagtgTAGCGCATCAGGTAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTCTATAAAGTGCCATTTTTTTATTGAGTAACCAGACGGCGTAATCGCGACACATCCTCTGTTTTTTCTGCCTCTGCCCGACTCCAGTCGTTTTGCa                                                                >  2:471920/7‑149 (MQ=38)
                                                                             aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               <  4:36562/102‑1 (MQ=255)
                                                                             aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               >  3:36562/1‑102 (MQ=255)
                                                                             aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                                                                    <  4:1334351/145‑1 (MQ=255)
                                                                             aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                                                                    >  3:1334351/1‑145 (MQ=255)
                                                                             aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa                                                                <  1:191512/149‑1 (MQ=255)
                                                                               aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                                                                    >  1:233596/1‑143 (MQ=255)
                                                                               aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGttt                                                                    <  2:233596/143‑1 (MQ=255)
                                                                                   gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGAc                                                                                                                            <  4:415537/83‑1 (MQ=255)
                                                                                   gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGAc                                                                                                                            >  3:415537/1‑83 (MQ=255)
                                                                                     aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                            <  2:1442518/113‑1 (MQ=255)
                                                                                     aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGc                                                        <  2:901251/149‑1 (MQ=255)
                                                                                          aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAt                                                   >  4:928419/1‑149 (MQ=255)
                                                                                              ggCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        <  4:220114/60‑1 (MQ=255)
                                                                                              ggCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        >  3:220114/1‑60 (MQ=255)
                                                                                               gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                                                                                                                                             <  1:950039/54‑1 (MQ=255)
                                                                                               gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                                                                                                                                             >  2:950039/1‑54 (MQ=255)
                                                                                               gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACGCATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcg                                              >  1:845788/1‑149 (MQ=255)
                                                                                                         gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGt                                                                                                   >  2:609222/1‑86 (MQ=255)
                                                                                                         gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGt                                                                                                   <  1:609222/86‑1 (MQ=255)
                                                                                                         gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTc                                                                             <  4:1067431/108‑1 (MQ=255)
                                                                                                              tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctg                                                                                     >  4:1028731/1‑95 (MQ=255)
                                                                                                              tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctg                                                                                     <  3:1028731/95‑1 (MQ=255)
                                                                                                                atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               >  3:613846/1‑67 (MQ=255)
                                                                                                                atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               <  4:613846/67‑1 (MQ=255)
                                                                                                                 tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAGCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAttt                                                 >  3:997230/1‑128 (MQ=255)
                                                                                                                 tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCa                                                                                                       <  1:1074001/74‑1 (MQ=255)
                                                                                                                     ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa                        >  4:796721/1‑149 (MQ=255)
                                                                                                                       cTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      <  2:1118542/149‑1 (MQ=255)
                                                                                                                            tGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                               <  2:1318022/119‑1 (MQ=255)
                                                                                                                            tGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                               >  1:1318022/1‑119 (MQ=255)
                                                                                                                             ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                              <  3:338068/103‑1 (MQ=255)
                                                                                                                             ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                              >  4:338068/1‑103 (MQ=255)
                                                                                                                              gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACGCATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAt               <  2:845788/149‑1 (MQ=255)
                                                                                                                                gcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                    >  1:183053/1‑142 (MQ=255)
                                                                                                                                gcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                    <  2:183053/142‑1 (MQ=255)
                                                                                                                                 cTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      <  1:196672/139‑1 (MQ=255)
                                                                                                                                 cTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      >  2:196672/1‑139 (MQ=255)
                                                                                                                                 cTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      <  1:997512/139‑1 (MQ=255)
                                                                                                                                 cTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      >  2:997512/1‑139 (MQ=255)
                                                                                                                                  tttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                    >  4:1103488/1‑140 (MQ=255)
                                                                                                                                  tttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                    <  3:1103488/140‑1 (MQ=255)
                                                                                                                                     ataAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGACAGACTCCAGGCGTTTTGCa                                                                <  1:1306947/93‑1 (MQ=255)
                                                                                                                                     ataAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGACAGACTCCAGGCGTTTTGCa                                                                >  2:1306947/1‑93 (MQ=255)
                                                                                                                                       aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                     >  4:567997/1‑86 (MQ=255)
                                                                                                                                       aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                     <  3:567997/86‑1 (MQ=255)
                                                                                                                                       aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      <  3:110241/133‑1 (MQ=255)
                                                                                                                                       aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa                      >  4:110241/1‑133 (MQ=255)
                                                                                                                                       aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTc      >  2:1214568/1‑149 (MQ=255)
                                                                                                                                        aaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCt     >  2:961193/1‑149 (MQ=255)
                                                                                                                                         aGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTg    >  1:230253/1‑149 (MQ=255)
                                                                                                                                           tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGc  >  3:1098541/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGC  >  NC_000913/1552249‑1552536

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: