Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,261 Δ1 bp 6.4% coding (123/711 nt) narJ → molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,2610T. 6.4% 403.9 / 19.2 110coding (123/711 nt)narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (52/51);  new base . (4/3);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

Only 100 of 110 total aligned reads displayed.

GCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCC  >  NC_000913/1285120‑1285403
                                                                                                                                             |                                                                                                                                               
gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                                                                                         >  3:1327530/1‑149 (MQ=255)
     aCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtt                                                                                                                                    <  4:1257590/149‑1 (MQ=255)
              cATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTc                                                                                                                                 <  4:491639/143‑1 (MQ=255)
              cATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTc                                                                                                                                 >  3:491639/1‑143 (MQ=255)
               aTCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGcgcg                                                                                                                          <  4:289869/149‑1 (MQ=255)
                     gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                    <  2:532235/149‑1 (MQ=255)
                     gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                    <  4:767765/149‑1 (MQ=255)
                     gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                    >  4:1335497/1‑149 (MQ=255)
                           aCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAAcgacga                                                                                                              <  3:1335497/149‑1 (MQ=255)
                             tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                   <  2:376475/142‑1 (MQ=255)
                             tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                   >  1:376475/1‑142 (MQ=255)
                               gTGATTGTATCGCGTCTCCTTGACTATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGa                                                                                                          >  4:89382/1‑149 (MQ=255)
                               gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgc                                                                                                                                               <  1:1117026/113‑1 (MQ=255)
                               gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgc                                                                                                                                               >  2:1117026/1‑113 (MQ=255)
                               gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAAcgacga                                                                                                              >  2:601529/1‑145 (MQ=255)
                               gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAAcgacga                                                                                                              <  1:601529/145‑1 (MQ=255)
                                        tCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGTATCCGCTCGa                                                                                                 <  3:89382/149‑1 (MQ=255)
                                                  ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                       >  2:1037968/1‑149 (MQ=255)
                                                          ccGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgc                                                                                                                                               <  4:1182682/86‑1 (MQ=255)
                                                          ccGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgc                                                                                                                                               >  3:1182682/1‑86 (MQ=255)
                                                              aTGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                           <  1:1394286/149‑1 (MQ=255)
                                                                  tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCt                                                                                                                               <  1:723793/93‑1 (MQ=255)
                                                                  tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCt                                                                                                                               >  2:723793/1‑93 (MQ=255)
                                                                  tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTg                                                                                                                              <  1:67022/94‑1 (MQ=255)
                                                                  tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                       >  3:520529/1‑133 (MQ=255)
                                                                  tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                       <  4:520529/133‑1 (MQ=255)
                                                                     cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtttt                                                                                                                                  <  3:193181/87‑1 (MQ=255)
                                                                     cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtttt                                                                                                                                  >  4:193181/1‑87 (MQ=255)
                                                                       taaggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                   >  2:146383/4‑100 (MQ=255)
                                                                       taaggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                   <  1:146383/97‑1 (MQ=255)
                                                                       tATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc                                                                  >  4:813081/1‑149 (MQ=255)
                                                                           gcagcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACa                                                                                >  1:761959/1‑131 (MQ=255)
                                                                           gcagcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACa                                                                                <  2:761959/131‑1 (MQ=255)
                                                                           gcagcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTg                                                              >  4:1064049/1‑149 (MQ=255)
                                                                               caTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGcgc                                                          >  4:437677/1‑147 (MQ=255)
                                                                                aTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                             >  1:192344/1‑145 (MQ=255)
                                                                                aTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                             <  2:192344/145‑1 (MQ=255)
                                                                                 tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgac                                                                                                                <  1:670807/93‑1 (MQ=255)
                                                                                 tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgac                                                                                                                >  2:670807/1‑93 (MQ=255)
                                                                                 tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTc                                                                     >  1:1072866/1‑136 (MQ=255)
                                                                                 tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTc                                                                     <  2:1072866/136‑1 (MQ=255)
                                                                                 tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg                                                        >  1:152665/1‑149 (MQ=255)
                                                                                   aacaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc                                                                  >  2:578073/1‑137 (MQ=255)
                                                                                   aacaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc                                                                  <  1:578073/137‑1 (MQ=255)
                                                                                   aacaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGcc                                                      >  2:725357/1‑149 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGTCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGAt                                                                                                             >  4:853607/1‑93 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                    <  3:839354/86‑1 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                    >  4:839354/1‑86 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGAt                                                                                                             <  3:853607/93‑1 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                           >  2:313351/1‑127 (MQ=255)
                                                                                    acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAcc                                                                 >  4:1211102/1‑137 (MQ=255)
                                                                                     caaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCa                                                                                     >  4:319542/1‑116 (MQ=255)
                                                                                     caaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCa                                                                                     <  3:319542/116‑1 (MQ=255)
                                                                                       agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCGCGGCCACgg                                                  <  3:437677/149‑2 (MQ=255)
                                                                                        gagaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                      <  2:35012/128‑1 (MQ=255)
                                                                                        gagaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                      >  1:35012/1‑128 (MQ=255)
                                                                                         agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGAt                                                                                                <  3:186035/101‑1 (MQ=255)
                                                                                         agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGAt                                                                                                >  4:186035/1‑101 (MQ=255)
                                                                                         agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCa                                                >  1:122306/1‑149 (MQ=255)
                                                                                             gTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                      <  4:1294156/123‑1 (MQ=255)
                                                                                             gTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                      >  3:1294156/1‑123 (MQ=255)
                                                                                                   ggCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                      <  2:152665/149‑1 (MQ=255)
                                                                                                        ttGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                 <  2:122306/149‑1 (MQ=255)
                                                                                                        ttGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                 <  1:1037968/149‑1 (MQ=255)
                                                                                                           ccGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              <  3:813081/149‑1 (MQ=255)
                                                                                                           ccGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              <  1:725357/149‑1 (MQ=255)
                                                                                                             gcgTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGATGATGTATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                            <  3:1064049/149‑1 (MQ=255)
                                                                                                                  gAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAttt                                                                                                                                   >  3:980498/1‑41 (MQ=255)
                                                                                                                  gAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                       >  4:1428275/1‑149 (MQ=255)
                                                                                                                       aTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                               <  3:1315530/120‑1 (MQ=255)
                                                                                                                       aTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                               >  4:1315530/1‑120 (MQ=255)
                                                                                                                       aTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                    >  2:690330/1‑131 (MQ=255)
                                                                                                                       aTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                    <  1:690330/131‑1 (MQ=255)
                                                                                                                       aTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa                  >  1:573704/1‑149 (MQ=255)
                                                                                                                          tGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              <  1:132285/134‑1 (MQ=255)
                                                                                                                             cAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc            >  4:171255/1‑149 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt                                                                       <  2:27552/89‑1 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt                                                                       >  1:27552/1‑89 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                     <  1:920598/123‑1 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                     >  2:920598/1‑123 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg           >  1:357241/1‑149 (MQ=255)
                                                                                                                              aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg           >  1:314810/1‑149 (MQ=255)
                                                                                                                                aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGAttt                                                                                                                      <  3:783409/40‑1 (MQ=255)
                                                                                                                                aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGAttt                                                                                                                      >  4:783409/1‑40 (MQ=255)
                                                                                                                                aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTATTGCTGTTTGAACATGTGCACGGCGAATCACGCGAc         >  3:973511/1‑149 (MQ=255)
                                                                                                                                    aaaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGtt                                           >  4:887292/3‑111 (MQ=255)
                                                                                                                                     gATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCGCGGCCACggt                                                 <  3:1160678/104‑3 (MQ=255)
                                                                                                                                     gATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCgg    >  3:496252/1‑149 (MQ=255)
                                                                                                                                     aaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTg                                                              <  3:368524/90‑1 (MQ=255)
                                                                                                                                     aaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTg                                                              >  4:368524/2‑91 (MQ=255)
                                                                                                                                     aaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAAt                 <  2:776955/135‑1 (MQ=255)
                                                                                                                                     aaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAAt                 >  1:776955/2‑136 (MQ=255)
                                                                                                                                      aTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc   >  2:1249351/1‑149 (MQ=255)
                                                                                                                                       tGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc  <  3:1428275/149‑1 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                           <  4:28590/73‑1 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                           >  3:28590/1‑73 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              <  3:1227296/118‑1 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              >  4:795496/1‑118 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              >  4:1227296/1‑118 (MQ=255)
                                                                                                                                          ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                              <  3:795496/118‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCC  >  NC_000913/1285120‑1285403

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