Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,541,077 G→C 15.0% A505G (GCC→GGC)  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,541,0770GC15.0% 432.8 / 68.2 167A505G (GCC→GGC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (81/60);  new base C (9/16);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 5.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 169 total aligned reads displayed.

ATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGC  >  NC_000913/1540981‑1541176
                                                                                                |                                                                                                   
   gAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                               <  2:2689559/100‑1 (MQ=255)
   gAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                               <  2:19549/100‑1 (MQ=255)
   gAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                               <  1:2222910/100‑1 (MQ=255)
   gAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTgggg                                                                                               <  2:2210623/100‑1 (MQ=255)
     gcgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtg                                                                                             >  1:3612066/1‑100 (MQ=255)
      ttccgatctgTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtga                                                                                            <  2:3504973/91‑2 (MQ=255)
      cgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                            <  2:441873/100‑1 (MQ=255)
      cgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                            >  2:2009452/1‑100 (MQ=255)
       gcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgan                                                                                           <  2:174551/100‑3 (MQ=255)
       gcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgag                                                                                           <  2:2199376/100‑3 (MQ=255)
       gcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgag                                                                                           <  1:2519441/100‑3 (MQ=255)
       gcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgag                                                                                           <  1:859158/100‑3 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                                                          >  2:3543163/1‑100 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                          <  2:983245/100‑4 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                          <  1:3692852/100‑4 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                          <  2:3300918/100‑4 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                          <  1:1825233/100‑4 (MQ=255)
        cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                          <  2:16831/100‑4 (MQ=255)
         ccttGCGTTTTATGGGCGGTATTGGCAAATTCTCGGGCGATGGTTTCAATATACTGGCGCGGCACGCAGGTAATTTGCTCTCCCCTGGCTgggg                                                                                               <  1:2426284/91‑1 (MQ=255)
         ccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagca                                                                                         <  1:3494278/100‑5 (MQ=255)
          cATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACgg                                                                                        >  2:400460/1‑100 (MQ=255)
            tGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGGGCGGCACGCCGGTAATTTGCTCACGCCAGGATGGGGTGTACg                                                                                         >  1:3134187/1‑97 (MQ=255)
             gCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                     <  2:2392154/100‑1 (MQ=255)
              cGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                    <  2:3433755/100‑1 (MQ=255)
               gTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTCCGGTTTg                                                                                    >  1:1022729/1‑99 (MQ=255)
               gTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGa                                                                                   <  1:1543654/100‑1 (MQ=255)
               gTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgaagatcggaa                                                                                   >  1:3504952/1‑90 (MQ=255)
                ttttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGATCAGCCCAGGCTGGGGTGTGCGGTTTGAt                                                                                  >  1:1484186/1‑100 (MQ=255)
                 tttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAtt                                                                                 <  2:1353121/100‑1 (MQ=255)
                  ttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAtt                                                                                 <  1:1260158/99‑1 (MQ=255)
                   tATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                               <  2:735571/100‑1 (MQ=255)
                    aTGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCa                                                                              <  2:886220/100‑1 (MQ=255)
                     tGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                             <  2:1036138/100‑1 (MQ=255)
                     tGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                             <  2:901814/100‑1 (MQ=255)
                     tGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                             <  1:3127408/100‑1 (MQ=255)
                      gggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGc                                                                            >  1:442104/1‑100 (MQ=255)
                      gggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGc                                                                            <  2:1613759/100‑1 (MQ=255)
                       ggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCg                                                                           <  1:1368216/100‑1 (MQ=255)
                       ggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCg                                                                           <  2:807641/100‑1 (MQ=255)
                             aTCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCt                                                                     <  2:3167759/100‑1 (MQ=255)
                             aTCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTAGGGTGTACGGTTTGATTTCAGCGTAAt                                                                       >  1:2300360/1‑98 (MQ=255)
                               cGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                   >  2:560106/1‑100 (MQ=255)
                                ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                   >  2:949883/1‑99 (MQ=255)
                                ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTc                                                                  >  2:455691/1‑100 (MQ=255)
                                ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTc                                                                  >  2:3627864/1‑100 (MQ=255)
                                  cAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCg                                                                 >  2:2246312/1‑99 (MQ=255)
                                  cAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGc                                                                <  1:3648570/100‑1 (MQ=255)
                                  cAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGc                                                                <  2:800314/100‑1 (MQ=255)
                                         aCGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGttt                                                         <  2:13582/100‑1 (MQ=255)
                                           gggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGttttn                                                       <  1:3055622/100‑2 (MQ=255)
                                            ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGttttaa                                                      >  1:2156400/1‑98 (MQ=255)
                                               gATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCt                                                   >  2:366447/1‑100 (MQ=255)
                                               gATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCt                                                   >  1:2706728/1‑100 (MQ=255)
                                                 tGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATGtt                                                  >  1:3626270/1‑99 (MQ=255)
                                                    tttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAg                                              <  2:3131526/100‑1 (MQ=255)
                                                    tttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAg                                              <  1:566320/100‑1 (MQ=255)
                                                    tttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCAATGCAGGGGTGTACGGTTTGATTTCAGCATAATCTTTCGCACTGTTTTCATCTTc                                                 >  1:919377/1‑97 (MQ=255)
                                                       cAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGccc                                           <  2:1563026/100‑1 (MQ=255)
                                                       cAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGccc                                           >  2:2276322/1‑100 (MQ=255)
                                                         atatACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCcgc                                         >  2:3022043/1‑100 (MQ=255)
                                                          tataCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCcgcg                                        >  2:2152370/1‑100 (MQ=255)
                                                              cTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCg                                    >  2:1466669/1‑100 (MQ=255)
                                                              cTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCg                                    >  1:1578995/1‑100 (MQ=255)
                                                                 gcgcgcCACGGCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATTTTTCGCACTGTTTTCATCGTCCAGCCCGCGATCGAgg                                 <  1:1181628/100‑1 (MQ=255)
                                                                 gcgcGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAgg                                 <  1:2915037/100‑1 (MQ=255)
                                                                  cgcgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGc                                >  2:2641444/1‑100 (MQ=255)
                                                                  cgcgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGc                                <  2:1438390/100‑1 (MQ=255)
                                                                   gcgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGcc                               <  2:1695218/100‑1 (MQ=255)
                                                                      gCACGCCGGTAATTTGCGCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAAGCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                            >  1:2432942/1‑100 (MQ=255)
                                                                        aCGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAAt                          >  1:962472/1‑100 (MQ=255)
                                                                         cGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCACGATCGAGGCCGTAAtt                         >  2:2939094/1‑100 (MQ=255)
                                                                          gCCGGTAATTTGCTCAACACAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCCTCTTCCAGCCCGCGATAGAGGCCGTaa                           <  1:2470321/97‑1 (MQ=255)
                                                                             ggTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                     <  2:257278/100‑1 (MQ=255)
                                                                             ggTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                     <  1:2152360/100‑1 (MQ=255)
                                                                              gTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCa                    >  2:59133/1‑100 (MQ=255)
                                                                               tAATTTGCTCACCCCAGGCTGGGGTGTACGGGTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAg                   <  2:1072383/100‑1 (MQ=255)
                                                                                 aaccgtacACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTTGCCAGTa                 <  2:1997302/93‑1 (MQ=255)
                                                                                   ttGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTAcc               >  2:1272148/1‑100 (MQ=255)
                                                                                   ccgtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCCTTCGCACTGTTTTCATCTTCCCGCCCGCGATCGAGGCCGTAATTCGCCAGTAcc               >  1:1128538/6‑100 (MQ=255)
                                                                                    tGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCa              >  2:344131/1‑100 (MQ=255)
                                                                                    cgtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCa              >  2:3204212/5‑100 (MQ=255)
                                                                                     gtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaa             >  2:2786425/4‑100 (MQ=255)
                                                                                     gCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaa             >  1:725512/1‑100 (MQ=255)
                                                                                     gCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGGAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaa             <  1:3777197/100‑1 (MQ=255)
                                                                                      cTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCATCCCACGATCGAGGCCGTAATTCGCCAGTACCaaa            >  1:328949/1‑100 (MQ=255)
                                                                                      cTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaaa            >  2:3730359/1‑100 (MQ=255)
                                                                                      cTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaaa            >  1:2935998/1‑100 (MQ=255)
                                                                                        cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTACTCTTTCGCACTGTTTTCATCTTTCAGCCCGCGATCGAGGCAGTAATTCGCCAGTACCAAATc          >  1:375849/1‑100 (MQ=255)
                                                                                        cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATc          <  1:1668090/100‑1 (MQ=255)
                                                                                          ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        >  1:3306359/1‑100 (MQ=255)
                                                                                          ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        >  2:3795062/1‑100 (MQ=255)
                                                                                          ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        >  2:650242/1‑100 (MQ=255)
                                                                                          ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        <  1:560099/100‑1 (MQ=255)
                                                                                          ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        >  2:778296/1‑100 (MQ=255)
                                                                                          ccccAGGATGGAGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt        >  1:2003346/1‑100 (MQ=255)
                                                                                          ccccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCa         >  1:664913/1‑99 (MQ=255)
                                                                                           cccAGGCTGGCGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATa       <  1:2509896/100‑1 (MQ=255)
                                                                                            ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa      >  2:3084376/1‑100 (MQ=255)
                                                                                             cAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa      >  1:1854173/1‑99 (MQ=255)
                                                                                                gCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGc  >  2:2128935/1‑100 (MQ=255)
                                                                                                |                                                                                                   
ATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGC  >  NC_000913/1540981‑1541176

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