Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,548,243 T→C V940A (GTC→GCC)  fdnG → formate dehydrogenase‑N, nitrate‑inducible, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,548,2430TC98.9% 280.1 / ‑4.6 91V940A (GTC→GCC) fdnGformate dehydrogenase‑N, nitrate‑inducible, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (49/41);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

TCAACGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGTCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACGG  >  NC_000913/1548153‑1548332
                                                                                          |                                                                                         
tCAACGCAATTGCTCAGCCGGAACAGTTTGTGGAAATTAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGcct                                                                                         >  2:838567/1‑93 (MQ=255)
tCAACGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGcct                                                                                         >  1:2434280/1‑93 (MQ=255)
tCAACGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGcct                                                                                         >  1:609330/1‑93 (MQ=255)
   aCGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCa                                                                                      <  1:280290/93‑1 (MQ=255)
   aCGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCa                                                                                      >  1:3062590/1‑93 (MQ=255)
   aCGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCa                                                                                      >  2:1059579/1‑93 (MQ=255)
       aaTTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCaa                                                                                  <  1:4029270/93‑1 (MQ=255)
        aTTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAg                                                                                 >  1:443299/1‑93 (MQ=255)
             tCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTgg                                                                            >  1:3445964/1‑93 (MQ=255)
             tCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTgg                                                                            >  1:1276963/1‑93 (MQ=255)
                gCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCtt                                                                         >  2:1485591/1‑93 (MQ=255)
                  cGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGGTTTa                                                                       >  1:2577572/1‑93 (MQ=255)
                  cGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCtttt                                                                       >  2:2820905/1‑92 (MQ=255)
                  cGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTa                                                                       >  2:873107/1‑93 (MQ=255)
                  cGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTa                                                                       <  1:538017/93‑1 (MQ=255)
                   ggAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTAt                                                                      <  1:3080076/93‑1 (MQ=255)
                   ggAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTAt                                                                      <  1:3698409/93‑1 (MQ=255)
                   ggAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTAt                                                                      >  1:965008/1‑93 (MQ=255)
                   ggAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTAt                                                                      <  2:1282985/93‑1 (MQ=255)
                    gAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATc                                                                     >  1:874282/1‑93 (MQ=255)
                    gAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATc                                                                     <  2:357877/93‑1 (MQ=255)
                     aaCAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATcc                                                                    <  1:4354717/93‑1 (MQ=255)
                     aaCAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATcc                                                                    <  1:4151563/93‑1 (MQ=255)
                        aGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcg                                                                 <  1:1440113/93‑1 (MQ=255)
                        aGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcg                                                                 >  2:4266648/1‑93 (MQ=255)
                         gTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcgc                                                                >  1:4593882/1‑93 (MQ=255)
                         gTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcgc                                                                <  2:4374762/93‑1 (MQ=255)
                         gTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcgc                                                                >  2:13505/1‑93 (MQ=255)
                         gTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcgc                                                                >  1:4224822/1‑93 (MQ=255)
                          tttGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCcgcgcg                                                               >  1:3202947/1‑93 (MQ=255)
                           ttGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCgg                                                              >  2:2456843/1‑93 (MQ=255)
                            tgtgGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGt                                                             <  2:2433641/93‑1 (MQ=255)
                             gtgGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTg                                                            <  2:969802/93‑1 (MQ=255)
                             gtgGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTg                                                            >  1:830417/1‑93 (MQ=255)
                              tgGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTgg                                                           <  2:2183808/93‑1 (MQ=255)
                               ggAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGc                                                          >  1:316941/1‑93 (MQ=255)
                                  aaTCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCtgt                                                       <  1:5697/93‑1 (MQ=255)
                                   aTCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCtgtg                                                      >  2:990117/1‑93 (MQ=255)
                                   aTCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCtgtg                                                      >  2:4515784/1‑93 (MQ=255)
                                    tCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTgg                                                     >  2:3687142/1‑93 (MQ=255)
                                    tCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTgg                                                     >  2:4088877/1‑93 (MQ=255)
                                      aGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGCGGCTGTGGTa                                                   >  1:3171512/1‑93 (MQ=255)
                                       gCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTaa                                                  <  1:3701041/93‑1 (MQ=255)
                                        cGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAAc                                                 <  2:766180/93‑1 (MQ=255)
                                         gAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAAcg                                                >  1:2336263/1‑93 (MQ=255)
                                          aaaCGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAAcgc                                               <  1:231474/93‑1 (MQ=255)
                                           aaCGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGTCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAAcgcg                                              >  1:3098478/1‑93 (MQ=255)
                                             cGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGcgtc                                            <  2:2712847/93‑1 (MQ=255)
                                             cGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGcgtc                                            <  1:1517582/93‑1 (MQ=255)
                                              gCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGcgtcg                                           >  1:4216642/1‑93 (MQ=255)
                                                tGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCgtcgtc                                         <  1:4358914/93‑1 (MQ=255)
                                                tGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCgtcgtc                                         >  1:3574313/1‑93 (MQ=255)
                                                tGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCgtcgtc                                         <  1:4506378/93‑1 (MQ=255)
                                                 ggcggcggcGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCt                                        >  1:1868940/1‑93 (MQ=255)
                                                 ggcggcggcGAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCt                                        <  1:1271350/93‑1 (MQ=255)
                                                 gacggcggcgAAAGGTGTTAATGATGGCGGTCGGGTCACTGCCTCCAGAAAGCGTGGTTTTATCCGGGCGGTGGCTGTGGTAACGCGTCGTCt                                        <  2:2126640/91‑1 (MQ=255)
                                                     gcggcgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGTCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGaaa                                    <  1:537159/93‑1 (MQ=255)
                                                       ggcgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAAcc                                  <  1:4369347/93‑1 (MQ=255)
                                                       ggcgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAAcc                                  <  2:3478673/93‑1 (MQ=255)
                                                        gcgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCg                                 >  2:1993953/1‑93 (MQ=255)
                                                        gcgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCg                                 >  2:2759873/1‑93 (MQ=255)
                                                         cgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGc                                <  1:4241649/93‑1 (MQ=255)
                                                         cgAAAGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGc                                <  2:547427/93‑1 (MQ=255)
                                                           aaaGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTg                              >  2:1572264/1‑93 (MQ=255)
                                                            aaGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGa                             >  2:2895120/1‑93 (MQ=255)
                                                             aGGCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGaa                            >  2:2972196/1‑93 (MQ=255)
                                                              agCATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGACTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAAt                           <  2:864550/92‑1 (MQ=255)
                                                                cATTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGt                         >  2:1329456/1‑93 (MQ=255)
                                                                 aTTAATAATGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTa                        >  1:4253135/1‑93 (MQ=255)
                                                                   taataaTGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTaaa                      >  1:134766/1‑93 (MQ=255)
                                                                    aataatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAAt                     <  2:3063753/93‑1 (MQ=255)
                                                                     ataatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATg                    >  2:956175/1‑93 (MQ=255)
                                                                     ataatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATg                    <  1:715304/93‑1 (MQ=255)
                                                                      taatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATgg                   >  1:1260525/1‑93 (MQ=255)
                                                                      taatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATgg                   <  1:3641645/93‑1 (MQ=255)
                                                                      taatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATgg                   >  1:3965854/1‑93 (MQ=255)
                                                                      taatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATgg                   <  2:4306091/93‑1 (MQ=255)
                                                                       aatGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGt                  <  2:3692054/93‑1 (MQ=255)
                                                                         tGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTca                <  2:1286894/93‑1 (MQ=255)
                                                                         tGGCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTca                <  2:724250/93‑1 (MQ=255)
                                                                          ggCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTcag               <  1:93604/93‑1 (MQ=255)
                                                                           gCGATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTcagc              <  2:1987733/93‑1 (MQ=255)
                                                                             gATCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGGGCGGTGGCTGTGGTAACGCGTCTTCTGAAGCCGCTGAATGTAAATGGTcagcag            <  2:3434384/93‑1 (MQ=255)
                                                                               tCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGt          >  2:3625254/1‑93 (MQ=255)
                                                                               tCGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGt          <  1:2589727/93‑1 (MQ=255)
                                                                                cGTGTCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGtt         <  1:1714958/93‑1 (MQ=255)
                                                                                 gtgtCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTg        >  1:1404318/1‑93 (MQ=255)
                                                                                 gtgtCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTg        >  1:1075234/1‑93 (MQ=255)
                                                                                 gtgtCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTg        >  2:4244545/1‑93 (MQ=255)
                                                                                 gtgtCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTg        >  2:2029474/1‑93 (MQ=255)
                                                                                  tgtCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGa       >  1:1411978/1‑93 (MQ=255)
                                                                                    tCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGaaa     >  2:2614234/1‑93 (MQ=255)
                                                                                    tCACTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGaaa     <  2:1569507/93‑1 (MQ=255)
                                                                                      aCTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACg   >  2:2642224/1‑93 (MQ=255)
                                                                                       cTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACgg  <  1:3560777/93‑1 (MQ=255)
                                                                                       cTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACgg  >  1:1922889/1‑93 (MQ=255)
                                                                                       cTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACgg  >  1:391510/1‑93 (MQ=255)
                                                                                       cTGCCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACgg  >  1:4647302/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCAACGCAATTGCTCAGCCGGAACAGTTTGTGGAAATCAGCGAAACGCTGGCGGCGGCGAAAGGCATTAATAATGGCGATCGTGTCACTGTCTCCAGCAAGCGTGGCTTTATCCGCGCGGTGGCTGTGGTAACGCGTCGTCTGAAACCGCTGAATGTAAATGGTCAGCAGGTTGAAACGG  >  NC_000913/1548153‑1548332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: