Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 815,667 C→A 6.0% intergenic (+120/+72) uvrB → / ← ybhK exision nuclease of nucleotide excision repair, DNA damage recognition component/putative CofD superfamily transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,6670CA 6.0% 344.7 / 6.8 98intergenic (+120/+72)uvrB/ybhKexision nuclease of nucleotide excision repair, DNA damage recognition component/putative CofD superfamily transferase
Reads supporting (aligned to +/- strand):  ref base C (51/41);  new base A (3/3);  total (54/44)
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

AGCTGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGG  >  NC_000913/815525‑815815
                                                                                                                                              |                                                                                                                                                    
aGCTGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGTTGTTTATGCCGGATGCGGCGTGAAcgcc                                                                                                                                                 >  4:1040498/1‑148 (MQ=25)
   tGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGTTGTTTATGCCGGATGCGGCGTGAAcgcctta                                                                                                                                              <  3:1040498/148‑4 (MQ=25)
   tGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                             >  2:972246/1‑149 (MQ=255)
   tGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                             >  1:690734/1‑149 (MQ=255)
     tttATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                           >  2:180742/1‑149 (MQ=255)
        aTCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                        <  2:630624/149‑1 (MQ=25)
        aTCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                        >  2:972490/1‑149 (MQ=255)
         tCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                        >  3:709124/1‑148 (MQ=255)
           gcgGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                     >  2:777642/1‑149 (MQ=255)
             ggCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaa                                                                                                                                   <  1:933765/149‑1 (MQ=255)
                  cGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac                                                                                                                               >  4:883778/1‑148 (MQ=255)
                  cGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac                                                                                                                               >  4:71488/1‑148 (MQ=255)
                     aaCAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAcgcg                                                                                                                            <  3:779964/148‑1 (MQ=255)
                       cAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCAc                                                                                                                         <  1:952085/149‑1 (MQ=255)
                       cAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTACGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCAc                                                                                                                         <  1:145519/149‑1 (MQ=255)
                         ggATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACtt                                                                                                                       <  2:399182/149‑1 (MQ=255)
                            tAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGt                                                                                                                    <  1:251363/149‑1 (MQ=255)
                            tAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGt                                                                                                                    >  1:473787/1‑149 (MQ=255)
                              gCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCt                                                                                                                  >  1:459775/1‑149 (MQ=255)
                                gaagaaGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCtt                                                                                                                 <  4:459361/148‑1 (MQ=255)
                                gaagaaGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCtt                                                                                                                 <  4:408486/148‑1 (MQ=255)
                                 aagaagACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTc                                                                                                                >  4:628577/1‑148 (MQ=255)
                                  agaagaCTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCaa                                                                                                              <  2:982261/149‑1 (MQ=255)
                                  agaagaCTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCa                                                                                                               <  4:382772/148‑1 (MQ=255)
                                        cTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGa                                                                                                         <  4:876322/148‑1 (MQ=255)
                                             gaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTCCAATAGGAAGCCgg                                                                                                   >  2:84338/1‑149 (MQ=255)
                                                aaaCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAt                                                                                                >  1:258416/1‑149 (MQ=255)
                                                aaaCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAt                                                                                                >  1:365726/1‑149 (MQ=255)
                                                 aaCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtt                                                                                               >  2:581574/1‑149 (MQ=255)
                                                 aaCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAt                                                                                                >  3:193144/1‑148 (MQ=255)
                                                  aCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtt                                                                                               >  3:1049818/1‑148 (MQ=255)
                                                  aCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtt                                                                                               >  4:757748/1‑148 (MQ=255)
                                                   cGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAttt                                                                                              >  3:336820/1‑148 (MQ=255)
                                                    ggAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtttt                                                                                             <  3:71488/148‑1 (MQ=255)
                                                    ggAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtttt                                                                                             >  3:844128/1‑148 (MQ=255)
                                                    ggAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTc                                                                                            >  1:47818/1‑149 (MQ=255)
                                                    ggAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTc                                                                                            >  2:812215/1‑149 (MQ=255)
                                                       aaaCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACttt                                                                                                                      >  1:10745/1‑120 (MQ=255)
                                                       aaaCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACttt                                                                                                                      <  2:10745/120‑1 (MQ=255)
                                                       aaaCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCtt                                                                                         <  2:47818/149‑1 (MQ=255)
                                                           aGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGa                                                                                     <  2:722985/149‑1 (MQ=255)
                                                           aGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGa                                                                                     <  1:545158/149‑1 (MQ=255)
                                                            gCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGa                                                                                     >  4:661737/1‑148 (MQ=255)
                                                                  tGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCg                                                                               >  4:391268/1‑148 (MQ=255)
                                                                  tGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGt                                                                              >  2:732247/1‑149 (MQ=255)
                                                                     gcTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTaa                                                                           >  1:678273/1‑149 (MQ=255)
                                                                     gcTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTa                                                                            <  3:661866/148‑1 (MQ=255)
                                                                      cTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTaa                                                                           <  4:216381/148‑1 (MQ=255)
                                                                       gtagcgcaTCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAAc                                                                          <  4:259263/141‑1 (MQ=255)
                                                                        aCGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAAcc                                                                         >  3:90467/1‑148 (MQ=255)
                                                                        aCGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCAGTTAAcc                                                                         >  3:909283/1‑148 (MQ=255)
                                                                             tATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAg                                                                    <  3:628577/148‑1 (MQ=255)
                                                                                 aGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTg                                                                <  3:76605/148‑1 (MQ=255)
                                                                                 aGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGt                                                               >  1:811757/1‑149 (MQ=255)
                                                                                   gCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTaa                                                             <  1:682261/149‑1 (MQ=255)
                                                                                      tACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACg                                                           >  4:24117/1‑148 (MQ=255)
                                                                                       aCATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGcc                                                                  >  4:492176/1‑140 (MQ=255)
                                                                                       aCATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGcc                                                                  <  3:492176/140‑1 (MQ=255)
                                                                                       aCATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGc                                                          <  4:193144/148‑1 (MQ=255)
                                                                                         aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtt                                                       <  1:45964/149‑1 (MQ=255)
                                                                                         aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtt                                                       <  1:972490/149‑1 (MQ=255)
                                                                                         aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtt                                                       <  2:290283/149‑1 (MQ=255)
                                                                                                ccGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAg                                                 >  4:827155/1‑148 (MQ=255)
                                                                                                  gCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAgcg                                               >  3:862040/1‑148 (MQ=255)
                                                                                                      tatatTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGtt                                           <  3:827155/148‑1 (MQ=255)
                                                                                                      tatatTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTa                                          <  2:258416/149‑1 (MQ=255)
                                                                                                          ttGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      >  1:194521/1‑149 (MQ=255)
                                                                                                          ttGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      >  2:705153/1‑149 (MQ=255)
                                                                                                           tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      <  4:755694/148‑1 (MQ=255)
                                                                                                           tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      >  4:717746/1‑148 (MQ=255)
                                                                                                           tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      >  3:323825/1‑148 (MQ=255)
                                                                                                           tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                      >  4:415169/1‑148 (MQ=255)
                                                                                                            gAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaac                                    >  1:739969/1‑149 (MQ=255)
                                                                                                            gAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaa                                     >  3:116216/1‑148 (MQ=255)
                                                                                                              aTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTagca                                   <  4:449172/148‑1 (MQ=255)
                                                                                                              aTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaacaa                                  >  1:145955/1‑149 (MQ=255)
                                                                                                                ttGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTaa                                                                      >  3:462516/1‑111 (MQ=255)
                                                                                                                ttGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTaa                                                                      <  4:462516/111‑1 (MQ=255)
                                                                                                                        tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtt                                                       >  4:1044900/1‑118 (MQ=255)
                                                                                                                        tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtt                                                       <  3:1044900/118‑1 (MQ=255)
                                                                                                                          tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATg                      <  1:732247/149‑1 (MQ=255)
                                                                                                                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATg                      <  3:661737/148‑1 (MQ=255)
                                                                                                                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGa                     <  1:342284/149‑1 (MQ=255)
                                                                                                                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGa                     <  1:202028/149‑1 (MQ=255)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGAt                 >  4:404445/1‑148 (MQ=255)
                                                                                                                                ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGAt                 >  4:396933/1‑148 (MQ=255)
                                                                                                                                ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGAt                 >  4:410133/1‑148 (MQ=255)
                                                                                                                                ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACAAt                 <  4:482381/148‑1 (MQ=255)
                                                                                                                                   gatgcggCGTGAACGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTa                                                              >  4:574245/6‑100 (MQ=37)
                                                                                                                                   gatgcggCGTGAACGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTa                                                              <  3:574245/95‑1 (MQ=37)
                                                                                                                                     tAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGa            >  4:41539/1‑148 (MQ=255)
                                                                                                                                       aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTa                                                                            >  3:83285/1‑82 (MQ=255)
                                                                                                                                       aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTa                                                                            <  4:83285/82‑1 (MQ=255)
                                                                                                                                       aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAat          >  4:521139/1‑148 (MQ=255)
                                                                                                                                         gCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATc       >  2:576187/1‑149 (MQ=255)
                                                                                                                                           gTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGc     >  1:697323/1‑149 (MQ=255)
                                                                                                                                            ttCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGc     <  3:757748/148‑1 (MQ=255)
                                                                                                                                            ttCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGc     <  4:396612/148‑1 (MQ=255)
                                                                                                                                              cACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTgg  <  2:678273/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AGCTGTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGG  >  NC_000913/815525‑815815

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