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

ATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGC  >  NC_000913/813423‑813718
                                                                                                                                                    |                                                                                                                                                    
atGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑g                                                                                                                                                     <  2:24294/149‑1 (MQ=255)
    tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCa                                                                                                                                                 <  3:183242/149‑1 (MQ=255)
     gatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAg                                                                                                                                                >  1:23570/1‑149 (MQ=255)
      atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGc                                                                                                                                               <  4:862200/149‑1 (MQ=255)
      atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGc                                                                                                                                               <  4:687769/149‑1 (MQ=255)
         aaCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCag                                                                                                                                            >  1:139406/1‑149 (MQ=255)
          aCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCaga                                                                                                                                           <  2:122457/149‑1 (MQ=255)
           cTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaa                                                                                                                                          <  4:137647/149‑2 (MQ=255)
             gTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGc                                                                                                                                        >  1:420594/1‑149 (MQ=255)
              tttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaaggt                                                                                                                                       <  2:486391/149‑5 (MQ=255)
               ttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAACTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaaggtt                                                                                                                                      <  4:81698/149‑6 (MQ=255)
                  tttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTc                                                                                                                                   <  4:341301/149‑1 (MQ=255)
                     aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGAc                                                                                                                                <  1:742536/149‑1 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCATCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACg                                                                                                                               >  3:670570/1‑149 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaaggtt                                                                                                                                      >  1:599961/1‑137 (MQ=255)
                      tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaaggtt                                                                                                                                      <  2:599961/142‑6 (MQ=255)
                       ccAGTATAATTTGTTGGCATAATTGAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGt                                                                                                                              <  2:508604/149‑1 (MQ=255)
                       ccAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGt                                                                                                                              >  3:538016/1‑149 (MQ=255)
                          gTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGtctc                                                                                                                           <  3:331697/149‑1 (MQ=255)
                             taATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGaa                                                                                                                        <  1:850972/149‑1 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      <  1:106931/149‑1 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      >  1:642521/1‑149 (MQ=255)
                               aTTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAaga                                                                                                                      <  2:478165/149‑1 (MQ=255)
                                 ttgttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAgg                                                                                                                    <  3:522027/149‑1 (MQ=255)
                                  tgttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                   <  2:355472/149‑1 (MQ=255)
                                         aTAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGaa                                                                                                            >  4:446604/1‑149 (MQ=255)
                                           aaTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAAGTAGCGACTCATGAGTAAACCGTTCAAACTAAAATACGCTTATAAACATTTTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGATGGaaaa                                                                                                          <  2:269606/149‑3 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         >  2:658227/1‑149 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         >  1:713804/1‑149 (MQ=255)
                                            aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                         <  4:72738/149‑1 (MQ=255)
                                             ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                        >  1:772384/1‑149 (MQ=255)
                                             ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGATTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                        >  2:885187/1‑149 (MQ=255)
                                             ttAAGTACCACGAGTAAAATTACATACCTGCCCGCCCACCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAACTTCTGGC‑GCTCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGCCTGGAAGATg                                                                                                        <  4:854844/149‑1 (MQ=255)
                                                 gTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCt                                                                                                    >  1:633877/1‑149 (MQ=255)
                                                  tACGACGAGTAAAATTACATACCTGTCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTg                                                                                                   <  2:248401/149‑1 (MQ=255)
                                                    cgacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  2:888886/149‑1 (MQ=255)
                                                     gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTg                                                                                                               >  1:505290/1‑134 (MQ=255)
                                                     gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTg                                                                                                               <  2:505290/134‑1 (MQ=255)
                                                           tAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAg                                                                                          >  1:414656/1‑149 (MQ=255)
                                                             aaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGAc                                                                                        >  4:434303/1‑149 (MQ=255)
                                                             aaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGAc                                                                                        >  4:93615/1‑149 (MQ=255)
                                                                ttACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAACCCGTTCAAACTGAATTCCGCCTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGtt                                                                                     >  2:216206/1‑149 (MQ=255)
                                                                  aCATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTAc                                                                                   <  2:23570/149‑1 (MQ=255)
                                                                    aTACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 >  4:70030/1‑133 (MQ=255)
                                                                    aTACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  3:70030/133‑1 (MQ=255)
                                                                       ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGc                                                                              <  4:706851/149‑1 (MQ=255)
                                                                         tGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGt                                                                            >  3:288643/1‑149 (MQ=255)
                                                                           cccgcccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGa                                                                          >  4:589353/1‑149 (MQ=255)
                                                                              gcccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTg                                                                       >  2:659098/1‑149 (MQ=255)
                                                                                 cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 >  1:103227/1‑120 (MQ=255)
                                                                                 cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgc                                                                                                 <  2:103227/120‑1 (MQ=255)
                                                                                  aaCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                                   >  2:23746/1‑149 (MQ=255)
                                                                                       cTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                              >  2:166069/1‑149 (MQ=255)
                                                                                          cAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                           <  2:605214/149‑1 (MQ=255)
                                                                                            ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                                <  3:266314/142‑1 (MQ=255)
                                                                                            ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                                >  4:266314/1‑142 (MQ=255)
                                                                                              tAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCtt                                                       <  2:597095/149‑1 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     >  1:469840/1‑149 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     <  2:469840/149‑1 (MQ=255)
                                                                                                gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                     <  3:450866/149‑1 (MQ=255)
                                                                                                 cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAc                                                    <  1:498062/149‑1 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaag                                                                                                                                         >  3:197525/1‑58 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaag                                                                                                                                         <  4:197525/60‑3 (MQ=37)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTTGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                                >  3:402819/1‑149 (MQ=255)
                                                                                                     tCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                                >  1:714817/1‑149 (MQ=255)
                                                                                                         gAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCa                                            >  2:667630/1‑149 (MQ=255)
                                                                                                         gAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCa                                            >  1:893679/1‑149 (MQ=255)
                                                                                                           gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAAt                                          <  2:642521/149‑1 (MQ=255)
                                                                                                           gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAAt                                          >  1:588765/1‑149 (MQ=255)
                                                                                                             aaaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGt                                        <  4:538016/149‑1 (MQ=255)
                                                                                                              aaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTc                                       >  2:811274/1‑149 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCaa                                           <  1:315938/144‑1 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCaa                                           >  2:315938/1‑144 (MQ=255)
                                                                                                               aCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCa                                      <  3:538396/149‑1 (MQ=255)
                                                                                                                 cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAtt                                    <  4:332308/149‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTc                                                                                                                                   >  3:637461/1‑52 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTc                                                                                                                                   <  4:637461/52‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGc                                  <  4:433915/149‑1 (MQ=255)
                                                                                                                   ttCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGc                                  >  2:20606/1‑149 (MQ=255)
                                                                                                                     cAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTg                                <  3:446604/149‑1 (MQ=255)
                                                                                                                       aaCTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAc                              >  2:527604/1‑149 (MQ=255)
                                                                                                                       aaCTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAc                              <  1:21108/149‑1 (MQ=255)
                                                                                                                         cTGAATTCCGCTTTTAAACCTTCGGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGGTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCt                            >  2:692575/1‑149 (MQ=255)
                                                                                                                           gAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTc                          >  4:168413/1‑149 (MQ=255)
                                                                                                                                 cGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGcc                    <  3:205254/149‑1 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  4:234837/1‑149 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  4:576051/1‑149 (MQ=255)
                                                                                                                                  gCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                   >  2:670091/1‑149 (MQ=255)
                                                                                                                                         aaCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg            <  3:168413/149‑1 (MQ=255)
                                                                                                                                          aCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATgg           <  3:303680/149‑1 (MQ=255)
                                                                                                                                                tGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTg     >  1:460547/1‑149 (MQ=255)
                                                                                                                                                   c‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTgcgc  <  2:772384/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGC  >  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: ‑