Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,365 G→T 5.6% E81* (GAA→TAA)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,274,3650GT 5.6% 269.8 / 3.3 83E81* (GAA→TAA) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (42/36);  new base T (3/2);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

ACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAATATCGGTGGTG  >  NC_000913/4274219‑4274512
                                                                                                                                                  |                                                                                                                                                   
aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                                                                                                                   <  4:200135/149‑1 (MQ=255)
 cATTACGCTGTCTACTTCCGAATCCTGGCGTGATAAAGAGACCTGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                                                                                                                  <  6:61109/149‑1 (MQ=255)
  aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                                                                                                                 >  1:321416/1‑149 (MQ=255)
    tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTTAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGt                                                                                                                                                >  6:164994/1‑148 (MQ=255)
     aCGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTGGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTc                                                                                                                                               <  6:301541/148‑1 (MQ=255)
         tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTg                                                                                                                                          >  7:340831/1‑149 (MQ=255)
           gCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                                                                                                                         <  6:248228/148‑1 (MQ=255)
             tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTa                                                                                                                                      >  2:221602/1‑149 (MQ=255)
                ttCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCAGCTGCGTAccc                                                                                                                                   >  5:70176/1‑149 (MQ=255)
                ttCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                                                                                                                   >  8:46323/1‑149 (MQ=255)
                ttCCGAATCCTGGCGTGATAAAGCGACCGGAGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGCCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                                                                                                                   >  6:102975/1‑149 (MQ=255)
                 tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCg                                                                                                                                  >  7:181498/1‑149 (MQ=255)
                  ccGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCg                                                                                                                                  >  5:272368/1‑148 (MQ=255)
                     aaTCCTGGCGTGATAAAGCGACCTGCGATATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGTCAGAAGTGTCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTaaa                                                                                                                              <  4:317609/149‑1 (MQ=255)
                      aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt                                                                                                                             <  2:321194/149‑1 (MQ=255)
                         cTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCAGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCATGTTTTTATCGAAGGTCAGCTGCGTACCCGTAAATgg                                                                                                                           <  4:62275/148‑1 (MQ=255)
                          tGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGCGCGAATATCTGCGTAAAGGTTCTCAGGTGTATATCGAAGGTCAGCTGCGTCCCCGTAAATGGAc                                                                                                                         >  8:329265/1‑149 (MQ=255)
                          tGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAc                                                                                                                         >  4:1395/1‑149 (MQ=255)
                              gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     <  3:119497/149‑1 (MQ=255)
                               tGATAAAGCGACCGGCGAGATGAACGAACAGCCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATc                                                                                                                    <  4:218082/149‑1 (MQ=255)
                               tGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     >  5:151866/1‑148 (MQ=255)
                                 aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     <  6:5409/146‑1 (MQ=255)
                                 aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     >  5:10197/1‑146 (MQ=255)
                                 aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     >  5:5409/1‑146 (MQ=255)
                                 aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGGTTATATCTAAGGTCAGCTGCGTACCCGTAAATGGACCGAt                                                                                                                     <  6:10197/146‑1 (MQ=255)
                                  tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt                                                                                                                 >  6:45099/1‑149 (MQ=255)
                                      gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                                                                                                             >  8:317332/1‑149 (MQ=255)
                                        gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGt                                                                                                            >  8:122828/1‑148 (MQ=255)
                                          ccGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAg                                                                                                         <  8:181498/149‑1 (MQ=255)
                                             gCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGAt                                                                                                      >  1:6575/1‑149 (MQ=255)
                                                   tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTa                                                                                                 >  7:300345/1‑148 (MQ=255)
                                                     aaaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTaca                                                                                               <  4:218854/148‑1 (MQ=255)
                                                            aGACTGAATTTCAACGCTGTGTTCTGTTCGTCAAACTGGCAGAAGTTGCTAGCGAATATCTGAGTAAAGGTTCTCTGGTTTATATCGAAGGTCAGATGCGTACCCGTAAATGTACCGATCAATCCGGTCAGGATCGCTACACCACAGaa                                                                                       <  1:181156/149‑1 (MQ=255)
                                                                tGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTa                                                                                                 <  3:123772/135‑1 (MQ=255)
                                                                tGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTa                                                                                                 >  4:123772/1‑135 (MQ=255)
                                                                tGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCg                                                                                   >  7:189547/1‑149 (MQ=255)
                                                                tGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCg                                                                                   >  1:56356/1‑149 (MQ=255)
                                                                 gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgt                                                                                  >  6:169193/1‑149 (MQ=255)
                                                                  aaTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgtg                                                                                 <  7:89901/149‑1 (MQ=255)
                                                                   aTGGCACCGCGTTGTGCTGTTCGGCAAACTTGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCAGTAAATGGACCGATCTATCCGGTCAGGATCGCTACACCACAGAAGTCgtgg                                                                                <  5:286179/149‑1 (MQ=255)
                                                                   aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAg                                                                                         >  7:348526/1‑140 (MQ=255)
                                                                   aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAg                                                                                         <  8:348526/140‑1 (MQ=255)
                                                                   aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgtgg                                                                                <  8:271299/149‑1 (MQ=255)
                                                                    tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgtggt                                                                               <  2:294712/149‑1 (MQ=255)
                                                                        aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAAc                                                                           <  7:122828/149‑1 (MQ=255)
                                                                           gcggcgtgCTGTTCGGCCAACTGGCCGAAGTGGCGAGCGAATCGCTCCGTAAAGGTTCTCAGGTGGATATCGAAGGTCAGCTGCGTACCCGTCCTTGGACCGCTCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGtt                                                                        <  4:168704/144‑1 (MQ=255)
                                                                            cgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTg                                                                       <  7:46323/149‑1 (MQ=255)
                                                                                   gtGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCAc                                                                >  5:9795/2‑149 (MQ=255)
                                                                                     gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCa                                                              >  3:244596/1‑149 (MQ=255)
                                                                                       tctgCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTATATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATAAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCAt                                                             <  5:65356/145‑1 (MQ=255)
                                                                                         ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGc                                                           <  3:74722/148‑1 (MQ=255)
                                                                                             aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGa                                                        >  1:311972/1‑147 (MQ=255)
                                                                                             aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGa                                                        <  2:311972/147‑1 (MQ=255)
                                                                                             aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATg                                                      <  5:164994/149‑1 (MQ=255)
                                                                                              aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACATTGGCGGCACCATGCAGATGc                                                     <  2:81361/149‑1 (MQ=255)
                                                                                              aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGAGCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGc                                                     <  4:141628/149‑1 (MQ=255)
                                                                                                     gAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGggtgg                                              >  8:145856/1‑149 (MQ=255)
                                                                                                      aaGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGCCCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGggtggt                                             >  7:167418/1‑149 (MQ=255)
                                                                                                       aGTGGCGAGCGAATATCTGCCTACAGGATCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAACTGGACCTATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTTAACGTTGGCGGCACCATGCAGATGCTGGGTGgtc                                            <  4:237729/149‑1 (MQ=255)
                                                                                                             gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAgg                                      >  4:74815/1‑149 (MQ=255)
                                                                                                                 gaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGGGAACGTTGGCGGCACCATGCAGATGCTGGGTGGGCGTCAGggtgg                                  >  3:92729/1‑149 (MQ=255)
                                                                                                                   atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCTTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGgtggtg                                <  8:18714/149‑1 (MQ=255)
                                                                                                                      tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTgg                                                  <  3:89467/128‑1 (MQ=255)
                                                                                                                      tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTgg                                                  >  4:89467/1‑128 (MQ=255)
                                                                                                                        tGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgtgg                                                                                <  6:200660/96‑1 (MQ=255)
                                                                                                                        tGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCgtgg                                                                                >  5:200660/1‑96 (MQ=255)
                                                                                                                            tAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCgg                       >  1:163689/1‑149 (MQ=255)
                                                                                                                            tAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCgg                       >  6:14897/1‑149 (MQ=255)
                                                                                                                            tAAAGGTTCTCAGGTTTATATCGAAGGTCAGCGGCGTACCCGTAAATGGACCGCTCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTAGGCGGCACCATGCAGATGCGGGGTGGTCGTCAGGGTGGTGGCGCTCCgg                       >  5:273239/1‑149 (MQ=255)
                                                                                                                              aaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGc                      <  2:56356/148‑1 (MQ=255)
                                                                                                                              aaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGc                      <  8:189547/148‑1 (MQ=255)
                                                                                                                                ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACCCAACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGAGGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAgg                   >  6:288564/1‑149 (MQ=255)
                                                                                                                                 gTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGt                  >  5:308467/1‑149 (MQ=255)
                                                                                                                                   tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTgg                >  4:52222/1‑149 (MQ=255)
                                                                                                                                    ctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCCCCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTgg                >  1:59955/1‑148 (MQ=255)
                                                                                                                                    ctcAGGGTTATATCGAAGGTCAGCTGCGTACCCGTAACTGGACCGATCAATCCGGTCAGGATCGCGACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTgg                <  5:264232/148‑1 (MQ=255)
                                                                                                                                     tcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGc               >  3:102329/1‑148 (MQ=255)
                                                                                                                                      cAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCaa             >  1:257043/1‑149 (MQ=255)
                                                                                                                                       aGGTTTTTATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCTGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCGCCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCaa             <  8:340831/148‑1 (MQ=255)
                                                                                                                                       aGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCaa             <  6:272368/148‑1 (MQ=255)
                                                                                                                                         gTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAatat          >  7:75316/1‑149 (MQ=255)
                                                                                                                                              tatCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAATATCggtg     >  5:40148/1‑149 (MQ=255)
                                                                                                                                                tCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAATATCggtggt   >  6:289608/1‑149 (MQ=255)
                                                                                                                                                  gAAGGTCAGCTGCGTACCCGTATATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAATATCggtggtg  <  2:310805/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAGTCGTGGTGAACGTTGGCGGCACCATGCAGATGCTGGGTGGTCGTCAGGGTGGTGGCGCTCCGGCAGGTGGCAATATCGGTGGTG  >  NC_000913/4274219‑4274512

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