Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,281,011 A→G H642R (CAC→CGC)  narG → nitrate reductase 1 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,281,0110AG83.8% 590.9 / 102.6 266H642R (CAC→CGC) narGnitrate reductase 1 alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (29/14);  new base G (135/88);  total (164/102)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 270 total aligned reads displayed.

TTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGA  >  NC_000913/1280919‑1281101
                                                                                             |                                                                                          
  tCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCt                                                                                          <  1:5880821/93‑1 (MQ=255)
   cTATACCCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCtt                                                                                         <  1:1018570/93‑1 (MQ=255)
   cTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACtt                                                                                         >  2:1633954/1‑93 (MQ=255)
    tataCCCCCTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTg                                                                                        <  1:3986105/93‑1 (MQ=255)
     ataACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGa                                                                                       <  1:7213685/93‑1 (MQ=255)
     ataACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGa                                                                                       <  1:5732974/93‑1 (MQ=255)
      taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGAt                                                                                      >  1:2593246/1‑93 (MQ=255)
       aaCCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATc                                                                                     >  2:1396280/1‑93 (MQ=255)
        aCCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCg                                                                                    >  2:5616203/1‑93 (MQ=255)
          cACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGAc                                                                                  >  2:287199/1‑93 (MQ=255)
          cACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAAAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGAc                                                                                  >  1:6216452/1‑93 (MQ=255)
           aCTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACt                                                                                 >  2:4564287/1‑93 (MQ=255)
           aCTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACt                                                                                 >  2:6127195/1‑93 (MQ=255)
           aCTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACt                                                                                 >  1:3064244/1‑93 (MQ=255)
            cTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACtt                                                                                <  1:6919087/93‑1 (MQ=255)
            cTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACtt                                                                                >  1:7327679/1‑93 (MQ=255)
            cTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACtt                                                                                <  1:927409/93‑1 (MQ=255)
            cTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACtt                                                                                >  2:1558569/1‑93 (MQ=255)
             tCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACttt                                                                               >  2:2768168/1‑93 (MQ=255)
             tCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACttt                                                                               >  2:55115/1‑93 (MQ=255)
              ccagccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTa                                                                              <  1:4426022/93‑1 (MQ=255)
                agccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAAc                                                                            <  1:6393267/93‑1 (MQ=255)
                 gccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACg                                                                           >  2:6729451/1‑93 (MQ=255)
                 gccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACg                                                                           <  2:657763/93‑1 (MQ=255)
                  ccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGt                                                                          >  1:1853383/1‑93 (MQ=255)
                  ccagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGt                                                                          >  1:3884088/1‑93 (MQ=255)
                   cagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTc                                                                         <  1:5696107/93‑1 (MQ=255)
                   cagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTc                                                                         >  1:3296621/1‑93 (MQ=255)
                    agTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTcc                                                                        >  2:3840246/1‑93 (MQ=255)
                    agTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTcc                                                                        >  1:6574943/1‑93 (MQ=255)
                     gTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCg                                                                       <  2:3546793/93‑1 (MQ=255)
                      tGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGt                                                                      >  1:294353/1‑93 (MQ=255)
                      tGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGt                                                                      >  1:3595164/1‑93 (MQ=255)
                        gCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGc                                                                    >  1:2366755/1‑93 (MQ=255)
                        gCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGc                                                                    >  2:7140010/1‑93 (MQ=255)
                           ttATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGGTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGa                                                                 >  1:4521234/1‑93 (MQ=255)
                           ttATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGa                                                                 >  2:479081/1‑93 (MQ=255)
                           ttATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGa                                                                 >  1:6264190/1‑93 (MQ=255)
                            tATGAAACCGTCACGGCGGAAAGAGTTGCTGTCACCGAGGGCGGACAAATCCCGCTAAACCGGACGCTTGATCGACTTTAAAGTCCGGGCGGa                                                                 >  1:7457785/1‑93 (MQ=255)
                            tATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGaa                                                                <  2:5771522/93‑1 (MQ=255)
                            tATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGaa                                                                >  2:4663664/1‑93 (MQ=255)
                             aTGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAAc                                                               <  1:4254121/93‑1 (MQ=255)
                              tGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACg                                                              >  2:6977152/1‑93 (MQ=255)
                              tGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACg                                                              >  1:1918371/1‑93 (MQ=255)
                                aaaCCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCa                                                            <  2:2204152/93‑1 (MQ=255)
                                 aaCCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCAt                                                           >  2:2523609/1‑93 (MQ=255)
                                 aaCCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCAt                                                           >  1:2518334/1‑93 (MQ=255)
                                  aCCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATg                                                          <  2:4060744/93‑1 (MQ=255)
                                   ccGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATgg                                                         <  2:4325050/93‑1 (MQ=255)
                                   ccGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATgg                                                         <  1:2216444/93‑1 (MQ=255)
                                      tCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGggct                                                      >  1:459540/1‑93 (MQ=255)
                                        aCGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGggctgg                                                    <  1:1696480/93‑1 (MQ=255)
                                         cggcggAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGggctggc                                                   <  2:3448106/93‑1 (MQ=255)
                                          ggcggAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGggctggct                                                  <  1:471575/93‑1 (MQ=255)
                                            cggAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGc                                                <  1:864763/93‑1 (MQ=255)
                                            cggAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGc                                                >  2:3959378/1‑93 (MQ=255)
                                              gAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCg                                              >  1:1197670/1‑93 (MQ=255)
                                                   agTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTg                                          <  2:3348091/93‑1 (MQ=255)
                                                    gTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGc                                         >  1:942513/1‑93 (MQ=255)
                                                     ttGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCa                                        <  1:6810910/93‑1 (MQ=255)
                                                         tGTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGc                                    >  2:246640/1‑93 (MQ=255)
                                                         tGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGc                                    >  2:6478957/1‑93 (MQ=255)
                                                          gTCACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCa                                   >  1:368106/1‑93 (MQ=255)
                                                           nnaCCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAg                                  <  1:3246019/91‑1 (MQ=255)
                                                            cACCGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGt                                 <  2:741662/93‑1 (MQ=255)
                                                              ccGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTa                               <  1:495293/93‑1 (MQ=255)
                                                               cGATGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAg                              <  1:5368926/93‑1 (MQ=255)
                                                                 aTGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGc                            >  2:4401949/1‑93 (MQ=255)
                                                                 aTGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGc                            >  2:2708858/1‑93 (MQ=255)
                                                                  tGGCGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCa                           >  1:6868091/1‑93 (MQ=255)
                                                                     cGGACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTa                        <  2:115325/93‑1 (MQ=255)
                                                                      ggACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTaa                       >  1:1433138/1‑93 (MQ=255)
                                                                      ggACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTaa                       <  2:3799960/93‑1 (MQ=255)
                                                                       gACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAAc                      >  2:3153624/1‑93 (MQ=255)
                                                                       gACAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAAc                      >  2:1478865/1‑93 (MQ=255)
                                                                        aCAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAAcc                     <  1:3383562/93‑1 (MQ=255)
                                                                        aCAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAAcc                     <  1:5536913/93‑1 (MQ=255)
                                                                        aCAAATCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAAcc                     <  1:6059115/93‑1 (MQ=255)
                                                                            aTCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCt                 <  1:7367804/93‑1 (MQ=255)
                                                                            aTCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCt                 >  1:1489943/1‑93 (MQ=255)
                                                                             tCCCGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTg                >  2:6414233/1‑93 (MQ=255)
                                                                               ccGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGAc              <  2:3771663/93‑1 (MQ=255)
                                                                               ccGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGAc              >  2:2057065/1‑93 (MQ=255)
                                                                                cGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACt             >  1:12607/1‑93 (MQ=255)
                                                                                cGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACt             >  1:4327768/1‑93 (MQ=255)
                                                                                cGCTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACt             >  1:5033166/1‑93 (MQ=255)
                                                                                  cTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTAt           >  2:822571/1‑93 (MQ=255)
                                                                                  cTATACCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTAt           >  2:4913860/1‑93 (MQ=255)
                                                                                  cTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTAt           >  2:5349619/1‑93 (MQ=255)
                                                                                   tataCCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATc          >  2:3399147/1‑93 (MQ=255)
                                                                                   tataCCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATc          >  1:6635718/1‑93 (MQ=255)
                                                                                    ataCCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCg         >  1:5176398/1‑93 (MQ=255)
                                                                                    ataCCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCg         <  1:2709112/93‑1 (MQ=255)
                                                                                     taCCGGACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGc        >  2:1803464/1‑93 (MQ=255)
                                                                                     taCCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGc        <  1:6479260/93‑1 (MQ=255)
                                                                                       ccGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTg      <  1:3000543/93‑1 (MQ=255)
                                                                                         ggACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGc    >  2:6846954/1‑93 (MQ=255)
                                                                                         ggACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGc    >  2:7004947/1‑93 (MQ=255)
                                                                                          gACGCTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCg   >  2:5903400/1‑93 (MQ=255)
                                                                                           acacTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGa  >  1:4089168/1‑93 (MQ=255)
                                                                                             |                                                                                          
TTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGA  >  NC_000913/1280919‑1281101

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