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

CAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC  >  NC_000913/1285118‑1285407
                                                                                                                                                    |                                                                                                                                              
cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                                                                                               >  4:638403/1‑149 (MQ=255)
cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                                                                                               <  2:48124/149‑1 (MQ=255)
 aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTg                                                                                                                                              >  2:413945/1‑149 (MQ=255)
        cGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtcc                                                                                                                                       <  2:572746/149‑3 (MQ=255)
                  tCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGa                                                                                                                             <  1:195953/149‑1 (MQ=255)
                     aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGAttt                                                                                                                          >  3:330094/1‑149 (MQ=255)
                       gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                        <  4:490877/149‑1 (MQ=255)
                        aTCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                       >  1:295289/1‑149 (MQ=255)
                           gAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                                                                                                         >  4:261798/1‑128 (MQ=255)
                           gAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                                                                                                         <  3:261798/128‑1 (MQ=255)
                            aaCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgacg                                                                                                                   >  1:131296/1‑149 (MQ=255)
                               tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATg                                                                                                                >  3:77113/1‑149 (MQ=255)
                                      tGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGc                                                                                                         <  2:621514/149‑1 (MQ=255)
                                         aTCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCg                                                                                                      >  3:43890/1‑149 (MQ=255)
                                              gTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGcc                                                                                                 <  3:551311/149‑1 (MQ=255)
                                                 ttcTTAAATACCCGGATGCTGCCTTAAGGCAGCATAAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAg                                                                                              >  2:372599/3‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  1:271784/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  3:30695/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  2:481984/1‑149 (MQ=255)
                                                    ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATATAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                           >  1:177511/1‑149 (MQ=255)
                                                      gAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCa                                                                                         <  1:645414/149‑1 (MQ=255)
                                                             cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGc                                                                                  <  1:31007/149‑1 (MQ=255)
                                                             cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGc                                                                                  <  4:24871/149‑1 (MQ=255)
                                                              ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                 >  4:55812/1‑149 (MQ=255)
                                                              ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                 >  4:262868/1‑149 (MQ=255)
                                                                aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGTGGATGCCCATGCGCTGGGCATTTTCCTGCGCGAt                                                                                                                            <  3:460384/104‑1 (MQ=255)
                                                                aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGAt                                                                                                                            >  4:460384/1‑104 (MQ=255)
                                                                     gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                          <  4:43890/149‑1 (MQ=255)
                                                                       cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCg                                                                        >  2:500580/1‑149 (MQ=255)
                                                                            ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGt                                                                   <  2:511315/149‑1 (MQ=255)
                                                                                  aTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgc                                                             >  2:426160/1‑149 (MQ=255)
                                                                                   tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc                                                                                     >  1:199500/1‑124 (MQ=255)
                                                                                   tCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc                                                                                     <  2:199500/124‑1 (MQ=255)
                                                                                       caaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCAc                                                        >  3:204349/1‑149 (MQ=255)
                                                                                        aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGACGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACg                                                       >  4:419086/1‑149 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                               <  4:9985/124‑1 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa                                                                               >  3:9985/1‑124 (MQ=255)
                                                                                         agagATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGt                                                      >  4:273487/1‑149 (MQ=255)
                                                                                                    aGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTg                                           >  4:419166/1‑149 (MQ=255)
                                                                                                       cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                        <  4:330094/149‑1 (MQ=255)
                                                                                                       cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGACGCGCCACGTCACTGTTGCTGttt                                        <  3:307006/149‑1 (MQ=255)
                                                                                                        gATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                       >  4:291821/1‑149 (MQ=255)
                                                                                                             ccGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                                  <  1:413945/149‑1 (MQ=255)
                                                                                                               gcgTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                <  4:30695/149‑1 (MQ=255)
                                                                                                               gcgTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                <  3:638403/149‑1 (MQ=255)
                                                                                                                  tCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                             >  3:397716/1‑149 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                         <  1:652358/135‑1 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGtt                                         >  2:652358/1‑135 (MQ=255)
                                                                                                                    gAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                           >  2:354841/1‑149 (MQ=255)
                                                                                                                           cTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                    <  2:295289/149‑1 (MQ=255)
                                                                                                                           cTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                    >  4:672388/1‑149 (MQ=255)
                                                                                                                              ccAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                 >  3:417885/1‑149 (MQ=255)
                                                                                                                                aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg               >  3:563352/1‑149 (MQ=255)
                                                                                                                                      aaaTGCCCAGCGCATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcgcg         >  2:27696/3‑149 (MQ=255)
                                                                                                                                        aTGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc       <  3:672388/149‑1 (MQ=255)
                                                                                                                                         tGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc      >  2:309478/1‑149 (MQ=255)
                                                                                                                                          gCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     >  3:24592/1‑149 (MQ=255)
                                                                                                                                           cccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAg    <  4:77113/149‑1 (MQ=255)
                                                                                                                                            ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                                 <  3:380699/87‑1 (MQ=255)
                                                                                                                                            ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTgg                                                                 >  4:380699/1‑87 (MQ=255)
                                                                                                                                            ccATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAgg   >  1:121514/1‑149 (MQ=255)
                                                                                                                                             cATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGtt                                               >  4:304167/1‑104 (MQ=255)
                                                                                                                                             cATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGtt                                               <  3:304167/104‑1 (MQ=255)
                                                                                                                                             cATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc  <  2:131296/149‑1 (MQ=255)
                                                                                                                                              aTGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     <  4:350546/145‑1 (MQ=255)
                                                                                                                                              aTGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa     >  3:350546/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC  >  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: