Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,377 G→A 6.2% intergenic (+386/+21) 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,3770GA 6.2% 263.3 / 5.3 80intergenic (+386/+21)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (46/29);  new base A (3/2);  total (49/31)
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.59e-01

ACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATC  >  NC_000913/1552232‑1552524
                                                                                                                                                 |                                                                                                                                                   
aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCtt                                                                                                                                                  <  4:587607/149‑1 (MQ=35)
  gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                                                                                                                                                >  4:530975/1‑149 (MQ=38)
          tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                                                                                                                                        <  4:442429/149‑1 (MQ=255)
                     gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                                                             <  4:865382/149‑1 (MQ=255)
                                aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAg                                                                                                                  >  3:739066/1‑149 (MQ=255)
                                aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtt                                                                                                                                    >  3:198803/1‑131 (MQ=255)
                                aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtt                                                                                                                                    <  4:198803/131‑1 (MQ=255)
                                 aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa                                                                                                                 <  3:530975/149‑1 (MQ=255)
                                  cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGAc                                                                                                                <  4:602903/149‑1 (MQ=255)
                                      ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg                                                                                                            >  1:841892/1‑149 (MQ=255)
                                         tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTaa                                                                                                         >  1:1110931/1‑149 (MQ=255)
                                                  aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAt                                                                                                                                     >  3:103335/1‑112 (MQ=255)
                                                  aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAt                                                                                                                                     <  4:103335/112‑1 (MQ=255)
                                                   aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGAcaca                                                                                               >  1:1128407/1‑149 (MQ=255)
                                                      cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATcc                                                                                            >  1:57375/1‑149 (MQ=255)
                                                         ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACt                                                                                         >  3:202440/1‑149 (MQ=255)
                                                           cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGt                                                                                       >  2:316751/1‑149 (MQ=255)
                                                            aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtt                                                                                      >  4:833603/1‑149 (MQ=255)
                                                             aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                             >  4:807641/1‑125 (MQ=255)
                                                             aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                             <  3:807641/125‑1 (MQ=255)
                                                                 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 <  3:833603/149‑1 (MQ=255)
                                                                 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 >  2:610511/1‑149 (MQ=255)
                                                                 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 >  3:886629/1‑149 (MQ=255)
                                                                 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 >  1:502168/1‑149 (MQ=255)
                                                                 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 >  4:923241/1‑149 (MQ=255)
                                                                  aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                <  2:1110931/149‑1 (MQ=255)
                                                                  aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                >  3:248696/1‑149 (MQ=255)
                                                                            aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCa                                                                      >  3:1068791/1‑149 (MQ=255)
                                                                               gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGAc                                                                   >  3:776463/1‑149 (MQ=255)
                                                                                   gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 >  2:510171/1‑131 (MQ=255)
                                                                                   gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                 <  1:510171/131‑1 (MQ=255)
                                                                                   gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCa                                                               <  2:57375/149‑1 (MQ=255)
                                                                                           cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCa                                                                      <  3:280566/134‑1 (MQ=255)
                                                                                           cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCa                                                                      >  4:280566/1‑134 (MQ=255)
                                                                                           cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                       >  1:49688/1‑149 (MQ=255)
                                                                                           cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                       >  4:667435/1‑149 (MQ=255)
                                                                                             gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGc                                                     >  2:532821/1‑149 (MQ=255)
                                                                                                aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                  <  4:1117636/149‑1 (MQ=255)
                                                                                                  gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCttt                                                                                                                                                 >  4:356938/1‑52 (MQ=255)
                                                                                                  gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCttt                                                                                                                                                 <  3:356938/52‑1 (MQ=255)
                                                                                                     tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAAGCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAg                                                   >  2:822400/1‑143 (MQ=255)
                                                                                                     tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAACTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAg                                                   <  1:822400/143‑1 (MQ=255)
                                                                                                      aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGc                                            >  3:117743/1‑149 (MQ=255)
                                                                                                           aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAAt                                                                                                        <  3:1023251/84‑1 (MQ=255)
                                                                                                           aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAAt                                                                                                        >  4:1023251/1‑84 (MQ=255)
                                                                                                           aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAt                                       >  3:54510/1‑149 (MQ=255)
                                                                                                             cAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAttt                                     <  1:532821/149‑1 (MQ=255)
                                                                                                              aGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTg                                    >  3:1065471/1‑149 (MQ=255)
                                                                                                               ggCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                   >  4:595921/1‑149 (MQ=255)
                                                                                                                 cAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcgg                                 >  3:914693/1‑149 (MQ=255)
                                                                                                                      ttGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGCAACCAGACGGCGTAATCGCGACACTTCCGCAGTATTATCTGCCTCTGCCAGACTACTGGCGTTTTGCAGACTCAGCCACTTTTGAGGCGTAc                            <  4:739066/149‑1 (MQ=255)
                                                                                                                        gCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTg                          >  1:659086/1‑149 (MQ=255)
                                                                                                                          gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                          >  2:537593/1‑83 (MQ=255)
                                                                                                                          gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                          <  1:537593/83‑1 (MQ=255)
                                                                                                                          gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGcc                        >  3:882311/1‑149 (MQ=255)
                                                                                                                          gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGcc                        >  1:803215/1‑149 (MQ=255)
                                                                                                                          gTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGcc                        >  1:606753/1‑149 (MQ=255)
                                                                                                                           tttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcg                                  <  3:852678/138‑1 (MQ=255)
                                                                                                                           tttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcg                                  >  4:852678/1‑138 (MQ=255)
                                                                                                                           tttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCg                       >  1:787962/1‑149 (MQ=255)
                                                                                                                            ttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGa                      <  2:803215/149‑1 (MQ=255)
                                                                                                                              gTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATc                    <  2:787962/149‑1 (MQ=255)
                                                                                                                               tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtttttt                                                                                  >  2:791456/1‑86 (MQ=255)
                                                                                                                               tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtttttt                                                                                  <  1:791456/86‑1 (MQ=255)
                                                                                                                               tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                         <  3:124893/111‑1 (MQ=255)
                                                                                                                               tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                         >  4:124893/1‑111 (MQ=255)
                                                                                                                               tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATca                   >  1:476242/1‑149 (MQ=255)
                                                                                                                                 atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacc                 >  4:98316/1‑149 (MQ=255)
                                                                                                                                   caGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCTACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcaccac               >  4:782295/1‑149 (MQ=255)
                                                                                                                                      ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacca                >  3:77516/1‑145 (MQ=255)
                                                                                                                                      ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa            >  3:358834/1‑149 (MQ=255)
                                                                                                                                      ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCACCCACATTTGCTGCGAACTGCCGATcagca                <  4:77516/145‑1 (MQ=255)
                                                                                                                                        cTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa          <  3:923241/149‑1 (MQ=255)
                                                                                                                                        cTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaaa          <  1:316751/149‑1 (MQ=255)
                                                                                                                                         tAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGc                         >  4:115275/1‑133 (MQ=255)
                                                                                                                                         tAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGc                         <  3:115275/133‑1 (MQ=255)
                                                                                                                                         tAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATca                   >  4:703666/1‑139 (MQ=255)
                                                                                                                                         tAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATca                   <  3:703666/139‑1 (MQ=255)
                                                                                                                                         aaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAg         >  1:46619/2‑149 (MQ=255)
                                                                                                                                          aaaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCTTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt        <  3:782295/149‑1 (MQ=255)
                                                                                                                                               gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAt   >  2:1113656/1‑149 (MQ=255)
                                                                                                                                               gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAt   >  1:317239/1‑149 (MQ=255)
                                                                                                                                                cgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATc  <  4:882311/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATC  >  NC_000913/1552232‑1552524

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: