Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,539,359 G→T P1078T (CCA→ACA)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,539,3590GT97.2% 600.2 / 3.4 179P1078T (CCA→ACA) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (96/78);  total (100/79)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 179 total aligned reads displayed.

GAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCT  >  NC_000913/1539261‑1539454
                                                                                                  |                                                                                               
gagCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGc                                                                                                >  2:2294460/1‑100 (MQ=255)
  gCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCgg                                                                                              >  1:4603626/1‑100 (MQ=255)
   cGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGt                                                                                             <  1:505986/100‑1 (MQ=255)
   cGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGt                                                                                             <  1:5188761/100‑1 (MQ=255)
     tCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTAt                                                                                           >  2:5222992/1‑100 (MQ=255)
      cAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATc                                                                                          <  1:4340748/100‑1 (MQ=255)
       aGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCt                                                                                         <  1:4519734/100‑1 (MQ=255)
        gCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTg                                                                                        >  1:1236701/1‑100 (MQ=255)
        gCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTg                                                                                        >  2:5301645/1‑100 (MQ=255)
          aTTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGgc                                                                                      >  1:107575/1‑100 (MQ=255)
             aGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCa                                                                                   >  1:1535807/1‑100 (MQ=255)
             aGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCa                                                                                   >  1:463573/1‑100 (MQ=255)
             aGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCa                                                                                   <  1:300338/100‑1 (MQ=255)
              gCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAt                                                                                  >  1:476040/1‑100 (MQ=255)
              gCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAt                                                                                  >  2:2093531/1‑100 (MQ=255)
               cAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAtc                                                                                 <  2:3981748/100‑1 (MQ=255)
               cAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAtc                                                                                 <  2:2009811/100‑1 (MQ=255)
                aGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAtct                                                                                >  1:4242670/1‑100 (MQ=255)
                aGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAtct                                                                                >  1:2435453/1‑100 (MQ=255)
                 ggTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCAtctc                                                                               >  2:2797267/1‑100 (MQ=255)
                   ttttCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCAc                                                                             >  1:881912/1‑100 (MQ=255)
                   ttttCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCAc                                                                             <  2:1513226/100‑1 (MQ=255)
                   tcttCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCAc                                                                             <  2:1200000/98‑1 (MQ=255)
                     ttCACTGTAGGTTGAGTGAATGCCCCAGTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTg                                                                           <  1:426664/100‑1 (MQ=255)
                        aCTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACg                                                                        >  2:461948/1‑100 (MQ=255)
                        aCTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACg                                                                        <  1:3398549/100‑1 (MQ=255)
                        aCTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACg                                                                        <  1:2530509/100‑1 (MQ=255)
                        aCTGTAGGTTGAGTGAATGACCCATTTCTGGTGCGGCGTCAGGAACTTAAGTGCCTTTTCCGGGAAGCCGTTCGTCGGTGTCTGGAGCATCTCCGTGACg                                                                        >  2:5032577/1‑100 (MQ=25)
                         cTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGc                                                                       >  1:614777/1‑100 (MQ=255)
                         cTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGc                                                                       >  2:384152/1‑100 (MQ=255)
                          tGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCt                                                                      <  2:1119590/100‑1 (MQ=255)
                          tGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCt                                                                      >  1:1276112/1‑100 (MQ=255)
                           gTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTa                                                                     >  1:1997965/1‑100 (MQ=255)
                            tAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTAc                                                                    <  1:3919753/100‑1 (MQ=255)
                            tAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTAc                                                                    <  1:2441480/100‑1 (MQ=255)
                            tAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTAc                                                                    >  1:2150311/1‑100 (MQ=255)
                             aGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACg                                                                   >  2:1951130/1‑100 (MQ=255)
                              ggTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACgg                                                                  >  2:4897870/1‑100 (MQ=255)
                              ggTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACgg                                                                  >  1:5462918/1‑100 (MQ=255)
                              ggTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACgg                                                                  >  1:4677147/1‑100 (MQ=255)
                              ggTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACgg                                                                  >  2:3441544/1‑100 (MQ=255)
                               gTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACggg                                                                 >  2:1374119/1‑100 (MQ=255)
                                  gagtgagtGCCCCAGTTGTGGTGCTGCGCCATGTTGTTAAGTGCTTTTTCCGGGAAGCCGATTGTCGGTATCTGGCGCATGTCACTGACGCTACGGgggt                                                              >  2:1231010/1‑100 (MQ=255)
                                  gAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGgtgt                                                              >  1:1093401/1‑100 (MQ=255)
                                   aGTGAATGCACCATTTCTGGTGCGGCGTCAGCAAGTTAAGTGATTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGCCGCTACGGGTGTc                                                             <  2:3898262/100‑1 (MQ=38)
                                     tGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCATTTTCCGGGAAGAAGTTTGTCGGTATCTGGCGCATCTCACTGACCCTAAGGGTGTCGa                                                           >  2:5246441/1‑100 (MQ=38)
                                      gAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGAt                                                          <  2:107672/100‑1 (MQ=255)
                                      gAATGCCACATTTCTGGTGCGGCGGCAGAAAGGTAAGTGCTATTTCGGGGAAGACGTGGGTGGGAATAGGGCGCATGTGCCTGGGGCTACGGgtgggtct                                                          >  2:4625023/1‑95 (MQ=255)
                                       aaTGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATa                                                         <  1:3624850/100‑1 (MQ=255)
                                        aTGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAg                                                        >  2:4387111/1‑100 (MQ=255)
                                          gCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTGGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggc                                                      >  2:2739904/1‑100 (MQ=255)
                                          gCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggc                                                      <  1:3448999/100‑1 (MQ=255)
                                           ccccATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcg                                                     >  1:5502333/1‑100 (MQ=255)
                                            cccATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgg                                                    <  1:2871021/100‑1 (MQ=255)
                                              cATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggc                                                  <  2:3468016/100‑1 (MQ=255)
                                              cATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggc                                                  <  1:4568153/100‑1 (MQ=255)
                                              cATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggc                                                  <  1:3331123/100‑1 (MQ=255)
                                              cATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggc                                                  <  2:4932618/100‑1 (MQ=255)
                                              cATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggc                                                  <  2:712175/100‑1 (MQ=255)
                                               aTTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGAAAGCCGTTTGTCGATATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgggcg                                                 >  2:5157619/1‑100 (MQ=39)
                                                  tCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTTGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGTGCGATa                                              <  2:2207173/100‑1 (MQ=255)
                                                  tCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATa                                              >  1:2707438/1‑100 (MQ=255)
                                                   cTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATaa                                             <  2:2826825/100‑1 (MQ=255)
                                                    tGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAg                                            <  1:148079/100‑1 (MQ=255)
                                                     ggTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGc                                           >  1:71081/1‑100 (MQ=255)
                                                       tGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGcca                                         >  2:2973937/1‑100 (MQ=255)
                                                             gTCAGGAAGTTAAGTTCTTTTTCCGGGCAGCCGTCTGTCGGCATCTGGCGCATCTCACGAACGCTACGGGGGTCGATAGGCGGACGATTAGCGACCAGGc                                   <  2:419759/100‑1 (MQ=255)
                                                             gTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGc                                   >  2:1733913/1‑100 (MQ=255)
                                                             gTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGc                                   >  1:4927286/1‑100 (MQ=255)
                                                               cAGGAAGTTAAGTGCTTTTTTTTTTTTTCTGTTTGTCGGTATCTGGCACATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATCCGCCACCAGGcgt                                 >  2:3800529/1‑98 (MQ=255)
                                                               cAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCtt                                 <  1:1409616/100‑1 (MQ=255)
                                                               cAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCtt                                 >  1:5501385/1‑100 (MQ=255)
                                                                aGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCttt                                >  2:3921505/1‑100 (MQ=255)
                                                                 ngAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                               >  1:410373/2‑100 (MQ=255)
                                                                 ggAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                               <  1:2768188/100‑1 (MQ=255)
                                                                 ggAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                               >  2:1199836/1‑100 (MQ=255)
                                                                 ggAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTc                               >  1:73477/1‑100 (MQ=255)
                                                                  gAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCa                              <  1:5291329/100‑1 (MQ=255)
                                                                       atcttGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCaaa                         <  2:4350997/96‑1 (MQ=255)
                                                                        aaGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGATTTAACAaaaa                        <  2:26467/100‑1 (MQ=255)
                                                                         aGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAg                       >  1:671463/1‑100 (MQ=255)
                                                                         aGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAg                       >  1:810425/1‑100 (MQ=255)
                                                                         aGTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAg                       >  1:1672335/1‑100 (MQ=255)
                                                                          gTGCTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGc                      >  1:1263794/1‑100 (MQ=255)
                                                                             cTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCTCATCTCACTGACGCTACGGGTGTCGATCGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACg                   >  2:2073029/1‑100 (MQ=38)
                                                                             cTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACg                   >  1:3386650/1‑100 (MQ=255)
                                                                             cTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACg                   <  2:76056/100‑1 (MQ=255)
                                                                             cTTTTTCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACg                   <  1:3795158/100‑1 (MQ=255)
                                                                              tttttCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGc                  <  1:1216000/100‑1 (MQ=255)
                                                                               ttttCCGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCa                 <  2:559768/100‑1 (MQ=255)
                                                                               ttttCCGGGAAGCCGTTTGTCGGTAGCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCa                 <  1:1819723/100‑1 (MQ=255)
                                                                               ttctCCGGGAAGCCGTTTGTGGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACTAACAGCACGCa                 <  2:3408518/97‑1 (MQ=255)
                                                                                tttCCGGGAAGCCGTTTNTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCACAAGCACGCAt                >  1:304844/1‑100 (MQ=255)
                                                                                    cGGGAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCAc            >  2:4499180/1‑100 (MQ=255)
                                                                                     gggAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACg           <  2:458494/100‑1 (MQ=255)
                                                                                      ggAAGCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACgn          <  1:4428592/100‑2 (MQ=255)
                                                                                          gCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGgatga      <  1:3534966/100‑1 (MQ=255)
                                                                                          gCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGgatga      >  1:1181226/1‑100 (MQ=255)
                                                                                          gCCGTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACACATCCACGgatga      <  2:3036009/100‑1 (MQ=255)
                                                                                             gTTTGTCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATcc   <  1:285006/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCT  >  NC_000913/1539261‑1539454

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