Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,768 C→A 9.4% intergenic (+242/+298) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,7680CA 9.4% 235.7 / 8.2 78intergenic (+242/+298)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base A (4/4);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGC  >  NC_000913/1286626‑1286915
                                                                                                                                              |                                                                                                                                                   
tgagtaTATTGTAGGCCCGGTAAACGTGGTGACATCGGGCAAAATTAGCTCGGAGTAATATGTTTAGATTATTTAAGGAAGCGATGGTTGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                <  2:91765/147‑1 (MQ=25)
ggaggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                >  2:167615/1‑148 (MQ=255)
   ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACga                                                                                                                                             <  5:92905/148‑2 (MQ=255)
     atatTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctc                                                                                                                                           >  2:330618/1‑148 (MQ=255)
        ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctct                                                                                                                                        >  4:272509/1‑148 (MQ=255)
         tGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                       >  3:311964/1‑148 (MQ=255)
         tGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACga                                                                                                                                             >  7:25010/1‑141 (MQ=255)
         tGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACga                                                                                                                                             <  8:25010/142‑2 (MQ=255)
         tGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                       >  7:438139/1‑148 (MQ=255)
           tAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                     <  2:8820/148‑1 (MQ=255)
             ggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   <  3:272509/148‑1 (MQ=255)
                 cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGTAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGgcg                                                                                                                                >  2:112053/1‑147 (MQ=255)
                 cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                               >  6:152061/1‑148 (MQ=255)
                 cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                               >  4:177139/1‑148 (MQ=255)
                  ggTAAACGTGGTGCCATCGGGCAAATTTAGCTCGGAGTAACGGGTTTTGATTATTTAAGGAAGCGATGGTGCTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCGCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGgcgca                                                                                                                              >  7:311330/1‑147 (MQ=25)
                  ggTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                              >  4:87695/1‑148 (MQ=255)
                    tAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             >  8:240090/1‑147 (MQ=255)
                    tAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             >  6:72001/1‑147 (MQ=255)
                       aCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGt                                                                                                                         <  3:157953/148‑1 (MQ=255)
                       aCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGt                                                                                                                         <  2:337764/148‑1 (MQ=255)
                          tggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACga                                                                                                                                             <  6:65702/125‑2 (MQ=255)
                          tggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACga                                                                                                                                             >  5:65702/1‑124 (MQ=255)
                               ccATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCtt                                                                                                                 >  2:88326/1‑148 (MQ=255)
                                 aTCGAGCAAAATTAGCTCGGAGTAACAGGGTTTGATTAGTTAAGGAAGCGATGGTGGTGGGGGCAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTgg                                                                                                               <  4:388774/148‑1 (MQ=255)
                                     ggCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCg                                                                                                            >  8:141405/1‑147 (MQ=255)
                                      gCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             <  7:439222/129‑1 (MQ=255)
                                      gCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             >  8:439222/1‑129 (MQ=255)
                                         aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAGGCt                                                                                                       <  2:9846/148‑1 (MQ=255)
                                         aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAGGCt                                                                                                       <  2:2787/148‑1 (MQ=255)
                                               gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                 <  3:87695/148‑1 (MQ=255)
                                               gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                 <  1:330618/148‑1 (MQ=255)
                                                   ggAGTATCAGGTGTTGCTTATTTAAGGCAGCGATGGTGGTGGGGGAAGGAGTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTGCTCATCCtt                                                                                             <  1:264107/148‑1 (MQ=14)
                                                    gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTTGGGGAAGGATTACTAAGAGCTGCGCGCTTCGCCCCTCGGGTCGTTGCATGCGGCAACGCTCTCTCTCTGGCGCTCAAGTCGAACCTTGGTAGAAGCtt                                                                                                      >  5:324369/1‑138 (MQ=255)
                                                    gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                            <  1:167615/148‑1 (MQ=255)
                                                        aaaaGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGAAGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTACCCg                                                                                        >  2:193241/4‑148 (MQ=17)
                                                        aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                >  3:101797/1‑140 (MQ=255)
                                                        aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  4:101797/140‑1 (MQ=255)
                                                        aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCg                                                                                        >  2:268926/1‑148 (MQ=255)
                                                          cAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                        >  7:397852/1‑130 (MQ=255)
                                                          cAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                        <  8:397852/130‑1 (MQ=255)
                                                           aGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctct                                                                                                                                        <  5:197252/97‑1 (MQ=255)
                                                           aGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGATTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctct                                                                                                                                        >  6:197252/1‑97 (MQ=255)
                                                               tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               >  5:105432/1‑134 (MQ=255)
                                                               tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  6:105432/134‑1 (MQ=255)
                                                               tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                >  5:120522/1‑133 (MQ=255)
                                                               tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  6:120522/133‑1 (MQ=255)
                                                                ttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                      <  2:111551/94‑1 (MQ=255)
                                                                ttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                      >  1:111551/1‑94 (MQ=255)
                                                                   attattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGaa                                                                             >  5:364607/1‑148 (MQ=255)
                                                                       tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAatat                                                                          <  5:55263/147‑1 (MQ=255)
                                                                       tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAtt                                                                         <  6:390934/148‑1 (MQ=255)
                                                                        ttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAgc                                                                                                                         >  7:420492/1‑98 (MQ=12)
                                                                        ttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAgc                                                                                                                         <  8:420492/99‑2 (MQ=12)
                                                                          aaGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                       <  4:125364/147‑1 (MQ=40)
                                                                          aaGGAAGCGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAATTCGAACCTTGGTCGAAGCTTCTAATCCTTCCCCGCTTGGGCAGAATATTTg                                                                       <  1:193241/147‑1 (MQ=255)
                                                                             gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTg                                                                   <  6:364607/148‑1 (MQ=40)
                                                                             gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTg                                                                   <  7:155578/148‑1 (MQ=40)
                                                                                 cGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGa                                                               >  4:397121/1‑148 (MQ=38)
                                                                                            ggggAATGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGACGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgaga                                                    <  5:186578/148‑1 (MQ=255)
                                                                                               gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAAtt                                                 >  1:371728/1‑148 (MQ=255)
                                                                                                 aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTcc                                               >  8:297436/1‑148 (MQ=255)
                                                                                                 aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTcc                                               <  1:325147/148‑1 (MQ=255)
                                                                                                        cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGTAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCtt                                        <  3:397121/148‑1 (MQ=255)
                                                                                                        cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCtt                                        <  2:371728/148‑1 (MQ=255)
                                                                                                          cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAAt                                                  >  5:315145/1‑136 (MQ=255)
                                                                                                          cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAAt                                                  <  6:315145/136‑1 (MQ=255)
                                                                                                            gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGa                                    <  7:240090/148‑1 (MQ=255)
                                                                                                             cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaa                                   <  5:152061/148‑1 (MQ=255)
                                                                                                                 gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGt                                >  6:114803/1‑147 (MQ=255)
                                                                                                                 gcgcgcTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTg                               <  5:72001/148‑1 (MQ=255)
                                                                                                                   gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAt                             <  5:114803/148‑1 (MQ=255)
                                                                                                                     gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtgg                           <  4:311964/148‑1 (MQ=255)
                                                                                                                           cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGTAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtgg                     <  1:112053/148‑1 (MQ=255)
                                                                                                                                   agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCt                                         >  2:249210/3‑120 (MQ=255)
                                                                                                                                   agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCt                                         <  1:249210/118‑1 (MQ=255)
                                                                                                                                      gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa          <  7:297436/148‑1 (MQ=255)
                                                                                                                                        tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaa                                   >  4:310525/1‑121 (MQ=255)
                                                                                                                                        tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaa                                   <  3:310525/121‑1 (MQ=255)
                                                                                                                                              cGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTCGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc  >  5:200454/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGC  >  NC_000913/1286626‑1286915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: