Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,265:1 +A 6.2% coding (127/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,2651.A 6.2% 265.9 / 12.7 64coding (127/711 nt)narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base . (35/25);  new base A (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

CAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC  >  NC_000913/1285118‑1285407
                                                                                                                                                    |                                                                                                                                              
cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑t                                                                                                                                               >  4:638403/1‑149 (MQ=255)
cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑t                                                                                                                                               <  2:48124/149‑1 (MQ=255)
 aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑Tg                                                                                                                                              >  2:413945/1‑149 (MQ=255)
        cGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAtcc                                                                                                                                       <  2:572746/149‑3 (MQ=255)
                  tCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGa                                                                                                                             <  1:195953/149‑1 (MQ=255)
                     aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGAttt                                                                                                                          >  3:330094/1‑149 (MQ=255)
                       gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTaa                                                                                                                        <  4:490877/149‑1 (MQ=255)
                        aTCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAac                                                                                                                       >  1:295289/1‑149 (MQ=255)
                           gAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                                                                                                         >  4:261798/1‑128 (MQ=255)
                           gAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                                                                                                         <  3:261798/128‑1 (MQ=255)
                            aaCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAacgacg                                                                                                                   >  1:131296/1‑149 (MQ=255)
                               tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATg                                                                                                                >  3:77113/1‑149 (MQ=255)
                                      tGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGc                                                                                                         <  2:621514/149‑1 (MQ=255)
                                         aTCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCg                                                                                                      >  3:43890/1‑149 (MQ=255)
                                              gTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGcc                                                                                                 <  3:551311/149‑1 (MQ=255)
                                                 ttcTTAAATACCCGGATGCTGCCTTAAGGCAGCATAAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAg                                                                                              >  2:372599/3‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  1:271784/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  3:30695/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  2:481984/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATATAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  1:177511/1‑149 (MQ=255)
                                                      gAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCa                                                                                         <  1:645414/149‑1 (MQ=255)
                                                             cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGc                                                                                  <  1:31007/149‑1 (MQ=255)
                                                             cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGc                                                                                  <  4:24871/149‑1 (MQ=255)
                                                              ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                 >  4:55812/1‑149 (MQ=255)
                                                              ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                 >  4:262868/1‑149 (MQ=255)
                                                                aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGTGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGAt                                                                                                                            <  3:460384/104‑1 (MQ=255)
                                                                aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGAt                                                                                                                            >  4:460384/1‑104 (MQ=255)
                                                                     gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                          <  4:43890/149‑1 (MQ=255)
                                                                       cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCg                                                                        >  2:500580/1‑149 (MQ=255)
                                                                            ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGt                                                                   <  2:511315/149‑1 (MQ=255)
                                                                                  aTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgc                                                             >  2:426160/1‑149 (MQ=255)
                                                                                   tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc                                                                                     >  1:199500/1‑124 (MQ=255)
                                                                                   tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc                                                                                     <  2:199500/124‑1 (MQ=255)
                                                                                       caaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCAc                                                        >  3:204349/1‑149 (MQ=255)
                                                                                        aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGACGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACg                                                       >  4:419086/1‑149 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                               <  4:9985/124‑1 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                               >  3:9985/1‑124 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGt                                                      >  4:273487/1‑149 (MQ=255)
                                                                                                    aGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTg                                           >  4:419166/1‑149 (MQ=255)
                                                                                                       cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                        <  4:330094/149‑1 (MQ=255)
                                                                                                       cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGACGCGCCACGTCACTGTTGCTGttt                                        <  3:307006/149‑1 (MQ=255)
                                                                                                        gATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                       >  4:291821/1‑149 (MQ=255)
                                                                                                             ccGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                                  <  1:413945/149‑1 (MQ=255)
                                                                                                               gcgTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                <  4:30695/149‑1 (MQ=255)
                                                                                                               gcgTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                <  3:638403/149‑1 (MQ=255)
                                                                                                                  tCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                             >  3:397716/1‑149 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                         <  1:652358/135‑1 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                         >  2:652358/1‑135 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                           >  2:354841/1‑149 (MQ=255)
                                                                                                                           cTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                    <  2:295289/149‑1 (MQ=255)
                                                                                                                           cTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                    >  4:672388/1‑149 (MQ=255)
                                                                                                                              ccAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                 >  3:417885/1‑149 (MQ=255)
                                                                                                                                aaaaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg               >  3:563352/1‑149 (MQ=255)
                                                                                                                                      aaaTGCCCA‑GCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcgcg         >  2:27696/3‑149 (MQ=255)
                                                                                                                                        aTGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc       <  3:672388/149‑1 (MQ=255)
                                                                                                                                         tGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc      >  2:309478/1‑149 (MQ=255)
                                                                                                                                          gCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     >  3:24592/1‑149 (MQ=255)
                                                                                                                                           cccATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAg    <  4:77113/149‑1 (MQ=255)
                                                                                                                                            ccATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                                 <  3:380699/87‑1 (MQ=255)
                                                                                                                                            ccATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                                 >  4:380699/1‑87 (MQ=255)
                                                                                                                                            ccATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAgg   >  1:121514/1‑149 (MQ=255)
                                                                                                                                             cATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGtt                                               >  4:304167/1‑104 (MQ=255)
                                                                                                                                             cATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGtt                                               <  3:304167/104‑1 (MQ=255)
                                                                                                                                             cATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc  <  2:131296/149‑1 (MQ=255)
                                                                                                                                              aTGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     <  4:350546/145‑1 (MQ=255)
                                                                                                                                              aTGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     >  3:350546/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC  >  NC_000913/1285118‑1285407

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