Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 813,570:1 +T 7.8% coding (45/2022 nt) uvrB → exision nuclease of nucleotide excision repair, DNA damage recognition component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913813,5701.T 7.8% 359.4 / 26.7 90coding (45/2022 nt)uvrBexision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base . (42/41);  new base T (2/5);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

ATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGC  >  NC_000913/813423‑813718
                                                                                                                                                    |                                                                                                                                                    
atGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCg                                                                                                                                                     <  2:24294/149‑1 (MQ=255)
    tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCa                                                                                                                                                 <  3:183242/149‑1 (MQ=255)
     gatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAg                                                                                                                                                >  1:23570/1‑149 (MQ=255)
      atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGc                                                                                                                                               <  4:862200/149‑1 (MQ=255)
      atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGc                                                                                                                                               <  4:687769/149‑1 (MQ=255)
         aaCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCag                                                                                                                                            >  1:139406/1‑149 (MQ=255)
          aCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCaga                                                                                                                                           <  2:122457/149‑1 (MQ=255)
           cTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaa                                                                                                                                          <  4:137647/149‑2 (MQ=255)
             gTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGc                                                                                                                                        >  1:420594/1‑149 (MQ=255)
              tttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaaggt                                                                                                                                       <  2:486391/149‑5 (MQ=255)
               ttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAACTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaaggtt                                                                                                                                      <  4:81698/149‑6 (MQ=255)
                  tttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTc                                                                                                                                   <  4:341301/149‑1 (MQ=255)
                     aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGAc                                                                                                                                <  1:742536/149‑1 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCATCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACg                                                                                                                               >  3:670570/1‑149 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaaggtt                                                                                                                                      >  1:599961/1‑137 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaaggtt                                                                                                                                      <  2:599961/142‑6 (MQ=255)
                       ccAGTATAATTTGTTGGCATAATTGAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGt                                                                                                                              <  2:508604/149‑1 (MQ=255)
                       ccAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGt                                                                                                                              >  3:538016/1‑149 (MQ=255)
                          gTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtctc                                                                                                                           <  3:331697/149‑1 (MQ=255)
                             taATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGaa                                                                                                                        <  1:850972/149‑1 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      <  1:106931/149‑1 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      >  1:642521/1‑149 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      <  2:478165/149‑1 (MQ=255)
                                 ttgttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAgg                                                                                                                    <  3:522027/149‑1 (MQ=255)
                                  tgttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                   <  2:355472/149‑1 (MQ=255)
                                         aTAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGaa                                                                                                            >  4:446604/1‑149 (MQ=255)
                                           aaTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAAGTAGCGACTCATGAGTAAACCGTTCAAACTAAAATACGCTTATAAACATTTTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGATGGaaaa                                                                                                          <  2:269606/149‑3 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         >  2:658227/1‑149 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         >  1:713804/1‑149 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         <  4:72738/149‑1 (MQ=255)
                                             ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                        >  1:772384/1‑149 (MQ=255)
                                             ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGATTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                        >  2:885187/1‑149 (MQ=255)
                                             ttAAGTACCACGAGTAAAATTACATACCTGCCCGCCCACCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAACTTCTGGCGCTCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGCCTGGAAGATg                                                                                                        <  4:854844/149‑1 (MQ=255)
                                                 gTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCt                                                                                                    >  1:633877/1‑149 (MQ=255)
                                                  tACGACGAGTAAAATTACATACCTGTCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTg                                                                                                   <  2:248401/149‑1 (MQ=255)
                                                    cgacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  2:888886/149‑1 (MQ=255)
                                                     gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTg                                                                                                               >  1:505290/1‑134 (MQ=255)
                                                     gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTg                                                                                                               <  2:505290/134‑1 (MQ=255)
                                                           tAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAg                                                                                          >  1:414656/1‑149 (MQ=255)
                                                             aaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGAc                                                                                        >  4:434303/1‑149 (MQ=255)
                                                             aaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGAc                                                                                        >  4:93615/1‑149 (MQ=255)
                                                                ttACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAACCCGTTCAAACTGAATTCCGCCTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGtt                                                                                     >  2:216206/1‑149 (MQ=255)
                                                                  aCATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTAc                                                                                   <  2:23570/149‑1 (MQ=255)
                                                                    aTACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 >  4:70030/1‑133 (MQ=255)
                                                                    aTACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  3:70030/133‑1 (MQ=255)
                                                                       ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGc                                                                              <  4:706851/149‑1 (MQ=255)
                                                                         tGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGt                                                                            >  3:288643/1‑149 (MQ=255)
                                                                           cccgcccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGa                                                                          >  4:589353/1‑149 (MQ=255)
                                                                              gcccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTg                                                                       >  2:659098/1‑149 (MQ=255)
                                                                                 cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 >  1:103227/1‑120 (MQ=255)
                                                                                 cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  2:103227/120‑1 (MQ=255)
                                                                                  aaCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                                   >  2:23746/1‑149 (MQ=255)
                                                                                       cTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                              >  2:166069/1‑149 (MQ=255)
                                                                                          cAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                           <  2:605214/149‑1 (MQ=255)
                                                                                            ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                                <  3:266314/142‑1 (MQ=255)
                                                                                            ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                                >  4:266314/1‑142 (MQ=255)
                                                                                              tAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCtt                                                       <  2:597095/149‑1 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     >  1:469840/1‑149 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     <  2:469840/149‑1 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     <  3:450866/149‑1 (MQ=255)
                                                                                                 cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAc                                                    <  1:498062/149‑1 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                         >  3:197525/1‑58 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                         <  4:197525/60‑3 (MQ=37)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTTGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                                >  3:402819/1‑149 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                                >  1:714817/1‑149 (MQ=255)
                                                                                                         gAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCa                                            >  2:667630/1‑149 (MQ=255)
                                                                                                         gAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCa                                            >  1:893679/1‑149 (MQ=255)
                                                                                                           gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAAt                                          <  2:642521/149‑1 (MQ=255)
                                                                                                           gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAAt                                          >  1:588765/1‑149 (MQ=255)
                                                                                                             aaaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGt                                        <  4:538016/149‑1 (MQ=255)
                                                                                                              aaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTc                                       >  2:811274/1‑149 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCaa                                           <  1:315938/144‑1 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCaa                                           >  2:315938/1‑144 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCa                                      <  3:538396/149‑1 (MQ=255)
                                                                                                                 cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAtt                                    <  4:332308/149‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTc                                                                                                                                   >  3:637461/1‑52 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTc                                                                                                                                   <  4:637461/52‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGc                                  <  4:433915/149‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGc                                  >  2:20606/1‑149 (MQ=255)
                                                                                                                     cAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTg                                <  3:446604/149‑1 (MQ=255)
                                                                                                                       aaCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAc                              >  2:527604/1‑149 (MQ=255)
                                                                                                                       aaCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAc                              <  1:21108/149‑1 (MQ=255)
                                                                                                                         cTGAATTCCGCTTTTAAACCTTCGGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGGTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCt                            >  2:692575/1‑149 (MQ=255)
                                                                                                                           gAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTc                          >  4:168413/1‑149 (MQ=255)
                                                                                                                                 cGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGcc                    <  3:205254/149‑1 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  4:234837/1‑149 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  4:576051/1‑149 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  2:670091/1‑149 (MQ=255)
                                                                                                                                         aaCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg            <  3:168413/149‑1 (MQ=255)
                                                                                                                                          aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATgg           <  3:303680/149‑1 (MQ=255)
                                                                                                                                                tGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTg     >  1:460547/1‑149 (MQ=255)
                                                                                                                                                   cGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTgcgc  <  2:772384/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGC  >  NC_000913/813423‑813718

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: