Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,105 C→G 5.4% A501G (GCC→GGC)  narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,1050CG 5.4% 456.9 / 12.6 129A501G (GCC→GGC) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (57/65);  new base G (3/4);  total (60/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-02

Only 100 of 132 total aligned reads displayed.

TACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAA  >  NC_000913/1284957‑1285249
                                                                                                                                                    |                                                                                                                                                
tACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGc                                                                                                                                                  >  2:126474/1‑149 (MQ=255)
   gAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCAt                                                                                                                                               >  2:405154/1‑149 (MQ=255)
      gATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGa                                                                                                                                            <  1:611812/149‑1 (MQ=255)
       aTCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAt                                                                                                                                           <  1:126474/149‑1 (MQ=255)
       aTCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAt                                                                                                                                           >  2:551730/1‑149 (MQ=255)
        tCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAtg                                                                                                                                          >  3:990502/1‑149 (MQ=255)
         cggtttGTGGTGCCGAGTAGTCATCGTGAACTCGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGGTTTACCTTTGGGGATGGCTGCCACGGTTCAGATACCAAATTCAATCTATTCAACAGCCGTCGTATCGATGCCATCGAtgt                                                                                                                                         <  3:812125/146‑1 (MQ=255)
         cGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATCGAtac                                                                                                                                         <  1:748583/149‑3 (MQ=255)
         cGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAtgt                                                                                                                                         <  1:405154/149‑1 (MQ=255)
         cGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAtgt                                                                                                                                         <  1:513583/149‑1 (MQ=255)
         cGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGAtgt                                                                                                                                         >  4:470844/1‑149 (MQ=255)
           cTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATCGAtacga                                                                                                                                       <  2:1408026/149‑5 (MQ=255)
            tttGTGGTGCCGAGTAGTCTTCGTGAACTGGCATGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGAc                                                                                                                                      >  3:1160968/1‑149 (MQ=255)
            tttGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGAc                                                                                                                                      >  4:146665/1‑149 (MQ=255)
                    gCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCaaaa                                                                                                                              <  1:773022/149‑1 (MQ=255)
                            gtCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCg                                                                                                                      >  2:1214968/1‑149 (MQ=255)
                             tCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGc                                                                                                                     >  2:395512/1‑149 (MQ=255)
                               aTCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCAt                                                                                                                   >  3:136949/1‑149 (MQ=255)
                                          gCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGaa                                                                                                        >  3:1401937/1‑149 (MQ=255)
                                          gCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGaa                                                                                                        <  4:336984/149‑1 (MQ=255)
                                           tccGGGAAGGCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAAc                                                                                                       >  4:607019/3‑149 (MQ=255)
                                           cACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAAc                                                                                                       <  1:1326793/149‑1 (MQ=255)
                                            aCGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGa                                                                                                                          <  3:610392/129‑1 (MQ=255)
                                            aCGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGa                                                                                                                          >  4:610392/1‑129 (MQ=255)
                                            aCGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACt                                                                                                      >  2:76489/1‑149 (MQ=255)
                                            aCGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACt                                                                                                      >  3:592051/1‑149 (MQ=255)
                                               ggAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTc                                                                                                     <  4:871245/147‑1 (MQ=255)
                                               ggAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTc                                                                                                     >  3:871245/1‑147 (MQ=255)
                                                gAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGAt                                                                                                            >  1:812136/1‑139 (MQ=255)
                                                gAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGAt                                                                                                            <  2:812136/139‑1 (MQ=255)
                                                 aaGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGa                                                                                                 <  1:395512/149‑1 (MQ=255)
                                                  aGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGAt                                                                                                <  4:369720/149‑1 (MQ=255)
                                                      ttCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCTTGATTGTa                                                                                            >  3:61282/1‑149 (MQ=255)
                                                      ttCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGa                                                                                                 >  1:967566/1‑144 (MQ=255)
                                                      ttCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGa                                                                                                 <  2:967566/144‑1 (MQ=255)
                                                        cccGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATc                                                                                          <  1:63156/149‑1 (MQ=255)
                                                          cGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATcgc                                                                                        <  2:814991/149‑1 (MQ=255)
                                                             aaaaaaaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGTATCCATGATCGAACTCGTGATTGTATCGCGtc                                                                                     <  4:592051/147‑1 (MQ=255)
                                                                aaaaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTcc                                                                                  <  3:146665/149‑1 (MQ=255)
                                                                 aaaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATc                                                                                                                                              <  4:431221/88‑1 (MQ=255)
                                                                 aaaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATc                                                                                                                                              >  3:431221/1‑88 (MQ=255)
                                                                  aaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCtt                                                                                <  4:136949/149‑1 (MQ=255)
                                                                  aaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCtt                                                                                <  1:1214968/149‑1 (MQ=255)
                                                                  aaTGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGATCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTCACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTg                                                                                                  <  2:96781/131‑1 (MQ=255)
                                                                     ggCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGc                                                                                                                                                  >  4:282125/1‑80 (MQ=255)
                                                                        tGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATcgc                                                                                        <  2:854233/135‑1 (MQ=255)
                                                                        tGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATcgc                                                                                        >  1:854233/1‑135 (MQ=255)
                                                                        tGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAatat                                                                          >  1:1087246/1‑149 (MQ=255)
                                                                             cTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGa                                                                     <  4:93209/149‑1 (MQ=255)
                                                                                   cTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATc                                                                                                           >  1:1092293/1‑105 (MQ=255)
                                                                                    tttGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGcc                                                              >  3:446242/1‑149 (MQ=255)
                                                                                     ttGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCt                                                             <  1:76489/149‑1 (MQ=255)
                                                                                        gTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                          <  4:1160968/149‑1 (MQ=255)
                                                                                         tGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                         <  2:762351/149‑1 (MQ=255)
                                                                                          gATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                        <  2:1087246/149‑1 (MQ=255)
                                                                                            tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGaa                                                                                                        >  4:285145/1‑99 (MQ=255)
                                                                                            tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGaa                                                                                                        <  3:285145/99‑1 (MQ=255)
                                                                                            aggCTGCCACGGTGCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAGCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgca                                                      >  4:945398/2‑149 (MQ=255)
                                                                                             ggCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGcc                                                              >  4:1889/1‑140 (MQ=255)
                                                                                             ggCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGcc                                                              <  3:1889/140‑1 (MQ=255)
                                                                                             ggCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcag                                                     <  4:446242/149‑1 (MQ=255)
                                                                                                  ccACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGgcag                                                     <  4:382361/144‑1 (MQ=255)
                                                                                                   cACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATcc                                                                        <  4:700122/124‑1 (MQ=255)
                                                                                                   cACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATcc                                                                        >  3:700122/1‑124 (MQ=255)
                                                                                                        ttCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcagc                                                    <  3:56365/139‑1 (MQ=255)
                                                                                                        ttCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcagc                                                    >  4:56365/1‑139 (MQ=255)
                                                                                                           aGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATc                                                                                          <  4:1287915/98‑1 (MQ=255)
                                                                                                           aGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGGCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATc                                                                                          >  3:1287915/1‑98 (MQ=255)
                                                                                                            gATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGa                                                                     <  4:224982/118‑1 (MQ=255)
                                                                                                            gATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGa                                                                     >  3:224982/1‑118 (MQ=255)
                                                                                                             aTACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgca                                                      >  4:864860/1‑132 (MQ=255)
                                                                                                             aTACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgca                                                      <  3:864860/132‑1 (MQ=255)
                                                                                                              tACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGtt                                    >  2:705456/1‑149 (MQ=255)
                                                                                                                ccAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGAtgctgc                                                               >  2:117520/1‑120 (MQ=255)
                                                                                                                ccAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGAtgctgc                                                               <  1:117520/120‑1 (MQ=255)
                                                                                                                 cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATcgc                                                                                        <  2:1063216/94‑1 (MQ=255)
                                                                                                                 cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATcgc                                                                                        >  1:1063216/1‑94 (MQ=255)
                                                                                                                 cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGa                                 >  1:957959/1‑149 (MQ=255)
                                                                                                                  aaaTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTAt                                                                                           <  3:341631/90‑1 (MQ=255)
                                                                                                                  aaaTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGt                                     >  4:1321035/1‑144 (MQ=255)
                                                                                                                  aaaTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGt                                     <  3:1321035/144‑1 (MQ=255)
                                                                                                                    aTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGa                                                                                                 >  2:167367/1‑82 (MQ=255)
                                                                                                                    aTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGa                                                                                                 <  1:167367/82‑1 (MQ=255)
                                                                                                                    aTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                         <  4:1296944/122‑1 (MQ=255)
                                                                                                                    aTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                         >  3:1296944/1‑122 (MQ=255)
                                                                                                                    aTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGttt                                   >  4:1245943/1‑144 (MQ=255)
                                                                                                                        aaTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGa                                                                              <  3:441195/97‑1 (MQ=255)
                                                                                                                        aaTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGAtt                          >  4:1072761/1‑149 (MQ=255)
                                                                                                                         aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTg                         >  4:800874/1‑149 (MQ=255)
                                                                                                                             gTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGcc                                                              >  1:1058929/1‑108 (MQ=255)
                                                                                                                             gTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGcc                                                              <  2:1058929/108‑1 (MQ=255)
                                                                                                                              ttCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAg                                <  3:573356/137‑1 (MQ=255)
                                                                                                                                cAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCTATTGCCGCGTc                  <  3:945398/149‑1 (MQ=255)
                                                                                                                                  aCAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCcgcg                    <  2:1125998/145‑1 (MQ=255)
                                                                                                                                    aGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGaaa              <  3:800874/149‑1 (MQ=255)
                                                                                                                                           acaTCGATGGCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGc       >  3:289869/3‑149 (MQ=255)
                                                                                                                                            tATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGcc      >  3:1257590/1‑149 (MQ=255)
                                                                                                                                                gATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGAtt                          >  3:246237/1‑125 (MQ=255)
                                                                                                                                                  tGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                          >  1:1384234/1‑107 (MQ=255)
                                                                                                                                                  tGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                          <  2:1384234/107‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAA  >  NC_000913/1284957‑1285249

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