Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,951 G→C 5.3% intergenic (+95/‑94) yihI → / → hemN activator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,9510GC 5.3% 267.1 / 5.1 75intergenic (+95/‑94)yihI/hemNactivator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

AGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATA  >  NC_000913/4051804‑4052093
                                                                                                                                                   |                                                                                                                                              
agGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTcgcg                                                                                                                                                >  4:787656/1‑148 (MQ=255)
    agaagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCgcgg                                                                                                                                               >  3:439788/1‑145 (MQ=255)
    agaagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCgcgg                                                                                                                                               <  4:439788/145‑1 (MQ=255)
       agaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaa                                                                                                                                         <  3:787656/148‑1 (MQ=255)
        gaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaaa                                                                                                                                        <  1:582662/148‑1 (MQ=255)
                   agaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTgg                                                                                                                                   <  2:229328/142‑1 (MQ=255)
                   agaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTgg                                                                                                                                   >  1:229328/1‑142 (MQ=255)
                     aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaaa                                                                                                                                        >  4:872293/1‑135 (MQ=255)
                     aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaaa                                                                                                                                        <  3:872293/135‑1 (MQ=255)
                           tgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                        >  3:125999/1‑145 (MQ=255)
                           tgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                        <  4:125999/145‑1 (MQ=255)
                                cGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCgcggcg                                                                                                                                             <  4:904837/119‑1 (MQ=255)
                                cGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCgcggcg                                                                                                                                             >  3:904837/1‑119 (MQ=255)
                                         cGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACg                                                                                                       <  3:443791/148‑1 (MQ=255)
                                                aCTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTcgcg                                                                                                                                                <  4:1058533/100‑1 (MQ=255)
                                                aCTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTcgcg                                                                                                                                                >  3:1058533/1‑100 (MQ=255)
                                                aCTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTcgcg                                                                                                                                                <  4:606790/100‑1 (MQ=255)
                                                aCTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTcgcg                                                                                                                                                >  3:606790/1‑100 (MQ=255)
                                                   aaCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATGTGGTGTGGTTATTCGTTAAACTACAGCGGTTATCGCGGCGACAAAAGTGGCTGCGAAACCGGCTGATGCCCCCAAACGGGCATccgcc                                                                                             >  2:377516/1‑148 (MQ=255)
                                                   aaCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgcc                                                                                             >  3:380419/1‑148 (MQ=255)
                                                   aaCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgcc                                                                                             >  1:778431/1‑148 (MQ=255)
                                                    aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                                     <  2:46104/123‑1 (MQ=255)
                                                    aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                                     >  1:46104/1‑123 (MQ=255)
                                                    aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACgg                                                                                                      >  1:702835/1‑138 (MQ=255)
                                                    aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACgg                                                                                                      <  2:702835/138‑1 (MQ=255)
                                                    aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgccg                                                                                            >  1:47115/1‑148 (MQ=255)
                                                                   gggCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccgc                                                                             <  2:778431/148‑1 (MQ=255)
                                                                            cGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                                     <  1:985073/99‑1 (MQ=255)
                                                                            cGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                                     >  2:985073/1‑99 (MQ=255)
                                                                               tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                                 <  2:47115/148‑1 (MQ=255)
                                                                                   ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                             >  4:361961/1‑148 (MQ=255)
                                                                                      aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACgg                                                                                                      <  3:668723/104‑1 (MQ=255)
                                                                                      aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACgg                                                                                                      >  4:668723/1‑104 (MQ=255)
                                                                                      aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                          <  4:380419/148‑1 (MQ=255)
                                                                                         aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTgg                                                                                                                                   <  1:1060330/72‑1 (MQ=255)
                                                                                         aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTgg                                                                                                                                   >  2:1060330/1‑72 (MQ=255)
                                                                                              tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                        <  4:20694/78‑1 (MQ=255)
                                                                                              tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                        >  3:20694/1‑78 (MQ=255)
                                                                                              tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                                 >  3:254954/1‑85 (MQ=255)
                                                                                              tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                                 <  4:254954/85‑1 (MQ=255)
                                                                                                         gTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                               >  4:567550/1‑92 (MQ=255)
                                                                                                         gTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                               <  3:567550/92‑1 (MQ=255)
                                                                                                          ttATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaaa                                                                                                                                        <  2:1073025/50‑1 (MQ=255)
                                                                                                          ttATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACaaa                                                                                                                                        >  1:1073025/1‑50 (MQ=255)
                                                                                                          ttATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCaa                                                                                                                           <  2:566392/63‑1 (MQ=255)
                                                                                                          ttATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCaa                                                                                                                           >  1:566392/1‑63 (MQ=255)
                                                                                                           tATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCa                                     >  2:647127/1‑148 (MQ=255)
                                                                                                            aTCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                    >  4:82130/1‑148 (MQ=255)
                                                                                                                   gtgGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACt                             >  3:73920/1‑148 (MQ=255)
                                                                                                                    tgGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccg                                                                              >  2:287020/1‑98 (MQ=255)
                                                                                                                    tgGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccg                                                                              <  1:287020/98‑1 (MQ=255)
                                                                                                                    tgGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTg                            <  1:377516/148‑1 (MQ=255)
                                                                                                                     ggTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                                 >  2:101395/1‑110 (MQ=255)
                                                                                                                     ggTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                                 <  1:101395/110‑1 (MQ=255)
                                                                                                                     ggTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                                 >  2:172581/1‑126 (MQ=255)
                                                                                                                     ggTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                                 <  1:172581/126‑1 (MQ=255)
                                                                                                                      gTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                               >  4:439544/1‑127 (MQ=255)
                                                                                                                      gTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                               <  3:439544/127‑1 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                          >  1:456481/1‑115 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                          <  2:456481/115‑1 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                          >  3:836351/1‑115 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                          <  4:836351/115‑1 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                    >  4:556035/1‑137 (MQ=255)
                                                                                                                       ttATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                    <  3:556035/137‑1 (MQ=255)
                                                                                                                           tCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGa                         >  2:1061085/1‑144 (MQ=255)
                                                                                                                           tCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGa                         <  1:1061085/144‑1 (MQ=255)
                                                                                                                           tCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTg                     <  1:647127/148‑1 (MQ=255)
                                                                                                                              ttAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGCCTGGGATCTGGcc                  <  3:82130/148‑1 (MQ=255)
                                                                                                                                 aaCTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg               <  4:73920/148‑1 (MQ=255)
                                                                                                                                  aCTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc                 >  3:884269/1‑145 (MQ=255)
                                                                                                                                  aCTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc                 <  4:884269/145‑1 (MQ=255)
                                                                                                                                  aCTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa              >  1:175957/1‑148 (MQ=255)
                                                                                                                                     aCAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACg                                                                                                       >  3:982352/1‑56 (MQ=255)
                                                                                                                                     aCAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGAACCGAAACg                                                                                                       <  4:982352/56‑1 (MQ=255)
                                                                                                                                         cGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGa                               >  1:863404/1‑124 (MQ=255)
                                                                                                                                         cGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGa                               <  2:863404/124‑1 (MQ=255)
                                                                                                                                         cGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTg                            <  1:796088/127‑1 (MQ=255)
                                                                                                                                         cGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTg                            >  2:796088/1‑127 (MQ=255)
                                                                                                                                         cGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGaa       <  3:361961/148‑1 (MQ=255)
                                                                                                                                          ttttGTCGCCGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                                      >  3:887852/3‑100 (MQ=255)
                                                                                                                                          ttttGTCGCCGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                                      <  4:887852/98‑1 (MQ=255)
                                                                                                                                           tttGTCGCCGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccg                                                                              >  1:560758/2‑75 (MQ=255)
                                                                                                                                           tttGTCGCCGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccg                                                                              <  2:560758/74‑1 (MQ=255)
                                                                                                                                              gTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAAtata  <  2:175957/148‑1 (MQ=255)
                                                                                                                                                cgcgGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGa                               <  1:343679/117‑1 (MQ=255)
                                                                                                                                                cgcgGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGa                               >  2:343679/1‑117 (MQ=255)
                                                                                                                                                  cggcgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCGAATCGACTGGGATCTGGCCCTGATCCAGAAAtata  <  1:536638/144‑1 (MQ=255)
                                                                                                                                                  cggcgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAAtata  >  2:536638/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATA  >  NC_000913/4051804‑4052093

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