Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,273,361 T→A 8.6% S171C (AGC→TGC)  uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,273,3610TA 8.6% 393.4 / 22.3 116S171C (AGC→TGC) uvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base T (63/43);  new base A (4/6);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 117 total aligned reads displayed.

ACTCGGCAAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTCGTGGT  >  NC_000913/4273214‑4273507
                                                                                                                                                   |                                                                                                                                                  
aCTCGGCAAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCtt                                                                                                                                                   <  1:791703/149‑1 (MQ=255)
   cGGCAAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGcc                                                                                                                                                <  1:268011/149‑1 (MQ=255)
     gCAAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAg                                                                                                                                              <  1:374909/149‑1 (MQ=255)
        aGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                           <  1:95091/149‑3 (MQ=255)
         gACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggct                                                                                                                                          <  2:182998/149‑4 (MQ=255)
         gACGCTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggct                                                                                                                                          <  1:380297/149‑4 (MQ=255)
          aCGTTGGGTAAGATCGTCACGAACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGTAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTtc                                                                                                                                          >  3:136435/1‑148 (MQ=255)
           cGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTtct                                                                                                                                         <  3:399973/148‑1 (MQ=255)
            gTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTcc                                                                                                                                       <  2:1060133/149‑1 (MQ=255)
             ttGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTcc                                                                                                                                       >  3:231548/1‑148 (MQ=255)
                   aaGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGtt                                                                                                                                 >  4:551332/1‑148 (MQ=255)
                    aGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCCCAGACTTCGGCATCAATACGAGCACGGATGTAGCGCTGGCTTGCCAGGTTCTCCAGCGtttt                                                                                                                               <  1:635350/149‑1 (MQ=255)
                      aTCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTgg                                                                                                                             <  1:523950/149‑1 (MQ=255)
                      aTCGTCACGCACCTTGAAGCGATCAACCACCACTTCAAGGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTg                                                                                                                              >  3:822728/1‑148 (MQ=255)
                       tCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            <  2:863567/149‑1 (MQ=255)
                         gTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCTTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            <  2:1036875/147‑1 (MQ=255)
                         gTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCTTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            >  1:1036875/1‑147 (MQ=255)
                             cGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGtt                                                                                                                                           <  3:221632/128‑1 (MQ=255)
                             cGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGtt                                                                                                                                           >  4:221632/1‑128 (MQ=255)
                                aCCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGcc                                                                                                                   >  2:988082/1‑149 (MQ=255)
                                   ttGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            <  1:157408/137‑1 (MQ=255)
                                   ttGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            >  2:157408/1‑137 (MQ=255)
                                   ttGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGATTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCttt                                                                                                                >  1:450925/1‑149 (MQ=255)
                                    tGAAGCGATCTACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCCGGTTCTCCAGCGTTTTGGTGTGTTCGCCttt                                                                                                                <  4:905596/148‑1 (MQ=255)
                                    tGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            >  1:544184/1‑136 (MQ=255)
                                    tGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGgt                                                                                                                            <  2:544184/136‑1 (MQ=255)
                                     gAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTg                                                                                                               <  3:551332/148‑1 (MQ=255)
                                       aGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGtt                                                                                                                       >  4:664565/1‑138 (MQ=255)
                                       aGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTgcg                                                                                                             >  3:226295/1‑148 (MQ=255)
                                       aGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggct                                                                                                                                          <  1:399667/119‑4 (MQ=255)
                                       aGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggct                                                                                                                                          >  2:399667/1‑116 (MQ=255)
                                        gCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGct                                                                                                           >  2:173675/1‑149 (MQ=255)
                                        gCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                           >  4:9093/1‑115 (MQ=255)
                                        gCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                           <  3:9093/117‑3 (MQ=255)
                                          gATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGctct                                                                                                         >  1:344732/1‑149 (MQ=255)
                                            tCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTg                                                                                                                              >  1:689552/1‑126 (MQ=255)
                                            tCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTg                                                                                                                              <  2:689552/126‑1 (MQ=255)
                                             cAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                           <  3:472671/112‑3 (MQ=255)
                                             cAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                           >  4:472671/1‑110 (MQ=255)
                                                 caccacTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGa                                                                                                  >  1:524662/1‑149 (MQ=255)
                                                 caccacTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCGCTTTAATg                                                                                                   >  3:970863/1‑148 (MQ=255)
                                                  accacTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGa                                                                                                  >  4:204643/1‑148 (MQ=255)
                                                  accacTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGa                                                                                                  >  3:936922/1‑148 (MQ=255)
                                                   ccacTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGAt                                                                                                 <  4:231548/148‑1 (MQ=255)
                                                      cTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTgg                                                                                              >  3:387817/1‑148 (MQ=255)
                                                             ggTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTa                                                                                      <  1:490281/149‑1 (MQ=255)
                                                              gTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAg                                                                                     <  2:689512/149‑1 (MQ=255)
                                                                  gTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCAt                                                                                  >  3:423249/1‑148 (MQ=255)
                                                                   tttCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATc                                                                                 <  4:106697/148‑1 (MQ=255)
                                                                     tCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGa                                                                              <  2:909805/149‑1 (MQ=255)
                                                                      cTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGAc                                                                             <  1:73675/149‑1 (MQ=255)
                                                                       ttttGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACg                                                                            >  2:410701/1‑149 (MQ=255)
                                                                        tttGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACg                                                                            >  4:424292/1‑148 (MQ=255)
                                                                         ttGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGc                                                                           <  4:254402/148‑1 (MQ=255)
                                                                          tGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCtt                                                                         <  2:163365/149‑1 (MQ=255)
                                                                          tGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCt                                                                          <  3:404339/148‑1 (MQ=255)
                                                                            cAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTg                                                                        >  3:372722/1‑148 (MQ=255)
                                                                                tCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCtt                                                                    >  3:799124/1‑148 (MQ=255)
                                                                                tCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCtt                                                                    >  3:754539/1‑148 (MQ=255)
                                                                                tCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCtt                                                                    >  4:722984/1‑148 (MQ=255)
                                                                                 ccAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTc                                                                   >  3:1018681/1‑148 (MQ=255)
                                                                                   agctttGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTaa                                                                                                     <  4:83576/109‑1 (MQ=255)
                                                                                   aGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTaa                                                                                                     >  3:83576/1‑112 (MQ=255)
                                                                                      tttGGCGGATAGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTGGCATCAGACGCTTGCCTTCCGGCt                                                              >  3:697975/1‑148 (MQ=255)
                                                                                       ttGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTg                                                             <  4:260108/148‑1 (MQ=255)
                                                                                                 ggAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGt                                                  >  2:817049/1‑149 (MQ=255)
                                                                                                  gAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGtt                                                 >  1:582426/1‑149 (MQ=255)
                                                                                                  gAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGt                                                  >  4:522786/1‑148 (MQ=255)
                                                                                                      gATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATc                                              >  4:707656/1‑148 (MQ=255)
                                                                                                        tCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATcca                                            <  3:567775/148‑1 (MQ=255)
                                                                                                         cGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATccacc                                          >  2:403140/1‑149 (MQ=255)
                                                                                                           cAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATccacca                                         >  4:556980/1‑148 (MQ=255)
                                                                                                           cAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGCTTTGGTGTGTTCGCATTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCGTCGACAGCACGTTATCCACCAt                                        >  1:363568/1‑149 (MQ=255)
                                                                                                               cTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTgg                                    <  1:825380/149‑1 (MQ=255)
                                                                                                                ttCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGc                                   >  2:241309/1‑149 (MQ=255)
                                                                                                                   gCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTg                                 <  3:522786/148‑1 (MQ=255)
                                                                                                                   gCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGa                                >  2:5527/1‑149 (MQ=255)
                                                                                                                    ccATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGAGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGa                                >  3:88779/1‑148 (MQ=255)
                                                                                                                     cATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCa                                                     <  3:603758/126‑1 (MQ=255)
                                                                                                                     cATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCa                                                     >  4:603758/1‑126 (MQ=255)
                                                                                                                       tCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTGTCCACCATCt                                      <  3:846551/139‑1 (MQ=255)
                                                                                                                       tCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCt                                      >  4:846551/1‑139 (MQ=255)
                                                                                                                       tCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGt                            <  1:328349/149‑1 (MQ=255)
                                                                                                                        cAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGt                            >  4:523844/1‑148 (MQ=255)
                                                                                                                         aaTACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGttt                          >  2:84573/1‑149 (MQ=255)
                                                                                                                         aaTACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGtt                           <  3:204643/148‑1 (MQ=255)
                                                                                                                                gCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGcgccgc                   >  2:482123/1‑149 (MQ=255)
                                                                                                                                   cGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAg                <  2:524662/149‑1 (MQ=255)
                                                                                                                                      aTGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCgg             >  1:866164/1‑149 (MQ=255)
                                                                                                                                      aTGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCg              <  4:697975/148‑1 (MQ=255)
                                                                                                                                      aTGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCg              <  4:822728/148‑1 (MQ=255)
                                                                                                                                       tGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTg                                                                        >  2:604166/1‑89 (MQ=255)
                                                                                                                                       tGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTg                                                                        <  1:604166/89‑1 (MQ=255)
                                                                                                                                          gaaccTGGCAAGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGAcg         >  2:395216/4‑149 (MQ=255)
                                                                                                                                            cccTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGAc                               <  4:420837/125‑1 (MQ=255)
                                                                                                                                            cccTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGAcgtc       <  1:173675/149‑1 (MQ=255)
                                                                                                                                            cccTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGAcgt        >  3:957967/1‑148 (MQ=255)
                                                                                                                                                ggCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATccacc                                          <  3:394463/110‑1 (MQ=255)
                                                                                                                                                ggCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTCGTg    >  3:372701/1‑148 (MQ=255)
                                                                                                                                                  cTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTCGTGGt  >  3:932686/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACTCGGCAAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTCGTGGT  >  NC_000913/4273214‑4273507

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